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The  Earths  Radiation  Belts 


Interaction  of  the  solar  wind  flow  with  the  earth's  magnetic  field  gives  rise  t  i 
a  cavity  in  the  inter|)lanetarv  medium  known  as  the  earth's  magnetosphere.  Within 
this  eavitv  there  exists  a  limited  region  where  the  motion  >f  energetic  partieles  is 
.  onfined  bv  the  earth's  magnetic  field.  This  region  comprises  the  earth's  radiation 
belts  as  depicted  in  figure  5.  1.  The  radiation  belt  region  contains  electrons, 
protons,  helium,  carbon,  oxvgen,  and  other  ions  with  energies  from  less  than  1  ke\ 
to  hundreds  of  MeV.  Particles  below  200  ke\  energy  c  institute  the  principal 
corpuscular  energy  density  and  form  the  extra-terrestrial  ring  current.  Confine¬ 
ment  (or  trapping)  of  these  particles  results  from  the  dipolar-like  topology  of  the 
geomagnetic  field,  which  is  characterized  bv  magnetic  field  lines  that  converge  at 
high  latitudes  towards  the  poles  resulting  in  a  relative  minimum  magnetic  field 
strength  region  in  the  vicinity  of  the  geomagnetic  equator,  f  igure  5.2  illustrates 
principal  aspects  of  a  charged  particle  trajectory  in  magnetic  mirror  field  geometry. 

The  magnetic  force  (f  v  s  H)  deflects  the  particle  veloci tv  vector  \  so  tiiat  the 

— > 

particle  spirals  around  the  magnetic  field  U.  The  convc rgencc  of  the  magnetic 
lines  of  torce  causes  a  tightening  ol  the  spiral  angle  and  eventually  a  reflection  of  the 
particle  from  the  high  magnetic  field  region  (mirror  point);  for  this  reason  the  earth' 
magnetic  field  is  capable  of  confining  charged  partic  les.  A  detailed  account  of  single 
nartiele  motion  in  magnetic  fields  is  found  in  K  ledrrer.  * 

(Received  for  publication  26  August  1!'82) 

1.  Roederer,  I.G.  ( 1 ' 1 7 0 )  Dynamics  of  (ieoniagneticallv  I  rapped  Radiati  >n. 

Springer  Verlag,  New  York, 
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l-iaure  5.2.  A  Charged  {’article 
Tra iect.»rv  in  a  Magnetic  Bottle’  . 
Conservation  >1'  the  lirst  a  liahath 
invariant  can  cause  the  spiraling 
particle  to  he  reflected  where  ihe 
magnetic  field  is  stronger.  This 
.  ;iuses  the  particle  to  he  trapped  hv 
i  he  magnetic  I'iel.i 


To  a  lair  approximation,  the  earth's  magnetic  i'iel.i  in  the  radiation  belt  region 


can  be  described  in  terms  of  a  magnetic  dipole  located  near  the  center  of  die  earth. 


.2 


The  dipole  moment  is  M  0.2  12  G  HJ.,  and  the  dipole  is  directed  so  that  the  man¬ 


tlet  it:  south  pole  on  the  earth's  surface  is  located  in  northern  Greenland  (geographic 
coordinates:  78.  5*  X,  20  1’  Gl;  on  the  earth,  the  northern  end  of  the  compass  nee  Hr 
ooints  to  this  location.  The  spatial  distribution  of  the  dipolar  magnetic  field  strengt 
beyond  the  surface  of  the  earth  is 
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cos  A 


where  H  is  the  radial  distance  measured  from  the  center  of  ihe  earth,  H  ,  is  the 
radius  of  the  earth,  H  0.2  12  G  is  the  equatorial  field  at  II  II  ,  and  A  is  Ihe 
magnetic  latitude.  A  detailed  account  of  the  earth's  magnetic  field  and  its  variabilit 
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is  found  in  I'hapman  and  Martels.”  Sin  -c  the  ;  :■  i  i  '  -  ml.  .<•  •  a  , 

a  radiation  belt  particle  c\i>erieni  cs  v  irvin..  n  nijn."'.  :.■!  j  st  ret.'*  b-  \ :  ■- 

trajectorv.  Kiel. I  variations  m  i  homth  scale  r  t '••.<•  •  t.-r  >:  "i.-  !•-.:•• .  n  1  ■ 

radius  cause  a  net  drift  ross  the  magnet  i  - 1 •  •  l  i  ••>.  •tie  17 1 ’.r.  :  , i  i  i :  :  •  -  -1 

around  the  earth  as  illustrated  in  I  inure  a.  :t.  litis  is  a  hr.  '  result  >!  the  n  l  i 
strength  being  greater  .  loser  to  the  earth.  causing  tite  n:i rt i ■  •  orbital  r  t  Ins  c 
curvature  to  be  less  there.  The  direction  »f  tlie  magnetic  force  depen  is  >■>  Tie 
sign  of  the  particle  charge;  electrons  drift  eastward  and  (positive)  ions  iritt  .test- 
ward.  1  hus  the  energetic  trapped  parti,  les  are  spread  out  in  a  belt  like  configura 
tion  around  the  earth,  forming  the  radiation  belts. 
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Figure  5.3.  Schematic  Representation  >f 
the  Gvration  and  Azimuthal  Drift  (Solid 
Curve)  of  an  Equatoriallv  Mirroring  Proton, 
With  Associated  Current  Patterns  (Dashed 
Curves  r 
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Belts,  Springer  \  erlag,  New  York.  _  ” 


Radiation  belt  oli'i  lruns  la  ni'  at  vim'v  high  speeds.  I-  ir  example.  eloi  -irons 
with  a  kinetic  energy  of  500  kc\  move  at  85  percent  of  the  speed  of  light.  Radiation 
licit  electrons  must  thus  be  studied  vising  relativistic  thoorv.  Ions,  because  thev 
are  substant  ia  1 1  v  heavier,  genera  lie  move  at  subrelativistic  velocities:  at  500  kc\ 
the  proton  speed  is  8  percent  of  the  speed  of  light,  while  the  heavier  ions  are  even 
slower  at  the  same  energy. 

The  composition  and  flux  intensities  of  the  earth's  radiation  belts  are  leter- 
niined  bv  the  strength  of  the  sources,  internal  transport  processes,  and  loss 
mechanisms.  The  outer  boundary  of  the  radiation  belt  trapping  region  occurs  at 
'he  point  where  the  magnetic  field  is  no  longer  able  to  maintain  stable  trapping, 
and  at  low  altitudes  the  earth's  atmosphere  forms  an  effective  boundary  for  radia- 
'  ion  bell  particles.  I  luring  geomagneticallv  quiet  conditions,  the  radiation  belt 
region  extends  from  the  top  of  the  atmosphere  along  magnetic  field  lines  to  an  equa- 
' orial  radial  distance  of  at  least  7  earth  radii.  Knergetie  ions  and  electrons  that 
encounter  the  dense  atmosphere  collide  with  the  atmospheric  constituents  and  are 
readily  lost  from  the  radiation  belts.  Particles  with  mirror  points  well  inside  the 
atmosphere  (nominally  below  ~  100  km  altitude)  are  said  to  be  within  the  atm  ispheric 
fb ounce)  loss  cone.  Kinematically,  the  angle  between  the  velocity  vector  of  such 
narticlcs  and  die  magnetic  field  direction  (the  pitch  angle)  at  the  equator  is  below  a 
ertain  value,  known  as  the  loss  nine  angle.  Figure  5.4  exemplifies  this. 


f  igure  5.4.  Illustration  of  Mag¬ 
netic  Mirroring  in  a  Dipolar 
Magnetic  Field.  The  single  par¬ 
ticle  trajectory  sh  nvn  m  solid  line 
is  for  a  particle  outside  the 
atmospheric  bounce  loss  .  one  and 
the  dashed  line  represents  the 
trajectory  of  a  particle  inside  the 
loss  cone.  The  latter  particle  will 
encounter  the  denser  parts  of  the 
earth's  atmosphere  (mirror  point 
height  nominally  below  100  km), 
and  will  thus  precipitate  from  the 
radiation  belts. 


The  earth's  radiation  environment  is  best  studied  bv  combined  experimental  anti 
theoretical  means.  On  one  hand,  it  is  impossible  to  encompass  the  entire  mag¬ 
netosphere  bv  experimental  techniques  or  even  to  measure  all  the  physical  param¬ 
eters  that  may  have  bearing  on  the  dynamical  phenomena;  on  the  other  hand,  this 
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environment  is  so  complex  that  there  can  be  little  hope  of  theoretically  predicting 
the  total  radiation  belt  behavior  solely  from  a  set  of  mathematical  postulates.  1  or 
example,  the  governing  diffusion  equations  describing  the  trapped  radiation  phenomen; 
may  be  known  in  analytic  form,  but  the  transport  coefficients  that  enter  into  them 
must  be  empirically  determined. 

This  report  discusses  these  and  other  applied  concepts.  The  current  state  of 
knowledge  of  the  geomagnetically  trapped  radiation  is  described  both  from  the 
theoretical  perspective  and  from  direct  observations.  We  demonstrate  how  this 
knowledge  is  used  to  construct  physical  models  of  the  radiation  belts.  Derivations 
of  all  governing  physical  relations  from  first  principles  is  beyond  the  scope  of  this 
report.  However,  we  provide  references  to  a  number  of  books  and  research  articles 
where  such  information  can  be  found.  Empirical  radiation  belt  models,  based  on 
data  compiled  from  many  spacecraft  and  a  brief  survey  of  man's  interaction  with 
geospace  are  also  presented. 

5.1  THEORETICAL  PRELIMINARIES 

The  difficult  mathematical  problem  of  the  motion  of  energetic  charged  particles 

in  a  dipolar  magnetic  field  was  extensively  studied  during  the  first  half  of  the 

twentieth  century.  A  general  analytic  solution  to  the  equation  of  motion  was  never 

found,  and  in  most  cases  particle  orbit  tracing  had  to  be  done  numericallv.  The 

4 

interested  reader  is  referred  to  Stormer.  Physical  approximations  that  lead  to 
great  simplification  have,  however,  been  found.  This  is  known  as  the  adiabatic 
theory  for  trapped  particles,0  and  the  earth's  radiation  belts  have  now  been  success¬ 
fully  described  in  terms  of  adiabatic  invariants  and  their  perturbations. 

5.1.1  Single  Particle  Motion 

An  ensemble  of  ions  and  electrons  moving  in  space  constitutes  a  plasma  that 
can  exhibit  many  modes  of  collective  as  well  as  single  particle  behavior.  In  the 
presence  of  electric  (E)  and  magnetic  (§)  fields  these  particles  are  subject  to  the 
electromagnetic  Lorentz  force,  E  =  q  (E  +  v  X  B),  where  q  and  v  are  the  particle 
charge  and  velocity  vector  respectively.  For  ions  q  =  Ze,  where  Z  is  the  ionic 
charge  state  and  e  is  the  unit  charge;  for  electrons  q  =  -e.  This  force  controls  the 
particle  motion,  and  collectively  the  ensemble  of  charged  particles  can  modify  the 
fields  through  induction,  charge  separation  and  electrical  currents  formed  by 


Sformer,  L.  (1955)  Polar  Aurora,  Clarendon,  Press,  Oxford. 

Alfven,  H. ,  and  Falthammar,  C.G.  (1963)  Cosmical  Electrodynamics, 
Clarendon  Press,  Oxford,  United  Kingdom. 


5.1.2  Adiabati.  Invariants 


In  general,  the  motion  of  charged  particles  is  such  that  momentum  and  energy 
can  be  transferred  between  the  different  particles,  and  between  the  particles  and 
the  fields  that  influence  their  motion.  Therefore,  it  is  not  always  possible  to 
identify  constants  of  motion.  However,  under  certain  conditions  these  energy  and 
momentum  exchanges  are  very  small,  and  it  is  possible  to  identify  specific  quan¬ 
tities  that  remain  virtually  unchanged  with  the  particle  motion.  These  are  called 
adiabatic  invariants. 

Associated  with  each  of  tie  three  quasi-periodic  modes  of  motion  is  n  adiab.itr 
invariant  related  to  the  Hamilton- Jacobi  action  variable: 


J.  =  j>  |P  +  9.  A]  •  (U  (5.  2) 

i  c 

where  ti  i  is  a  vector  line  element  along  the  path  >f  integration.  Here  H  is  particle 

-J 

momentum,  A  is  the  magnetic  vector  potential  (that  is.  B  -  v  «  A)  and  e  is  the 
speed  of  light.  The  integration  is  extended  over  the  particle  orbit  for  gvro  motion, 
bounce  motion,  and  azimuthal  drift  motion  (for  i  •  1,  2,  J  respectively).  If  the 
particle's  trajectory  closed  exactly  or  itself,  then  the  action  variables  .1.  would  be 
absolute  constants  of  motion.  Finite  spatial  and  temporal  v  ariations  in  15  prevent 
perfect  closure,  and  thus  the  J.s  are  at  best  approximate  constants. 

5.  1.  2.  1  FIRST  ADIABATIC’  IX\  A  lit  A  XT 

is  obtained  by  evaluating  the  integral  in  Fq.  (5.  2)  over  the  particle  g\  ro 
motion  only,  that  is,  over  the  particle  orbit  projection  in  a  plane  perpendicular 
to  B.  Using  subscripts  n  and  ±  to  denote  directions  parallel  and  perpendicular 

-4  3 

to  B,  and  by  virtue  of  Stokes'  theorem  (Schulz  and  l.anzerotti  ]  one  derives 


J!  =  2^g  Pi  +  Z  V  Pg  B 


(5.  :i) 


where  15=  |b|,  and  q !  13  is  the  particle  gyro  (or  cyclotron)  radius.  From 

tliis,  one  defines  the  first  adiabatic  invariant 
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*  1  _  p  sin  a 

2  m  B  -  2  m  fi  ’ 

o  o 


(a.  4) 


also  known  as  the  relativistic  magnetic  moment.  Here  ni(|  is  the  particle  rest  mass 
and  a  -  arc  sin  (it  *  p/pB)  is  the  pitch  angle  illustrated  in  Figure  5.6. 
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Figure  5.6.  A  Particle  With  Its 
Velocity  Vector  Inclined  to  the 
Magnetic  Field  at  an  Angle  a. 
This  angle  is  called  the  particle 
pitch  angle 


For  non- relativistic  particles, 

_  1  2 

£.1  7  mo  v_l 

M  =  =  - B - 


(5.  5a) 


where  £j_  is  the  particle  kinetic  energy  associated  with  the  directions  perpendicular 
to  the  local  magnetic  field  direction,  while  for  relativistic  particles, 
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(5.  5b) 


where  the  relativistic  factory  =  1  /\\  -  vie.  p  is  an  approximate  constant  of 
motion  when  both  of  the  following  conditions  are  fulfilled: 


—  The  spatial  scale  of  B-field  variation  is  much  larger  than  the  particle  gyro- 
radius 


JL  >  > 
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p  =  p  /B  |q|  =  7  m  v  sina/B  Iql  . 


—  The  time  scale  of  change  of  the  B-field  is  much  larger  than  the  particle  gyro- 
period 


T  >  >  r  g  =  2n  PRI  v±  =  2v  y  mQ/  B  lq| 


5.  1.  2.  2  SECOND  ADIABATIC  I XV  A  HI  A  XT 


J 2  is  obtained  by  evaluating  the  integral  in  Eq.  (5.  2'  over  the  bounce  trajectory 
and  averaged  over  the  gyro  motion,  or  equivalently  along  the  magnetic  field  line 
(guiding  center  field  line)  around  which  the  particle  gyrates,  thereby  defining  the 
second  adiabatic  invariant: 


AAA 


■aaafidBfius 


if).  8) 


J  =  =  i  4  ‘  lli?  =  /  pn d; 


where  d is  an  element  of  length  along  that  field  line  segment  and  is  the  curvi¬ 
linear  distance  of  the  mirror-points  from  the  equator  measured  along  the  guiding 
center  magnetic  field  line.  Since  equatorially  mirroring  particles  (<»  -  rr/2'  do  not 

have  any  bounce  motion,  it  follows  that  J  =  0  for-  such  particles. 

Provided  the  particle  mirror-points  are  above  the  dense  atmosphere,  I  will 
remain  an  approximate  constant  when  the  time  scale  of  11-field  variations  is  much 
larger  than  the  particle  bounce  time  between  the  conjugate  mirror-points 
■  £ 
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Constancy  of  the  first  adiabatic  invariant  p.  implies  that 
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sin  a 
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sin  a 


(5.  10) 


over  the  bounce  motion  between  the  mirror  points.  Here,  the  subscript  zero  denotes 
equatorial  quantities  and  B  is  the  magnetic  field  induction  at  one  of  the  mirror 
points  (where  a  =  u/2).  Equation  (5.  10)  is  known  as  the  mirror  equation.  Using 
Eq,  (5.  10)  one  finds 


m  > 


dMl  '  TT"  sin  ao> - jr  f  (1  ■  TT~) 
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(5.  11) 


and  in  dipolar  coordinates  [Eq.  (5.  1 1 ) |  becomes 
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where  T(ao)  is  the  bounce  time  integral  given  by 
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(5.  13) 


and  A  («o>  is  the  magnetic  latitude  of  the  mirrorpoint,  which  depends  on  the  equa¬ 


torial  pitch  angle  a^. 


To  a  fair  approximation 


Tuy  >  '  1.  30  -  0.  56  sin  iy  1  Hamlin  et  al^  I  i5.H) 

o  o  ' 

or  alternatively 

T(cy  )  *  1.3802  -  0.3198  (sin  a  +  |stn  a  ]'  “>  [Schulz  and  1  .anze rotti  *1  . 

( a .  1  a  i 

Other  approximations  are  given  by  Davidson.  ‘  Gradients  of  Tuy  '  should  not, 
however,  be  derived  from  such  approximations.  I  rom  Kqs.  (’>.  1)  and  (a.  101  it 
follows  that 

B  B 

Bm  =  — r =  ,3.2 -  •  (ri-  lb) 

sin  «  L.  sin  a 
o  o 


where  in  dipolar  coordinates,  L  =  R/R^  (measured  at  the  equator),  and 


sin  «o  = 


cos65  A  (a  ) 
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(5.  17) 


Although  Eq.  (5.  17)  cannot  be  solved  explicitly  for  the  mirrorpoint  latitude  -xni,0'o'> 
a  numerical  solution  is  easily  obtained,  or  one  may  approximate,  as  in  Hamlin  et  al: 


cos  A  uy  '  »  | sin  a 
m  o  1  o 
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(5.  18) 


It  should  be  emphasized  that  using  a  dipolar  magnetic  field  representation 
explicitly  disregards  any  azimuthal  asymmetries  of  the  geomagnetic  field.  Such 
asymmetries  do  exist  and  become  significant  beyond  l.=  5.  Under  such  conditions 
a  different  magnetic  field  representation  should  be  used,  and  this  is  outlined  in 
Section  5.6.  1.4. 

5.  1.  2.  3  THIRD  ADIABATIC  INVARIANT 

is  obtained  by  evaluating  the  integral  in  Eq.  (5,  2)  over  the  particle  drift 
motion  around  the  earth,  and  averaged  over  gyro  and  bounce  motion 
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where  4>  is  the  magnetic  flux  enclosed  by  the  azimuthal  drift  orbit  and  1  is  linear 
azimuthal  distance.  Using  Stokes'  theorem  yields 

4>  -  A  ■  1  T  =  /  (i  •  d  S  (5.  20) 


where  S  is  a  surface  bounded  by  the  azimuthal  drift  path.  In  a  dipolar  magnetic 
field  one  calculates  |  Koedere r]  * , 


where  1.  is  the  lUcllwain 
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1.  parameter. 


J  will  remain  approximately  constant  when  the  time  scale  of  15-field  change  is 

much  longer  than  the  azimuthal  drift  time  7  ,  =  <f>  ds/V  ,  around  the  earth.  Ueter- 

0  drift  a 

mination  of  the  azimuthal  drift  velocity  is  discussed  in  the  next  section. 


5.1.3  Particle  Drift  Motion 

In  a  uniform  magnetic  field,  charged  particles  execute  a  spiral  motion  such  that 
the  angle  between  the  particle  velocity  vector  and  the  magnetic  field  direction  (the 
pitch  angle)  remains  constant.  When  the  magnetic  field  lines  converge,  the  particle 
will  respond  to  an  effective  net  magnetic  force  from  higher  to  lower  magnetic  field 
strengths.  This  is  illustrated  in  Figure  5.7.  The  physical  reason  for  this  force  is 
that  the  particle  gyro  motion  produces  an  elementary  current  (which  may  be  inter¬ 
preted  as  a  magnetic  dipole  current  loop).  For  each  such  loop  the  effective  current 
is 


i  -  dq  =  qp-i 
P  dt  2-pgmoy  ’ 


(5.  22) 


where  p  ^  is  the  mean  gvroradius  over  the  loop, 
loop  enclosing  an  area  A  is 
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The  magnetic  moment  of  a  current 


(5. 23  I 
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trapped  particles,  J.  Cleophys.  Res.  66:3681. 
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Figure  5.7.  (a)  Particle  Motion  in  a  Uniform  Magnetic 

Field  (Uniform  Spiral  Motion);  (b)  The  Tightening  of 
Spiral  Motion  in  a  Converging  Magnetic  Field;  and  (c) 
Illustration  of  Magnetic  Force  With  Gyroaveraged  Net 
Component  in  the  -  VB  Direction  in  a  Converging 
Magnetic  Field 


where  A  =  tt  p  ,  which  is  the  particle  magnetic  moment  itself.  The  particle  will 

§  -*  -*  |  -*  | 

therefore,  averaged  over  its  gyro  motion,  be  subject  to  a  net  force  F  =  -  M  "n  |  H 

in  the  direction  along  the  field  lines  away  from  the  higher  field  region. 

In  general,  the  magnetic  field  may  also  have  an  intensity  gradient  across  the 
field  lines.  This  is  illustrated  in  Figure  5.8.  Charged  particles  moving  in  such 
a  magnetic  field  will  have  a  smaller  gyroradius  in  the  higher  field  region  and  a  larger 
gyro  radius  in  the  lower  field  region.  As  a  consequence,  there  will  be  a  net  drift 
velocity  perpendicular  to  the  magnetic  field  direction.  In  this  figure,  a  positive 
charge  would  drift  downwards  (into  the  paper),  and  a  negative  charge  upwards  (out 
of  the  paper).  Defining  the  angular  gyrofrequeney 
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one  can  express  the  instantaneous  vector  gyroradius  as 


-♦  p  X  B  p  y  B 

pg~  moyftB  =  2 


and  the  drift  velocity  is  then  the  time  rate  of  change  of  p^ 
d  c  dp  - 

vd -df  '  At  hf '■  u’ 


(5.  25) 
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Figure  5.  8.  (a)  A  Uniform  B-l'ield  Represented  bv  Fvenly 

Spaced  Field  Fines  and  (b>  A  Magnetic  f  ield  With  Increasing 
Strength  Perpendicular  to  B,  Represented  as  Denser  Field 
Fines  With  Increasing  B 


where  it  is  assumed  that  the  magnetic  field  is  constant  in  time.  With  I-  =  being 
the  net  force  due  to  the  cross- B  gradient,  one  obtains  the  gradient  drift 
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which  non- relativistieallv  is  just 
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(5.  27  bl 


The  earth's  magnetic  field  is  also  curved  (that  is.  the  dipolar-like  field  line: 
form  loops  from  pole  to  pole),  and  the  field  line  radius  of  curvature  is  given  by 
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where  the  individual  field  lines  are  described  by  the  dipole  relation 


R  =  R  cos  A 
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with  R  =  R..I-,  Thus  a  charged  particle  moving  in  that  field  will  experience  a 
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centrifugal  force 
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which  non- relativistic. >  II .  becomes 
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Although  the  effect  of  the  earth'd  gravitational  field  id  rather  small  compared 
to  other  forces  on  radiation  belt  particles,  it  can  easily  be  included: 
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where  g  is  the  vector  gravitational  acceleration. 

The  effect  of  a  weak,  externally  imposed  electric  field  is  also  easily  taken  into 
account: 
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The  t  leetric  field  drift  is  independent  of  particle  charge  and  mass  as  long  as  either 
is  non-zero.  Thus,  under  the  influence  of  an  electrostatic  field,  the  electric  field 
.  or. portent  O'-rnendii  ular  to  the  n  agnetic  field  direction  causes  ions  and  electrons 
to  drift  together  (plasma  flow),  while  under  the  influence  of  an  inhomogeneous  mag¬ 
netic  field,  oppositely  charged  particles  drift  in  opposite  directions  (causing  current 
flow). 

The  total  particle  drift  velocitv  is  then  the  superposition  of  tiie  contributing 
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In  the  dipolar  magnetic  field  representation,  an  approximate  formula  for  the  drift 
period  is  given  by  Davidson 


1. 43  K, 
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10  sec  for  electrons,  K  =  5.055  sec  for  protons, 
for  ions  of  mass  nr  and  charge  state  Z  . 


and 


The  cartoon  in  Figure  5.  9  illustrates  the  principal  drift  effects  associated  with 
tile  different  drift  mechanisms.  In  Kqs.  <5.  36)  and  (5.  37)  the  tern  s  are  listed  in 
order  of  their  importance  in  the  radiation  belts.  Above  '10  k.eV  the  magnetic 
gradient-curvature  drift  is  generally  strongest,  apd  static  electric  field  and  gravity 
effects  are  usually  neglected  in  radiation  belt  studies.  The  gradient-curvature  drift 
carries  energetic  electrons  towards  the  east  and  ions  to  the  west.  Thus  there  will 
be  a  net  westward  electrical  current  encircling  the  earth.  This  is  the  extraterres¬ 
trial  ring  current.  These  findings  are  summarized  in  Figure  5.  10. 


5.2  TRAPPKO  RADIATION  SOl’RCKS 

Precisely  where  the  radiation  belt  particles  come  from  and  how  they  are 
accelerated  to  energies  in  the  ke\  and  Mel  range  are  still  areas  of  research  for 
which  a  comprehensive  answer  is  not  yet  available. 
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Figure  5.9.  Summary  Illustration 
of  Drift  of  Doth  Positive  and 
Negative  Charged  Particles 
in  (a)  a  Uniform  15-1  ield.  No 
External  Force;  (b)  a  Uniform 
B-Field  With  a  Perpendicular 
Electric  Field;  (c)  a  Uniform 
B-Field  With  an  External  Force 
That  is  Independent  of  Electric 
Charge,  Such  as  Gravity;  and 
(d)  a  B-Field  U  ith  a  Gradient 


Figure  5.  10.  A  Descriptive  Drawing  of  the  Three  Types  of  Motion  of  Particles 
Trapped  in  the  Earth's  Magnetic  f  ield 


5.2.1  Qualitative  Description 
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A  number  of  sources  are  considered  responsible,  and  the  effectiveness  of  each 
probably  also  varies  with  time: 

til  Particles  from  the  sun,  including  solar  wind  particles  ..nil  energetic  solar 
particle  emissions,  possibly  via  magnetotail  storage. 

(2)  Particles  from  the  earth's  ionosphere,  including  the  polar  wind  flow  into 
the  magnetotail,  particles  flowing  up  the  magnetic  field  lines  to  form  the  plasma- 
sphere  and  particles  accelerated  out  of  the  auroral  ionosphere. 

(3)  Cosmic  ray  albedo  neutron  decay  (GI1AND)  within  the  trapping  region. 

(4)  Particles  arriving  at  the  earth  having  been  accelerated  in  interplanetary 
shock  waves  or  in  the  magnetospheres  of  other  planets. 

(5)  Low  energy  components  of  galactic  cosmic  rays. 

(6)  In  situ  acceleration  of  pre-existing  lower  energy  trapped  particles  within 
the  radiation  belts. 

The  solar  wind  flows  past  the  earth's  magnetosphere  virtually  at  all  times. 

Some  of  these  particles  may  find  their  way  through  the  outer-  regions  of  the  mag- 

9 

netosphere  to  the  stable  trapping  region;  (Hovestadt  et  aT)  this  process  may  be  par¬ 
ticularly  effective  during  periods  of  southward  heliospheric  magnetic  field.  Direct 
transient  injections  of  solar  energetic  particles  probably  also  occur,  particularly  in 
conjunction  with  magnetic  storms. 

Ionospheric  particles  diffusing  out  of  the  polar  ionosphere  (polar  wind)  escape 
into  the  magnetotail  region  from  which  some  mav  become  energized  and  injected 
into  the  trapping  region.  Auroral  electric  fields  are  intermittent  and  can  have  a 
significant  component  parallel  to  the  magnetic  field,  and  ions  and  electrons  from  the 
topside  auroral  ionosphere  can  be  accelerated  to  multi -ke\  energies.  The  wave 
fields  associated  with  plasma  waves  may  also  cause  particle  acceleration.  This 
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could  be  a  source  of  H  ,  He  ,  ()  ,  and  electrons  provided  other  processes  act  to 
trap  the  particles. 

Cosmic  rays  impacting  the  earth's  atmosphere  undergo  nuclear  reactions,  and 
a  flux  of  neutrons  escapes  from  the  top  of  the  atmosphere.  Tree  neutrons  are 
unstable  and  decay  into  a  proton,  electron  and  neutrino  triplet  on  a  time  scale  of 
~  1000  seconds.  If  the  electrically  charged  decay  products  find  themselves  within 
the  radiation  belts,  they  will  immediately  be  subject  to  the  magnetic  force  and  may 
become  trapped.  Empirically,  this  is  an  important  source  of  multi-IUeV  protons 
in  the  innermost  part  of  the  inner  radiation  zone  (I.  *  1.5).  There  is  little  direct 
information  about  the  efficiency  of  direct  extraterrestrial  energetic  particle  trapping 

Hovestadt,  D.  ,  Gloeckler,  (4. ,  Tan,  C.  Y.  ,  Disk,  I..  A.,  Ipavich,  1  .  A I .  , 

klecker,  li.  ,  <  ('Gallagher,  .1..!.,  and  Scholer,  M.  tlf'TH)  Evidence  for  solar 
wind  origin  of  energetic  ions  in  earth's  radiation  belts,  Geophvs.  lies.  I.etts. 
5:1055.  ‘  — 
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in  the  radiation  belts.  One  may  surmise,  however,  that  time  variability  of  the 
geomagnetic  field  is  needed  for  trapping  to  occur  or  that  incident  extraterrestrial 
energetic  ions  in  low  charge  states  (for  example  (.  '  and  i>  i  may  chare*  exchange 
to  higher  charge  states  (for  example,  C  +  and  0^*1  within  the  magnetosphere  so 
that  their  gyroradii  become  small  enough  (/>  -  p^c  (HI  for  trapping. 

The  plasma  sheet  in  the  earth's  magnetotail  is  considered  a  possible  reservoir 
for  radiation  belt  particles.  However,  the  plasma  sheet  particles  themselves  are 
likely  to  be  a  mixture  of  particles  from  several  of  the  prime  sources  mentioned 
above.  During  magnetospheric  substorms  plasma  sheet  particles  convect  inward 
toward  the  earth,  and  in  the  process  can  be  accelerated  and  may  become  t  -apped 
in  the  radiation  belts.  Unfortunately,  little  is  known  about  the  specific  details  of 
the  time-dependent  trapping  process  and  it  is  not  yet  possible  to  make  a  quantitative 
evaluation  of  the  strength  and  characteristics  of  this  source. 

Current  research  also  points  to  the  earth's  ionosphere  as  an  important  contribu¬ 
tor  to  the  lower  energy  particle  population  i below  a  few  tens  of  keV)  perhaps  with  a 
roughly  equal  contribution  from  solar  wind  particles.  In  contrast,  the  high  energy 
particles  (above  '  500  ke\  '  appear  to  have  an  extraterrestrial  source.  The  former 
conclusion  is  derived  from  observations  of  dominant  oxygen  fluxes  at  times,  while 
the  latter  stems  from  observed  ca rbon-to-oxygen  ratios  of  order  unity,  l  or  the 
very  important  intermediate  energy  range  where  most  of  the  radiation  bell  energy 
density  is  found,  there  is  no  experimental  result  indicating  the  source. 

A  simple  theory  that  seeks  to  explain  observed  storm  time  enhancements  of 
radiation  belt  particle  fluxes  at  tens  and  hundreds  of  kcV  energies  has  been  advocate 
by  Lyons  and  Williams.  10  Particles  existing  in  the  outer  radiation  .’.one  may 
suddenly,  during  the  storm  main  phase,  be  subiect  to  an  electric  field  that  trans¬ 
ports  them  towards  lower  L- shells  on  a  time  scale  that  preserves  /.  and  .1  but 
violates  the  constancy  of  4*.  For  equatoriall.v  mirroring  particles,  a  radial  displace 
ment  from  I..  =  5  to  L  =  3  increases  the  particle  energy  by  a  factor  of  '  5.  further¬ 
more,  if  the  particle  spectrum  follows  a  power  law  distribution  j'P'  =  <L  I.  }  U  with 
u  =  3,  for  example,  then  the  apparent  particle  flux  enhancement  seen  at  fixed 
energies  will  be  a  factor  of  125.  Seen  at  a  fixed  1  -shell,  the  observable  tiux 
increase  will  be  even  greater  if  the  pre-storm  radial  distribution  '.alls  off  tovard.- 
lower  I, -shells  (as  may  be  the  case  at  lower  radiation  belt  energies),  and  smalle- 
if  this  flux  gradient  is  negative.  Given  the  knowledge  of  the  stern,  perturb -ti.  n 
electric  field  (magnitude,  direction,  azimuthal  extent,  and  duration'  together  wit:-, 
observations  of  the  pre-storm  radiation  belt  structure,  this  in  situ  Sourc*  is  in 
principle  assessable. 


10.  Lyons,  L.  I{.  ,  and  W  illiams,  D.  J.  ll!180)  \  source  for  the  geomagnetic  store 
main  phase  ring  current,  Geonhvs.  Pcs.  0  5:5  2.). 
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.">.2. 2  Simplifying  Assumptions 

For  most  of  the  radiation  belt  source  mechanisms  accurate  quantitative  informa¬ 
tion  is  still  lacking. 

ldeallv,  one  would  like  to  know; 

(1)  The  source  strength  for  different  particle  species  as  function  of  energy 
and  pitch  angle, 

(2)  The  effective  source  locations  within  and  on  the  boundaries  of  the  trapping 
region,  and 

(3)  The  source  strength  as  function  of  the  different  geophysical  conditions 
during  quiet  and  disturbed  times. 

Unfortunately,  contemporary  research  has  not  yet  yielded  quantitative  answers 
to  these  requirements.  Without  this  information,  how  can  we  understand  and  model 
the  earth's  radiation  belts’.’ 

For  quiet  time  conditions,  one  can  solve  the  steady  state  radiation  belt  trans¬ 
port  equations  for  the  interior  of  the  radiation  belts  subject  to  suitable  outer  zone 
boundary  conditions  on  the  trapped  fluxes.  This  amounts  to  the  assumption  that  the 
radiation  belt  source  is  capable  of  supplying  particles  to  the  outer  radiation  zone 
boundary  at  a  rate  sufficient  to  offset  losses  within  the  trapping  region.  The 
existence  of  long-term  approximate  stability  of  the  radiation  belts  as  a  whole  during 
extended  quiet  periods  supports  this  contention,  and  fortunately,  trapped  flux  ob¬ 
servations  from  geostationary  spacecraft,  such  as  ATS-6,  of  the  outer  zone  flux 
levels  at  L  =  6.6  makes  this  a  feasible  solution. 

This  is  not  a  satisfactory  situation  for  magnetic  storms  and  other  disturbances. 
As  a  consequence  most  radiation  belt  modeling  has  been  done  for  steady  state,  quiet 
time  conditions.  Time  dependent  radiation  belt  modeling  would  require  time  de¬ 
pendent  boundary  conditions  averaged  over  local  time.  It  is  conceivable  that  data 
from  several  geostationary  satellites  may  be  used  for  this  purpose,  but  this  has 
not  yet  been  done. 


5.3  TRANSPORT  PROCESSKS  IN  THE  RADIATION  BELTS 

As  we  have  seen,  in  the  static  geomagnetic  field,  radiation  belt  particles 
execute  the  three  periodic  motions:  gyration  around  the  magnetic  field  lines,  bounce 
motion  between  mirror  points,  and  azimuthal  drift  around  the  earth.  The  latter 
type  of  periodicity  is  caused  by  the  gradient-curvature  drift  due  to  the  inhomogeneitv 
of  the  magnetic  field.  Effects  of  electric  fields  and  gravity  cause  departures  from 
this  simple  picture,  but  both  forces  are  primarily  important  at  low  energies,  typic¬ 
ally  below  a  few'  tens  of  keV. 
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5.3.1  Convection 

Alagnetuspheric  convection  results  from  externally  imposed  electric  fields. 

The  solar  wind  flows  past  the  earth  at  a  velocity  in  the  range  200  -  600  km  see. 

-* 

This  implies  a  convection  electric  field  E  -  —  ■  H  across  the  earth's  mag¬ 
netosphere  directed  from  dawn  to  dusk.  Combined  with  the  electric  field  induced 
by  the  rotation  of  the  earth,  a  characteristic  magnetosnherie  convection  pattern  is 
set  up.  The  radial  co rotational  electric  field  is  induced  by  the  earth's  magnetic 
dipole  field  corotating  with  the  earth.  Figure  5.  11  shows  model  electric  equi- 
potential  lines  around  the  earth  due  to  these  electric  fields  in  a  time-independent 
situation.  1'articles  of  quite  low  energies  (  E  1  ke\  l  are  primarilv  controlled  bv  the 
electric  field  drift,  and  their  motion  approximately  follows  the  equipotential  lines 
as  indicated  by  the  arrows  in  f  igure  5.  11.  Notice  the  topologically  distinct  regions: 
near  the  earth  where  the  convective  motion  follows  oval  paths  around  the  earth,  and 
at  greater  distances  where  the  drift  paths  are  open  to  the  magnetopause.  Departures 
from  this  overall  configuration  of  course  take  place  during  disturbed  conditions  when 

the  imposed  "convection"  electric  field  is  time  variable  and  the  low  energe  particle 
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distributions  are  not  in  equilibrium. 


Figure  5.  11.  Kquipotential 
Contours  for  an  Electric 
Field  in  the  Earth's  Equa¬ 
torial  Plane  (Dashed  Lines). 
These  are  also  drift  paths 
for  very  low  energy  particles. 
The  electric  field  is  a  super¬ 
position  of  a  corotational 
E-field  due  to  the  rotation 
of  the  earth  and  its  imbedded 
magnetic  field,  and  a  uni¬ 
form  dawn-dusk  electric 
field.  The  separatrix 
(solid  equipotential  curve) 
is  the  low  energy  particle 
laver  that  separates  the 
'pen  and  closed  drift  paths. 
The  outer  oval  represents  the 
magnetopause  in  the 
equatorial  plane 
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(Due  to  the  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  127.  ) 


5.3.2  Stochastic  Processes:  Phenomenology 


At  higher  energies,  particularly  above  a  few  tens  of  keV,  the  dominant  drift 
of  trapped  particles  is  due  to  the  gradient  and  curvature  effects  of  the  geomagnetic 
field.  The  prime  drift  motion  is  therefore  circular  (with  gyro  and  bounce  motion 
superimposed)  around  the  earth.  Departures  from  this  pattern  are  due  to  the 
fluctuations  of  the  geoelectric  and  geomagnetic  fields  induced  by  variations  in  the 
solar  wind  flow  and  internal  magnetospheric  processes.  Radiation  belt  particles 
are  also  subject  to  interactions  w'ith  plasma  waves  and  suffer  collisions  with 
exospheric  neutral  atoms  and  low  >  aergy  plasma  particles.  Common  to  these  pro¬ 
cesses  is  their  randomness  in  occurrence,  and  their  effects  are  described  by 
stochastic  analysis  that  can  be  reduced  to  diffusion  theory:  radial  diffusion  and 
pitch  angle  diffusion.  The  physical  ideas  are  illustrated  in  f  igure  5.  12.  Radial 
liffusion  transports  radiation  belt  particles  across  the  di polar -like  magnetic  field 
lines  in  the  radial  direction,  and  pitch  angle  diffusion  alters  the  particle  pitch  angle 
(or  equivalently,  the  mirrorpoint  location).  In  both  cases  the  earth's  atmosphere 
is  a  sink:  for  radial  diffusion  by  transport  to  very  !•  w  I. -shells,  and  for  pitch  tingle 
diffusion  bv  lowering  the  mirrorpoints  into  the  atmosphere.  In  addition  to  diffusive 
processes,  energy  degradation  from  collisions  with  exospheric  particles  also  occurs. 


SOURCE 


PITCH  ANGLE 
DIFFUSION 


f  igure  5.  12.  A  Conceptual 
Representation  of  Pitch-tingle  and 
Radial  Diffusion  in  the  Barth's 
Radiation  Belts.  Diffusion  occurs 
in  either  direction,  but  in  most 
cases  there  is  a  net  diffusion 
flux  in  the  indicated  direction 
towtirds  the  earth  s  atmosphere 


5.3.3  Effects  of  Field  Fluctuations 

Adiabatic  invariants  are  useful  substitutes  for  particle  constants  of  motion.  It 
should  be  kept  in  mind,  however,  that  the  earth  s  magnetic  field  is  never  perfectly 
static,  i  ield  fluctuations  are  associated  with  micropulsations,  magnetospheric 
substorms,  geomagnetic  storms,  and  other  phenomena;  and  field  oscillations  vary 


15.  Jacobs,  J.A.  (  1970)  Geomagnetic  Micropulsations.  Springer  Verlag,  New  York. 

16.  Akasofu,  S. -I.  (1968)  Polar  and  Magnetospheric  Substorms,  D.  Reidel 

Publishers,  Dordrecht,  Holland. 
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considerably  in  magnitude,  frequencies,  and  principal  location.  The  idiahatm  in¬ 
variants  are  said  to  be  violated  when  electric  or  magnetic  v;  iations  take  place 
near  or  above  the  adiabatic  motion  frequency  in  question,  lvven  slow  field  valua¬ 
tions  may  violate  the  third  adiabatic  invariant  4  'T  -  minutes  to  hours)  while 
gl -violation  requires  wave-like  fluctuations  on  a  time  scale  of  milliseconds.  Macro 
S topically,  the  earth's  radiation  belts  are  subject  to  field  fluctuations  til- it  occur 
at  quasi-random  times.  Their  effects  are  hesi  desc  ribed  by  stochastic  methods' 
that  treat  the  mean  deviations  in  the  adiabatic  invariants  and  tile  associated  distribu 
tion  function. 

Wo  deli  u-  a  particle  distribution  function  f  -  f(ii,  .1,  4>;  ll  such  ilia'  in  a  \  diinie 
element  of  parameter  s|>aee  (d,i,  -1.1,  d<l»).  the  number  >f  particles  is  riven  bv 

d.\  -  f{  x .  .1,  <1»:  t)  da  d.l  d4>  :p.) 

at  a  time  t. 

let  _P  ip.  ,  ,1,  <l>,  Ap  ,  A.),  A4>)  be  the  probability  that  a  mean  change  Ap  ,  A. I 
and  A4-  takes  place  in  the  adiabatic  invariants  per  unit  time  (  I  Ap  I  p 
I  A J  !  v  <  .1  and  |a4>!  <<  4>).  The  distribution  function  f(p  ,  J,  4-;  t)  averaged  over 
gyrophase,  bounce  phase  and  azimuthal  drift  phase  is  then  governed  by  the  relation 

f(p  ,  J,  <t>;  t  >  =  ///  d(A p)  d(A  J >  d(A4>>  f( p  -  Ap ),  J-AJ, 

4>-A4>;  t  -  At)  P  (p  -  Ap,  J-AJ,  4>-A4>;  Ap,  AJ,  A4>)  (a.  10) 

and  one  expands  f  anti  1’  in  Tavlor  series  around  the  unperturbed  quantities; 


+  higher  order  terms 

where  the  stochastic  coefficients  are  defined  by 

<Ai  •  =  /  /  /  d(Ap  I  d<  AJ)  d(A4>)  JPp,  J,  4>;  A/i,  AJ,  A4>)  Ai 

<  AiA  j  .  / //d(Ap)d(A.lld(A4*)  b((i,  <t>-  Ap ,  A  J .  A4>1  A  iA  ,j  (5.  42) 

17.  Chandrasekhar,  S.  (19135)  Plasma  Physics  (Compiled  by  S.  K.  Trehan), 

1’hoenix  Science  Series,  University  of  Chicago  Press,  Chicago. 


with  i  and  j  being  /i,  J  or  $  in  all  permutations. 

The  transport  equation  [Eq.  (5.41)]  can  be  greatly  simplified  by  recognizing 
that  violation  of  one  adiabatic  invariant  by  a  process  is  almost  always  uncorrelated 
with  the  process  violating  another.  In  that  case  all  cross  coefficients  must  vanish: 

<A|iAJ>  =  <A/iA<i»  =  <AJA4>>  =  0  .  (5.43) 

Furthermore,  <Ai>  and -j  <(Aii2  (i  =  4,  J,  <t>)  are  related.  In  the  absence  of 
external  sources  and  losses,  diffusion  would  proceed  to  transport  particles  away 
from  overabundant  regions  of  parameter  space  until  all  gradients  in  the  distribution 
function  had  vanished,  and  for  each  diffusion  mode 

2 

< A i  >  -  z.  0  (i  =  4,  J.  <i>)  .  (5.  44) 


Equation  (5.41)  then  simplifies  to 
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where 
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(5.  46) 


Equation  (5.45)  is  the  pure  diffusion  equation  describing  diffusive  transport  in 
the  earth's  radiation  belts.  It  is  valid  whenever  the  perturbations  are  small  (but 
accumulative).  Transformation  to  other  variables,  such  as  Oj,  <5.,,  0,^  is  facilitate 
by  the  Jacobian  G  =  G(4,  J,  4>;  bj,  0^,  0^)  (Reference  3)  such  that 
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and  F 


F(0  j,  <ft2.  03;  t)  . 


5.3.4  Radial  Diffusion 


Field  fluctuations  on  a  time  scale  comparable  to  the  azimuthal  drift  time  around 
the  earth  can  violate  the  third  adiabatic  invariant  4>but  preserve  the  4  and  J 
invariants.  This  is  the  limit  of  pure  radial  diffusion.  To  study  this  process  in 


\  • 


radiation  belt  physics  requires  knowledge  >i  (be  IM'Hismn  rui'Hirii'iit,  as  well  as  ihe 
particle  sources  and  losses.  It  is  convenient  to  use  the  dipole  I. -shell  coordinates 
defined  through 


2  ‘d; 

4>  =  2  -  K“  ~ 


<  a.  4 


where  B^.  -  0.312  G  is  the  equatoidal  B-vahie  at  the  surface  of  the  earth  where 
I.  =  1.  In  this  case  the  Jacobian  coordinate  transformation  is  just 

G(I.;<t>)  =  =  -  277  RJ\  Bj..  — Jy-  .  (a.  501 

and  the  radial  diffusion  equation  becomes 

y f  2  9  r  ,  -2  a f  1  .  „ 

Ft  =  L  FT  [Ul.l.  L  FTj+-'-L* 

where  and  _L  represent  particle  source  and  loss  functions.  IJ.  j  is  the  pure  radial 
diffusion  coefficient  at  constant  p  and  J  values.  It  has  been  estimated  that 


'Vl.  -  UI.I.lao)  "  ULL.(2)d(0!o) 


(5.  5  2) 


where  d(a())  is  a  function  only  of  the  particle  equatorial  pitch  angle  and  / ^  (^  ) 
is  the  radial  diffusion  coefficient  for  equatorially  mirroring  (a^  =  -j)  particles 
where  d(«o  =  =  1.  [tor  a  description  of  d(aQ)  see  Reference  18.) 

Both  geomagnetic  and  geoelectric  field  fluctuations  contribute  to  .  l  or 
geomagnetic  fluctuations  it  can  be  shown  that 


..(Ml,  .  .  2  ttUVI)  ,  ..  10 

B>j  j  (p,  I.)  -  u»d  P  (wd)L. 


(5.  53) 


where  P^'*  is  the  fluctuation  power  spectral  density  evaluated  at  the  azimuthal  drift  frequency 


.  Empirically  P^^(uj)  u>  1  where  the  value  of  r  most  often  is 

•  J 


=  -j-n —  =  - —  .  empirically  r  ^  co  where  the  value  < 

r  =  2  ±  1.  When  r  =  2,  one  obtains  the  very  simple  expression 

d(M)  =  k(M)  L  10  ( 

where  is  a  factor  dependent  on  the  fluctuation  magnitudes. 


(5.  54) 


18.  Schulz,  M.  ( 1975)  Geomagnetically  trapped  radiation.  Space  Sci.  Rev.  17:481. 
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l  or  gcoelectrie  lielu  tluctuations  it  can  be  shown  that 
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nhi’rc  1’  **''  is  'he  i!  :h  spatial  l  iiuni't'  component  of  the  electric  fid. I  iinvi  r 
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iisiti)ii,  I  or  ictatls,  see  Falthamma r,  anil  l >rm\  ul  1. 


rile  actual  calculation  of  L):  J  is  fuirlv  complicated  and  the  reader  is  referret 

11  .  2] 
to  the  research  literature.  However,  tor  typical  substorm  conditions,  Cornwall 

derived  the  simple  relation 


,dl->  .  K  1  b' 
J !  1  K 


Zir 


where  is  the  magnetic  moment  in  i\le\  G  and  Z.  is  the  particle  charge  state 
numbe  r. 

Both  and  are  dependent  on  the  geophysical  activity.  Likelv  values 

of  K*1'1*  fall  in  the  range  2  <  10  *°to2  <  10  (1. -shells)^  per  dav  and  mav  be 

-6-4  2 

found  in  the  range  10  to  10  (L-shells)  per  day.  The  total  radial  diffusion 
coefficient  is  then 

,,  ,.(M)  , dhl)  ,r  „„ 

1 '  1 . 1 .  ULL+ULL-  (5-°7 

It  should  be  emphasized  that  the  relations  in  fclqs.  (5.54)  and  (5.56)  represent 
simplified  conditions  that  may  not  always  be  realized.  If.  for  example,  p^1*  Qr  F 
do  not  follow  w>  ^  dependences,  then  L)|^  will  involve  dependence  on  ft,  and 
mav  have  a  different  form. 


5.3.5  Pitch  Angle  Diffusion 

The  presence  of  plasma  and  electromagnetic  waves  in  the  radiation  belts  im¬ 
plies  fast  low-amplitude  field  fluctuations,  and  some  of  these  waves  (such  as  the 
KLF  whistler  mode)  can  violate  the  first  adiabatic  invariant  ft.  These  fast 


19.  Falthammar,  C.G.  (1968)  Radial  diffusion  by  violation  of  the  third  adiabatic 

invariant,  in  Earth's  Particles  and  Fields  (Kd.  B.  M.  McCormac),  Reinhold 
Publishing,  New  York,  p.  157. 

20.  Cornwall,  J.A1.  (  1968)  Diffusion  processes  influenced  by  conjugation  point  wave 

phenomena.  Radio  Sci.  ^3^740. 

21.  Cornwall,  J.  Al.  (1972)  Radial  diffusion  of  ionized  helium  and  protons:  a  probe 

for  magnetospheric  dynamics,  J.  Geophys.  Res.  ^77: 1756. 


fluctuations  will,  in  principle,  also  violate  4-  and  J;  however,  for  these  adiabatic- 
invariants  the  effect  is  likely  to  be  at  least  in  part  averaged  out.  It  is  convenient 
to  convert  from  jl  to  equatorial  pitch  angle  (<y  1  coordinates 


2.3.2 
p  L  sin  a 
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wherebv  the  applicable  Jacobian  is 
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at  constant  particle  momentum.  In  the  pure  pitch  angle  diffusion  limit, 
unchanged  and  only  the  particle  direction  of  motion  changes. 

The  pure  pitch  angle  diffusion  equation  becomes 
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where  tS  and  are  the  source  ;ind  loss  functions  appropriate  for  the  pitch  angle 

diffusion  process.  D ^  ^  is  the  bounce  averaged  pure  pitch  angle  diffusion  eoeffi- 
o  o 

cient  related  to  the  local  pitch  angle  diffusion  coefficient  IJ  bv 
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1  he  actual  calculation  or  It  Iron  ibsorvc  i  .J  i-:-  a  v.  ave  di  st  ri  tuitions  in  die 
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magnetosphere  is  quite  complicated,  d  he  re  a  let  is  r.-p-rre  1  t  >  the  research  litera- 

■>•>  23  >.) 

ture  such  as  l.vons  id  al,  ’  "  Rettcrer  et  al.-  .m  l  reference.-  therein. 

In  general,  pitch  angle  diffusion  can  mne  alt-nit  b.-  violation  if  p  mlv  (that  is 
change  in  ),  bv  violation  of  J  onlv  (change  in  p,,  i.  or  bv  i  •  ombine  I  violation  of 
both  p  and  J.  The  pure  pitch  angle  diffusion  limp  is  i  leallv  reali-/e  1  when  both 
p  and  J  are  violated  so  that  the  ratio  >•'  p  p  .  hanges  while  T  remains  <  onslant 


22.  l.vons,  L.  11.  ,  Thorne.  H.  M.  ,  and  Kennel,  C.  1-',  (  T'7  1 )  Rlect ron  pitch  angle- 

diffusion  driven  bv  oblique  whistler  mode  turbulence,  .1.  Plasma.  Phvs.  8:58! 
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23.  l.vons,  1..R.,  Thorne.  K.M.,  and  Kennel,  C.  I  .  (1072)  Pitch  .angle  diffusion  of 

radiation  belt  electrons  within  the  plasruasphere,  .1.  Geophvs.  Res.  77:3155, 

24.  Retterer,  .I.M.,  .lasperse,  J.R..  ind  Chang,  1  .  S.  (1383)  A  new  approach  to 

pitch  angle  scattering  in  die  rn  ignetospliets-.  J.  (ieophvs.  Res.  88:201. 


When  this  happens,  there  is  tn>  energy  exchange  iietvvcrn  waves  and  parti,  les. 
Fhvsicallv,  this  implies  that  the  principal  interaction  is  between  the  partiele  and 
the  magnetic  field  of  the  wave. 


5.3.6  Energy  Diffusion 

Particles  can  become  energized  when  interacting  with  waves  (wave  damping) 
or  can  lose  energy  to  the  waves  (wave  instabi lit v).  One  mav  write  a  pure  energv 
diffusion  equation  in  the  form 
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(5.  621 


However,  such  an  equation  has  not  been  used  much  in  radiation  licit  physics,  sine 
almost  inevitably  <>  also  changes.  The  more  general  case  of  coupled  energy  and 
angular  diflusion  and  the  associated  diffusion  coefficients  needs  further  research. 
For  formalism  relating  to  such  coupled  processes,  see  Schulz  and  I.anzerotti. -i 


5.4  LOSS  MECHANISMS 


Energetic  particles  residing  iti  the  radiation  belts  are  subject  to  coIlisionaJ 
interactions  with  coexisting  particle  populations.  The  most  important  of  these  are 
'he  earth's  main  atmosphere,  the  atomic  hydrogen  exosphere  and  the  plasm  as  phe  re 
Near  the  earth  (that  is,  at  verv  low  1. -shells  or  for  small  equatorial  pitch  angles! 
such  collisions  constitute  a  dominant  energetic  particle  loss  mechanism.  Hut  even 
in  the  central  parts  of  the  radiation  belts.  Coulomb  collisions  and  charge  exchange 
.  an  he  quite  signilicant.  Filch  angle  scattering  of  particles  into  the  atmospheric 
bounce  loss  cone  (where  partiele -particle  collisions  are  dominant)  is  also  of  great 
significance,  especially  for  radiation  belt  electrons. 


5.4.1  Exosphere 

"1  he  terrestrial  exosphere  (or  geocorona)  is  a  continuation  of  the  atmosphere  to 
great  altitudes  where  collisions  are  infrequent  and  the  constituents  follow  ballistic 

•  trajectories.  The  principal  constituent  is  thought  to  be  atomic  hydrogen  with  a 

4  3  3  •>  •; 

ten-  i  tv  ranging  trim  -  10  atoms 'cm  at  10  km  altitude  to  ~  10“  atom  s  ■' cm '  at 

d  ‘  10  km  altitude.  I  able  o.  1  gives  the  mean  atomic  hvdrogen  number  density  |H| 

as  function  >!'  I. -shell  at  t he  equator*-  and  it  is  thought  to  be  an  averaim  ropresmta- 

M  >n  tot-  an  exospheric  temperature  of  ‘'dhoOK.  The  exospheric  temperatures  and 

densities  will  of  course  change  with  solar  and  £»eomai*net  if  activity. 

2;).  I  inslev,  H.  A.  (1M76)  Kvidencc  that  the  recovery  phase  rim*  current  consists  of 
helium  ions,  .1.  (Jeophvs.  Kcs.  8  1:61^2. 

- - - — - -  vA/v 


Table  5.  1.  Number  Density  of  Neutral  Hydrogen  at 
the  Equator  for  an  Exospheric  Temperature  of  950  K 
and  for  Average  Geomagnetic  Conditions""’ 


T  =  950  K 

L-Shell 

[H]  (1  cm3i 

L-Sholl 

[11]  ( 1  /cm3' 

1.  1 

15000. 

1 .  5 

16000. 

2.  0 

3700. 

2.  5 

1500. 

3.  0 

800. 

3.  5 

470. 

4.  0 

300. 

4.  5 

210. 

5.  0 

148. 

5.  5 

120. 

6.  0 

98. 

6,  5 

83. 

Also  overlapping  the  radiation  belts  is  the  terrestrial  plasm  as  phe  re  consisting 
of  thermal  ions  and  electrons  and  contained  within  an  I. -shell  range  roughlv  below 
about  L  -  4  to  6  bv  the  effect  of  the  corotational  electric  field  of  the  earth.  Figure  5.  1 
shows  the  electric  equipotential  lines  in  the  equatorial  plane  formed  bv  the  combina¬ 
tion  of  the  corotational  electric  field  and  the  solar-wind-induced  dawn -dusk 
"convection"  electric  field.  The  separatrix  between  closed  (around  the  earth)  and 
open  equipotential  lines  is  related  to  the  static  plasmapause,  although  in  a  dynamic 
time -variable  situation  no  simple  relation  between  the  two  exists.  Empirically, 
there  is  a  much  higher  densitv  of  "cold"  (  1000  -  10000  K)  plasmaspheric  particles 
below  the  plasmapause  than  beyond  it.  Based  on  data  deduced  from  ducted  V  L.K 
wave  propagation  experiments,  Cornwall^1  estimated  the  following  average  plasma - 
spheric  particle  densities: 
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cm 


for  L  >  L 
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(5.  63) 


where  L  =  4.  1  and  K  =  4. 64. 
o 

It  is  well  known  that  the  plasmasphere  deviates  frequently  and  strongly  from 
this  functional  form.  It  is  nevertheless  thought  that  this  formula  represents 
average  long-term  cold  plasma  densities  appropriate  for  long-term  steady  state 
radiation  belts  studies.  Further  improvements,  including  storm  and  substorm 
variability  will  be  needed  when  time -dependent  radiation  belt  models  are  developed. 
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5.4.2  Coulomb  Collisions 


Coulomb  collisions  arc  inelastic  interactions  between  charged  particles.  A 
radiation  belt  particle  colliding  with  an  exospheric  neutral  hydrogen  atom,  for 
example,  will  interact  with  the  internal  atomic  electric  field  whenever  the  impact 
parameter  is  less  than  the  atomic  radius,  or  with  the  electric  field  of  a  thermal 
(plasma)  proton  or  electron  out  to  the  particle's  Debve  shielding  distance.  The 
encounter  will  result  in  energy  transfer  from  the  energetic  (incident)  particle,  and 
in  dellection  (angular  scattering)  ot  both  particles.  Changes  in  ionization  states  of 
cither  or  both  particles  can  also  occur.  On  the  average,  angular  scattering  is 
important  tor  radiation  belt  electrons,  but  usually  not  so  important  for  the  much 
heavier  radiation  belt  ions.  The  eollisional  scattering  process  is  qualitatively 
illustrated  in  Figure  5.  13. 


Figure  o.  13.  Deflection  of  a  Positive  and  Negative 
Particle  in  the  Internal  Atomic  Electric  Field  of 
the  I  argot  Atom,  b  and  b~  are  impact  parameters 
and  t\\  is  the  atomic  radius.  If  the  impact 
parameter  is  greater  than  r^.  no  interaction  occurs. 

1'or  plasma  particles  r^  is  replaced  by  the  Debye 
shielding  distance  *  j 

Radiation  belt  particles  are  much  more  energetic  than  exospheric  particles, 
and  energy  is  transferred  from  the  energetic  ion  to  the  atomic  bound  electron(s) 
or  to  free  plasma  electrons.  I  he  changes  in  the  energy  spectrum  and  directional 
characteristics  of  particles  traversing  a  material  substance  whose  effective  thick¬ 
ness  is  greater  than  the  collision  mean  free  path  is  a  difficult  mathematical  problem. 
Great  simplification  is,  however,  obtained  when  the  differential  collision  cross 
sections  are  small  enough  to  overwhelmingly  favor  very  smalt  energy  losses  and 
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directional  changed  in  each  collision.  Fortunately,  tins  is  true  lor  Coulomb  colli¬ 
sions  in  which  the  classical  Rutherford  cross  section  is  valid.  See  Rossi  and 
26  27 

Olbert"  or  Jackson  for  details. 

One  finds  that  the  average  energy  loss  rate  for-  an  energetic  particle  passing 
through  a  gas  of  atomic  particles  is  given  by 

3P  =  "4  z  j  2^  [nl  >e  (nitie^  ','i  )  F( B,  7-t )  (a.  64) 


where  is  she  net  charge  state  number  of  the  incident  particle,  Z  is  the  nuclear- 

charge  number  of  the  target  gas  atoms,  m  is  the  electron  mass,  r  is  the  classical 

e  e 

electron  radius,  e  is  the  speed  of  light,  ,1  =  \  / e,  and  F(  is  a  slowly  varying 

function  determined  quantum -mechanically  bv 


F(0,Zt>  =  -2.6+  In  J7 2  ni2e^(  1-/J  2)*^2  l2(Zt> 


for  energetic  electrons, 


(5.  65) 


F(pj ,  Z{)  =  -2gi2+ln  4m2c^pi  1 -pi  2)2  I2(Zt>  for  energetic  ions,  (5.66) 

where  I(Z^)  -  13.  5  Z^  (in  eV)  is  the  approximate  ionization  potential  for  the  gas  atoms.2* 
l  et  f  =  f Ip.  ,  J,  4>;  t)  denote  the  distribution  function  for  equatorial  radiation 
belt  particles.  The  changes  due  to  stochastic  energy  loss  from  the  energetic  par¬ 
ticles  can  be  described  by 


5T  /Coulomb  =  "In  <<*fi/At>f)  - -J-j- (<A  J/At  >f) 


(5.  67) 


and  the  stochastic  time  averages  <Apt  /At>  and  <AJ/At>  are  then  simply  given  by 


< Ap  /At  >  = 


3pt  W  dC 


*E"M  dt 


(5.  68) 


-A  j  /a  t  >  =  ( )  (af-)- 


(5,  69) 


26.  Rossi,  B.  ,  and  Olbert,  S.  (1970)  Introduction  to  the  Physics  of  Space, 

McGraw-Hill  Book  Co.  ,  New  York. 

27.  Jackson,  J.  D.  (1975)  Classical  Electrodynamics,  2nd  edition,  John  Wiley 

and  Sons,  Inc,  ,  New  York. 


Cornwall^1  noted  that  it  is  possible  to  treat  the  slowly  varying  logarithmic 


i 

term  in  the  expression  for  as  an  approximate  constant.  Bv  including  contribu¬ 


tions  from  energetic  particle  plasma  electron  collisions  as  well  as  from  energetic 
particle  exospheric  atom  collisions,  one  derives  the  expression  valid  for  equatori- 


allv  mirroring  (a  -  7.12)  partii  les: 
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The  Coulomb  collision  loss  factor  is  given  bv 
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where  m..  1.  67  /  10  g  is  the  hvdrogen  atom  mass  and  e  (4.80288  :  0.00009) 


10 


10  psti  is  the  unit  charge.  is  the  ionic  charge  state  number  for  intrident 
■>  1  >f{  1 

energetic  ions.“  ’  “  For  a  fuller  treatment  of  the  collisional  process  the  reader 

9<i  vg 

is  referred  to  the  treatise  bv  Mott  and  Massev  or  Kossi  and  Olbert.  “ 


5.4.3  C  harge  Kxchange 

Collisi  mal  encounters  mav  also  change  the  ionization  state  of  the  colliding 
particles.  This  mav  email  ionization  of  the  target”  atom  and/or  alteration  of  the 
net  ionic  charge  of  'he  incident  particle.  The  latter  is  important  for  radiation  belt 
ions  since  ale  ost  all  of  *he  physical  processes  depend  directlv  on  their  charge 
state,  d  he  simples:  of  the  charge  exchange  reactions  is  the  one  that  neutralizes 
radiation  l»dt  protons  (fl  c 


fa.  72 ) 


where  underlining  den  >tes  energeth  particles.  The  above  reaition  is  an  example 
>f  the  electron  capture  .»r  1  cess  in  which  the  in-  idem  proton  picks  up  (or  captures) 
the  irhital  elect  r  »n  I  rom  the  "hernial  hvdrogen  atom,  which  then  becomes  a  low 
energy  proton.  Once  neutralized,  1  lie  incident  proton  (now  fast  neutral  hvdrogen) 
is  no  longer  subject  to  the  magnetic  deflecting  lorce  and  escapes  from  tin-  trapping 
region.  Mat  roseopiraltv  these  events  occur  at  random,  and  fast  neutral  atoms  thus 
exit  'he  radiation  belts  in  all  directions.  Since  the  speed  of  these  particles  greatly 

28.  Spjeldvik,  W.  N.  (  1977)  Kquilibrium  struc  ture  of  equatorial  mirroring  radiation 
belt  protons,  J.  Geophvs.  lies.  82:2801. 

29.  Molt,  \.  t  ,  ,  arul  Massev,  H.S.W.  (1952)  d  he  Theorv  of  Atomic  Collisions, 

Clarendon  Press,  Oxford,  I  nitcd  Kingdom. 


exceeds  the  earth's  gravitational  escape  speed,  manv  disappear  to  outer  space. 

Some  of  these  fast  neutral  atoms  move  towards  the  earth  where  they  produce 
secondary  interactions  upon  entering  the  atmosphere. 

In  a  similar  manner,  the  other  singlv  charged  ions  can  be  neutralized  in  such 
collisions: 

He*  •  H  -*  He  •  h' 

C*  +  H  ->  C  *  H*  (5.75) 

Ot  *  H  -*  O  •  H'  . 

These  newly  generated  thermal  hvdrogen  ions  (protons)  make  a  contribution  to  the 
earth's  plasma  envelope;  however,  their  rate  of  formation  is  probably  less  titan  the 
rate  at  which  such  ions  are  supplied  from  the  topside  ionosphere.  One  may  note 
that  the  total  charge  is  conserved  in  the  radiation  belts  under  the  charge  exchange 
reactions. 

There  is  also  another  class  of  ion  charge  exchange;  this  results  from  electron 
stripping  reactions  by  which  the  energetic  tons  (necessarily  heavier  than  hvdrogen) 
lose  one  or  more  of  the  remaining  bound  electrons.  For  example,  there  is  a  finite 
probability  that  an  ion  (sav  O  )  in  a  collision  with  a  thermal  hvdrogen  atom  loses 
several  of  its  electrons: 

O  ■  H  -*  O5  ■  H  •  4e'  .  (5.  74) 

In  this  example  four  electrons  are  stripped  off.  Whether  the  fast  O  ion  actually 
loses  one  or  more  of  its  electrons  or  captures  one  from  the  hvdrogen  atom  strongly 
depends  on  the  energy  of  the  incident  ion  and  the  details  of  the  collision.  Since  this 
type  of  reaction  preserves  a  non-zero  ionic  state,  the  ion  remains  trapped  within 
the  radiation  belts.  For  this  reason  one  distinguishes  between  internal  charge 
exchange  reactions  such  as  the  stripping  reactions: 


o’- 

ii 

-*  O'* ' 

•  H 

•  c 

o3* 

•  H 

-  o7' 

■  H 

•  4c 

(5.  75) 

etc. 

with  their  reverse  reactions  such  as 


02+  +  H  "*  0+  +  H  (5.  76) 

on  one  hand,  and  reactions  that  produce  immediate  particle  loss  (lowest  charge 

state  to  neutral)  noted  above  on  the  other.  Notice  also  that  to  become  neutralized, 
8+ 

an  O  ion  requires  a  minimum  of  eight  separate  collisions  with  hydrogen  atoms 
(since  only  one  electron  may  be  captured  in  each  collision).  This  should  be  con¬ 
trasted  with  electron  stripping  reactions  in  which  the  multiple  charge  state  changes 
occur  in  a  single  collisional  encounter.  The  principal  features  of  the  charge  ex¬ 
change  chemistry  are  illustrated  in  Figure  5.  14. 


Figure  5.  14.  Flow  Diagram  for  Energetic  Oxvgen  Ion  Charge  Exchange.  This 
figure  illustrates  that  while  multiple  charge  state  increases  are  possible  in  a 
single  encounter,  only  single  step  charge  state  reductions  can  occur  in  a 
hydrogen  atom  gas.  Oxygen  ions  in  lower  charge  states  may  be  further 
ionized  (through  single  or  multiple  electron  stripping)  in  single  collisional 
encounters  with  the  exospheric  gas  atoms.  The  ions  may  lose  net  charge  in 
collisions  by  capturing  the  bound  electron  from  a  thermal  hydrogen  atom. 

To  become  neutralized,  a  fully  ionized  oxygen  ion  requires  a  minimum  of 
eight  separate  collisions 

The  probability  that  a  given  charge  exchange  process  actually  takes  place  in  a 
collision  is  expressed  as  a  reaction  cross  section.  Such  cross  sections  derive 
from  laboratory  studies  of  collision  processes;  however,  for  manv  ions  the  per¬ 
tinent  cross  sections  have  not  vet  been  measured  over  a  sufficiently  large  range  of 
particle  energies.  Figure  5.  15  shows  the  cross  sections  for  the  charge  exchange  of 
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protons  incident  on  atomic  hvdrogt-n  from  a  compilation  bv  Spjeldvik."  Notice- 
that  the  charge-  exchange  cross  se-ction  is  quite  high  (ajQ  >  10  cm")  be-low  about 
20  keV,  but  falls  off  sharply  with  higher  energies. 


Figure  5.  15.  Proton  Charge  Exchange  Cross  Section  as  Function 
of  Proton  Energy.  It  is  worth  noting  that  this  cross  section  is 

p  .  15  9 

high  (CT-.g^  10  cm")  helow  a  few  tens  of  keV,  causing  short 
proton  lu  lifetimes  at  those  energies,  and  it  falls  off  rapidly 
towards  higher  energies.  The  values  between  0.  4  and  1  keV  are 
from  the  experimental  work  of  Kite  et  ala®,  the  data  from  1  keV 
are  due  to  a  compilation  by  Ctaflin^1  and  above  1000  keV  the 
theoretical  results  from  Brinkmann  and  Kramers'^ 


For  ions  heavier  than  protons,  multiple  potential  charge  states  are  available. 
F’or  helium  ions,  one  must  consider  not  only  the  cross  section  for  the  neutralization 
reaction  but  charge  state  changes:  state  1  state  2  and  state  2  ■*  state  1  as  well. 
Thus,  for  helium  there  are  three  important  cross  sections  to  be  included.  A  point 
worth  noting  is  that  the  process  transforming  He"  to  He^  outweighs  the  charge 
state  reducing  reactions  in  the  high  energy  part  of  the  radiation  belts,  essentially 
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above  800  keV.  This  is  explicitly  depi cted  in  Figure  IB,  which  sh  nvs  the  three 
■barge  exchange  cross  sections  tor  radiation  belt  helium  ions. 


Figure  5.  16.  Helium  Ion  Charge  Exchange  Cross  Section  as 
l- unction  of  Energy.  Solid  lines  show  the  values  compiled  bv 
ClafUm’1  from  many  experiments,  and  the  dotted  lines  are 
extrapolations.  The  dashed  lines  are  computed  from  the 
theoretical  work  of  Brinkmann  and  Kramers.  The  three 
pertinent  charge  exchange  cross  sections  are; 

He  — ►He  (ion  neutralization),  He2* — ►He*  and  He2* ►He+. 

The  extrapolated  curve  above  1000  keY  (for  He* — ►  He2*)  is 
progressively  unreliable  towards  higher  energies'^ 


Charge  exchange  cross  sections  for  the  heavier  ions  such  as  carbon  and 
ixvgen  have  been  measured  over  a  small  fraction  of  the  radiation  belt  energy  range. 
1  he  situation  is  particularly  severe  for  carbon  ions  where  almost  no  measurements 
have  been  made.  For  oxvgen  ions  there  exist  a  number  of  measurements,  but 
unfortunately  the  laboratory  work  used  particles  other  than  atomic  hydrogen  as 
targets,  figure  5.  17  shows  estimated  ion  charge  loss  (electron  capture)  cross 
sections  based  on  a  compilation  of  a  number  of  measurements  using  H,  air  and  O  as 

■id.  Spieldvik,  W.X.,  and  Fritz,  T.  A.  (11)78)  Energetic  ionized  helium  in  the  quiet 
tinie  radiation  belts:  Theory  and  comparison  with  observation,  J.  Geophys. 
lies,  83:654. 
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target  particles.  Estimates  of  the  charge  gain  (electron  loss)  cross  sections 
have  been  made  from  limited  data  available.  Examples  are  illustrated  in 
Figure  5.  18.  For  the  remaining  reactions  between  the  multiple  charge  states 

no  laboratory  measurements  have  been  reported,  and  one  must  use  crude 

,  34 

estimates. 


Figure  5.  17.  Charge  Exchange  Cross  Sections  for  (Atomic) 
Oxygen  Ions  in  an  Atomic  Hydrogen  Gas  (Charge  l.oss/Electron 

Capture  Cross  Sections  o'  — *-  i>+  Estimated  From 

Sparse  Available  Data)  Stebbings  et  al.  ^5  give  data  with 
M-targets  up  to  8  keV,  Lo  ej  aP6  give  data  with  O  targets, 
and  MacDonald  and  Martin^  give  data  with.  He  targets.  Also 
shown  are  results  with  rarefied  air  targets^.  The  cross 
sections  for  reactions  with  i>  1  have  not  been  measured 
below  ~  10  MeV^4 


The  mathematical  description  of  the  effects  of  the  charge  exchange  processes 
on  radiation  belt  ion  distribution  functions  is  through  loss  and  gain  terms.  For 
radiation  belt  protons  one  has  the  expression 


Charge  Exchange 


'A10f 


(5.77) 


(Due  to  tne  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  127.) 
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l-  igure  5.  18.  Charge  Exchange  Cross  -Sections  for  (Atomic) 
Oxygen  Ions  in  an  Atomic  Hydrogen  Gas  (Charge  Gain/ Electron 

Loss  Cross  Sections  ()+ — ►0^+^+  +  j  e  .  )  The  data  shown 
are  from  Lo  et  ar'  for  O  targets.  The  solid  lines  are  crude 
estimates  made  by  Spjeldvik  and  I-'ritz.  More  reliable 
cross  sections  for  oxygen  ions  and  other  heavy  ions  need  to 
be  established 


since  the  loss  rate  is  proportional  to  the  number  of  protons  present  (or  more  pre¬ 
cisely,  the  distribution  function  f),  and  also  proportional  to  the  factor 


=  <  vCTj  |  H]  >  where  v  is  the  ion  speed  and  cr.  the  charge  exchange  cross  section 
for  transformation  state  i  **  state  j.  To  describe  the  heavier  ions  where  more  than 


one  charge  state  applies,  one  considers  a  distribution  function  for  each  charge  state. 
For  helium  ions,  let  f^  and  f^  be  the  distribution  functions  for  lle+  and  He^+  re¬ 


spectively.  The  charge  exchange  processes  arc  then  represented  by 
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and  this  provides  a  coupling  between  the  distribution  functions. 

In  general,  for  an  ion  species  with  s  available  charge  states  one  obtains  charge 
exchange  expressions  of  the  form: 


Charge  Exchange 
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and  when  the  exosphere  predominantly  consists  of  atomic  hydrogen  (as  in  the  case 
of  the  earth)  the  first  summation  contains  contributions  from  j  -  i  -  1  only. 


5.4.4  Wave-Particle  Interactions 

Plasma  waves  play  an  important  role  in  radiation  belt  physics.  Different  con¬ 
ditions  for  wave  propagation,  growth,  and  decay  in  the  magnetosphere  have  led  to 

delineation  into  numerous  wave  modes,  l  or  a  detailed  classification  of  these  modes 
39  40 

see  Stix  or  Kennel  et  al.  In  a  previous  section,  the  radiation  belt  transport 
equation  was  given;  here  we  shall  outline  some  of  the  concepts  concerning  the  inter¬ 
actions. 

Angular  scattering  can  result  from  interactions  with  electromagnetic  waves. 

As  an  example,  consider  gyro  (or  cyclotron)  resonance:  A  radiation  belt  particle 
spiraling  around  a  magnetic  line  of  force  (the  guiding  center  locus  field  line)  will 
gyrate  at  a  rate  determined  by  the  magnetic  field  strength,  mass,  and  charge  of  the 
particle.  An  electromagnetic  wave  (with  electric  and  magnetic  wave  vectors)  also 
propagating  along  that  field  will  rotate  according  to  its  wave  frequency.  When  both 
the  sense  of  rotation  and  the  rotation  frequency  match  for  both  wave  and  particle, 
the  particle  will  be  subject  to  an  essentially  constant  "wave"  field  for  the  duration 
of  the  encounter.  The  particle  can  exchange  energy  with  the  wave  through  the  elec¬ 
trical  interaction  and/or  deflection  can  occur  through  the  magnetic  interaction. 
Higher  order  resonances  can  also  take  place,  for  example,  if  the  rotation  rates 
differ  by  a  factor  of  two.  When  the  particle-wave  interaction  is  primarily  via  the 
magnetic  wave  vector  (as  in  the  case  of  the  radiation  belt  electron  -  whistler  mode 
hiss  wave  interaction)  the  result  is  primarily  angular  scattering;  this  is  *he  pure 
pitch  angle  scattering  limit.  Other  resonances  including  bounce  resonance  and 
drift  resonance  can  also  be  important. 

Plasma  waves  in  the  radiation  belts  may  remain  in  the  area  where  they  were 

generated  (very  small  group  velocity)  or  propagate  afar  (large  group  velocity).  They 

have  frequency  and  wave  length,  and  there  are  different  polarization  properties. 

3() 

Propagation  properties  are  determined  by  a  dispersion  relation  "  that  in  part 
depends  on  the  density  of  the  plasma  in  which  the  wave  exists  and  on  the  geomagnetic 
field.  Stability  or  instability  of  the  waves  is  frequently  determined  bv  the  energetic 
particles  with  which  the  waves  can  interact. 


39.  Stix,  T.  11.  (1962)  The  Theory  of  Plasma  \\  a ves,  McGraw-Hill  Book  Go.  ,  New  Yoi 
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l’lasma  waves  are  said  to  be  unstable  to  growth  if  interactions  with  the  charged 
particles  transfer  energy  to  the  waves;  if  the  transfer  is  from  the  waves  to  the  par¬ 
ticles  the  waves  are  said  to  be  damped.  If  plasma  waves  interact  with  radiation 
belt  particles  in  such  a  way  as  to  have  their  principal  interactions  locally,  then  it 
is  possible  to  self-consistently  give  a  theoretical  treatment  of  both  wave  and  particle 
properties.  On  the  other  hand,  if  the  waves  have  significant  spatial  propagation  so 
that  they,  for  example,  gain  energy  from  one  particle  population  and  propagate  to 
interact  further  with  another  particle  population  elsewhere,  then  the  latter  process 
is  termed  "parasitic  ".  Both  types  of  processes  are  of  major  significance  in  the 
radiation  belts. 

5.4.4.  1  PITCH  ANGLE  SCATTERING  INTO  THE  LOSS  CONE 

Interactions  between  energetic  particles  and  plasma  waves  can  significantly 
affect  the  trapped  particle  population's  energy  and  angular  distribution.  If  the  pitch 
angle  is  altered  so  that  the  particle  finds  itself  within  the  atmospheric  bounce  loss 
cone,  it  will  have  a  high  probability  of  becoming  lost  upon  entering  the  atmosphere, 
lor  ions  this  probability  is  virtually  IOC  percent  while  for  energetic  electrons  a 
certain  fraction  of  the  precipitated  particles  are  backscattered  up  into  the  mag¬ 
netosphere  only  to  encounter  the  conjugate  hemisphere  during  the  succeeding  bounce 
motion. 

The  directional  change  in  the  angular  scattering  process  due  to  plasma  wave  or 
collisional  interactions  is  random:  to  lowest  order  it  is  just  as  likely  for  a  single 
particle  to  be  scattered  from  lower  to  higher  pitch  angles  as  in  the  opposite  direc¬ 
tion.  However,  w'hen  the  particles  are  anisotropically  distributed  in  pitch  angle 

(for  example,  with  g— —  >  0),  the  number  of  particles  scattered  from  the  higher 
a  o 

particle  density  region  towards  the  lower  density  region  is  greater  than  the  number 
scattered  the  other  way.  As  a  consequence,  the  stochastic  process  is  biased  by  the 
particle  pitch  angle  distribution  and  a  net  ditfusion  flux  occurs.  This  pitch  angle 
diffusion  flux  is  just 


)  „  sin  2 a  T (a  ) 

a  a  o  o  9  a 

o  o  o 


(5.  81) 


so  that  the  diffusive  transport  term  in  Eq.  (5.60)  may  be  written  as  simply 


(5.  82) 


Atmospheric  particle  losses  within  the  bounce  loss  cone  o  (j  generally  cause 
f(oro  =  1°  be  near  zero  (except  under  conditions  of  extremely  strong  scattering). 
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This  gencr.dlv  f.ivtii-s  a  positive  pitch  angle  anisotropy  (loss  c one  distribution),  and 
particles  can  he  lost  t<>  tin-  atmosphere  at  am  1. -shell  due  to  the  pitch  angle  scatter¬ 
ing  process.  1  »f  cieirse,  of  o  ,y  may  be  negative  under  certain  conditions,  and 
•bus  net  revers.-  Ir'iisum  can  take  nlace.  Lxamples  are  an  at r" >s pile ri c  source  cone 
in  the  aurora!  /one  s  .  consequence  of  auroral  electric  fields,  effects  of  I. -shell 
splitting  due  to  tievi  .tions  from  the  dipolar  azimuthal  symmetry  of  the  magnetic 
field,  *  or  (i  .rtich  mic  tions  during  disturbed  times. 

Analv.->i.->  of  resonant  Aave-pn  rticle  interactions  violating  the  first  adiabatic 
invariant  have  a.  ner  .llv  considered  waves  with  frequencies  near  the  gy  rofrequency 
of  tile  energetic  parti,  les.  1'tie  condition  for  resonance  between  waves  of  angular 
frequency  w  and  particles  at  the  gv  rofrequency  12  is  given  by 


‘  kU  VH 
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where  k||  and  vm  are  tile  parallel  (to  B)  nave  vector  and  particle  velocity  respective 
ly.  l  or  a  given  wave  mode  there  is  a  dispersion  relation  linking,  ut  and  k.  The 
evelotron  harmonic  resonances  have  the  (Doppler  shifted)  wave  frequency  equal  to 
a  harmonic  (n  =  ±  1,  n  =  ±  2,  n  =  ±  3,  . .  .  )  of  the  particle  gyrofrequency,  and  the 
classical  i.andau  resonance  (n  =  0)  has  the  wave  parallel  phase  velocity  v^  =  a>/k||. 
The  Landau  resonance  results  in  diffusion  solely  in  v  ,  conserving  v.  .  Thus 

M 

the  pitch  angle  change  Aa  ut  a  given  location  along  the  particle  trajectory  is  related 
to  the  parallel  velocity  change  AV|(  bv 

2 

.  sin  of  .  / r  o 4  \ 

Aa  ; Av.  .  (5.  84) 


The  cyclotron  harmonic  resonance  can  likewise  produce  pitch  angle  changes 
resulting  from  the  energy  exchange,  or  the  interaction  can  be  primarily  with  the 
magnetic  wave  vector  causing  pitch  angle  scattering  more  directly  with  little  energy 
exchange.  The  latter  process  is  particularly  important  for  electrons. 

In  general,  the  study  of  wave-particle  interactions  requires  an  extensive 
mathematical  treatment,  and  the  necessary  derivations  and  analysis  are  beyond  the 
scope  of  this  report.  The  interested  reader  is  referred  to  a  number  of  works  on 
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the  subject  |Stix.  “  Sagdeev  and  Galeev,  Lerche,  Kennel  and  Engelman, 

‘->2  ‘>3  24 

Lyons  et  al“  ’  “  and  Retterer  et  al  ]. 

5. 4. 4.  2  SGATTE  RING  OF  ENERGETIC  ELECTRONS 

Fitch  angle  scattering  is  particularly  important  for  radiation  belt  electrons. 

22  23 

Here,  we  bypass  the  extensive  mathematical  treatment  found  in  Lyons  et  al. 

The  different  resonances  are  illustrated  in  Figure  5.  19,  which  depicts  the  regions 
of  cyclotron  resonance  in  velocity  space  resulting  from  waves  distributed  over  a 
band  of  parallel  wave  vectors  Ak(..  Resonance  at  each  cyclotron  harmonic  occurs 
over  a  band  of  parallel  velocities  Av(|  ;  for  simplicity,  relativistic  effects  for 
electrons  >  500  keV  are  not  included.  No  cyclotron  harmonic  interaction  occurs  for 

v  <  v,.  .  ,  and  the  value  of  v,.  depends  on  the  actual  k.,  band  over  which  the 

II  nun  II  nun  r  11 

wave  energy  is  distributed.  This  is  locally  true  where  the  geomagnetic  field  may 
be  considered  uniform  on  the  scale  of  the  gvroradius. 

The  earth's  magnetic  field  is,  however,  quite  inhomogeneous  when  considered 
in  its  totality,  and  the  wave  energy  is,  in  general,  unevenly  distributed  over  the 
space  of  the  radiation  belts.  During  the  particle  bounce  motion,  as  the  particle 
moves  away  from  the  geomagnetic  equator  along  its  trajectory,  the  increasing  mag¬ 
netic  field  strength  causes  both  particle  pitch  angle  and  parallel  velocity  to  change. 
The  increase  in  the  local  pitch  angle  as  a  particle  moves  away  from  the  equatorial 
plane  (see  Eq.  (5.  10 ) )  implies  that  cyclotron  resonance  can  occur  for  all  equatorial 
parallel  particle  energies  greater  than  a  minimum  value  E’u  .  . 

Figure  5.  20  illustrates  radiation  belt  electron  and  ELF  whistler  mode  wave 
propagation.  The  significant  wave-particle  interactions  for  energetic  electrons 
are  shown.  Note  that  the  wave  propagation  does  not  necessarily  follow  magnetic 
field  lines  and,  therefore,  the  particles  may  interact  with  waves  generated  over  a 
significant  volume  of  the  magnetosphere.  For  comparison,  the  inner  and  outer 
radiation  zones  (for  electrons)  are  indicated. 

Satellite  measurements  have  shown  that  a  band  of  whistler  mode  waves  centered 
around  a  few  hundred  Hertz  exists  essentially  continually  within  the  plasmasphere. 
Because  of  the  persistence  and  audio  frequency  range,  this  type  of  wave  is  called 
hiss.  Ray  tracing  studies  have  shown  that  the  plasmaspheric  hiss  can  readily 
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Figure  5. 21.  Bounce  Averaged  Klectron 

Pitcli  Angle  Diffusion  Coefficient  Da  a 
a  o  o 

Calculated  for  All  Cyclotron  Harmonic 
Resonances  and  the  l.andau  Resonance. 
Note  that  at  each  energy  there  is  a 
region  of  very  low  pitch  angle  diffusion 
coefficient  (for  20  keV  electrons,  near 
87°);  this  "bottleneck"  in  the  pitch 
angle  diffusion  coefficient  is  the  cause 
of  the  "bumps"  in  the  actual  electron 
pitch  angle  distributions  giving  rise 
to  the  so-called  bell-shaped 
distributions'^ 


Figure  5.  22.  Radiation  Belt  Electron 
Precipitation  Lifetimes  Due  to  Electron 
Pitch  Angle  Scattering  Into  the  Atmos¬ 
pheric  Bounce  Loss  Cone.  These  life¬ 
times  are  given  for  an  average  whistler 
wave  amplitude  of  35  my  .  Notice  that 
there  is  generally  a  minimum  in  these 
lifetimes  at  each  energy,  and  that  this 
minimum  is  displaced  towards  lower 
L-shells  with  higher  electron  energy. 

At  2000  keV  this  minimum  occurs  at 
I.  -  2.  7  while  at  500  keV  it  is  sub¬ 
stantially  broader  and  occurs  at 
L—  4.  2.  This  minimum  electron 
lifetime  is  an  important  determinant 
of  the  electron  "slot"  location, 
defining  the  separation  between  the 
two  radiation  belts 


few  principal  aspects  of  the  radiation  belt  saturation  process  are  outlined;  for  a 
more  detailed  mathematical  treatment  the  reader  is  referred  to  Kennel  and  Petschek^ 
or  Schulz. 
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It  has  been  demonstrated  that  there  is  a  limiting  value  I'  that  the  radiation  belt 

integral  omni  -directional  flux  l  cannot  exceed  without  provoking  a  cyclotron  wave 

instability.  Schulz  and  Lanzerotti"*  estimated  the  equatorial  flux  limit  to  be  about 

['  -  10^  L  ^  particles/ cm'^  sec  integrated  over  all  energies  and  pitch  angles.  The 

linear  wave  growth  rate  is  denoted  y  which  is  proportional  to  the  ratio  I/I'  and 

therefore,  the  growth  rate  for  wave  energy  is  2  }  .  An  incipient  wave  undergoes 

a  partial  reflection  (reflection  coefficient  R  S  1)  upon  traveling  a  distance 

d  -  l.  R  ;  the  remaining  fraction,  1-R,  of  the  wave  is  lost  from  the  radiation  belts. 

The  time  interval  between  wave  reflection  is  T  -  1.  R.,  v  where  v  da’  dk  is 

!'•  g  g 

ihe  group  w av e  velocity.  The  condition  for  marginal  stability  is  that  the  waves  on 
the  average  do  not  grow  further  in  time;  that  is,  the  decrease  in  the  wave  amplitude 
upon  reflection  is  restored  in  one  traversal  between  bounces. 


R  exp  [  >■)  1.  R  / v  ] 

«  !'•  K 


(5.  85) 


This  defines  the  marginal  growth  rate  to  be 


<v(j'2I.  RK)  In  R  . 


(5.  86) 


If  I  exceeds  1'  (that  is,  if  )•  exceeds  7  '  )  the  consequence  is  a  net  growth  of  wave 
energy,  and  the  stronger  waves  scatter  the  particles  in  pitch  angle  resulting  in 
precipitation  into  the  atmosphere  and  reduction  in  the  particle  flux  until  I  no  longer 
exceeds  I'.  If  an  external  particle  injection  source  is  strong  enough  to  more  than 
offset  this  maximum  loss  rate  then  the  radiation  belt  particle  fluxes  will  increase 
bevond  the  stabtv  trapped  flux  limit  I;  this  may  indeed  happen  for  impulsive  sources 
associated  with  disturbed  conditions.  The  aftermath  of  such  an  injection  will,  how¬ 
ever,  be  a  rapid  reduction  to  the  stably  trapped  flux  intensity  level.  On  the  other 
hand,  if  [  is  below  V,  particle  losses  due  to  the  waves  will  diminish  due  to  the 
lower  wave  intensity. 

Figure  5.23  shows  a  comparison  of  the  theoretical  radiation  belt  flux  limit  with 

46 

data.  Notice  that  most  of  the  time  the  outer  radiation  zone  particle  fluxes  are 
just  below  the  theoretical  characteristic  flux  limit. 
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LOCAL  TIME 


Figure  5.  23.  Comparison  of  the  Predicted  Whistler  Self-Excitation 
Limit  With  E  >  40  keV  Electron  Data  Recorded  With  the  Explorer  14 
Spacecraft^® 


5.  4.  4.  4  STRONG  DIFFUSION  LIMIT 


There  is,  however,  another  important  limit  to  consider,  namely  the  upper 
limit  on  the  rate  at  which  radiation  belt  particles  can  become  lost  to  the  atmosphere 
by  scattering  into  the  bounce  loss  cone.  Consider  the  effect  of  very  rapid  pitch 
angle  scattering.  The  particle  pitch  angle  distribution  will  then  become  essentially 
isotropic;  that  is,  the  flux  per  differential  solid  angle  interval  will  be  the  same  at 
all  pitch  angles,  even  within  the  loss  cone.  Within  the  loss  cone  (a  5  °0LC^ 
particles  will  be  lost  from  the  radiation  belts  on  a  time  scale  T^,  and  for  the  strong 
dif.usion  limit  to  be  reached,  the  scattering  mechanism  must  be  faster  than  this. 

The  solid  angle  of  the  Loss  cone  is  just  Q  j  ^  =  2~  (1-cos  and  the  maximum 

particle  flux  loss  rate  6  is  then 


'b 


The  only  exception  to  this  concept  occurs  when  the  loss  cone  becomes  an  atmos¬ 
pheric  source  cone  due  to  the  auroral  electric  field  parallel  to  the  magnetic  field 
and  a  source  cone  distribution  exiting  one  hemisphere  can  become  lost  through 
direct  precipitation  into  the  conjugate  hemisphere. 

5.5  THEORETICAL  RADIATION  BELT  MODELS 

Two  types  of  radiation  belt  models  now  exist:  empirical  models  based  on  com¬ 
piled  observations  and  theoretical  models  derived  from  our  knowledge  of  radiation 
belt  physics.  The  ideal  empirical  models  represent  the  standards  with  which  the 
theoretical  models  must  be  compared,  and  they  are  also  useful  for  engineering  pur¬ 
poses.  However,  thev  give  only  a  smoothed  statistical  picture  of  the  time  period 
when  the  data  were  obtained,  and  the  data  collection  process  is  subject  to  experi¬ 
mental  errors  and  misinterpretations.  Theoretical  models,  on  the  other  hand,  ran 
be  used  to  simulate  and  predict  radiation  belt  behavior  under  a  variety  of  conditions; 
quiet  times,  disturbed  times,  or  magnetic:  storms.  They  can  be  used  to  studv  ionic 
species  and  charge  state  distributions  for  which  no  experimental  information  is 
available.  Such  models  are  only  as  valid  as  the  physical  approximation  thev  are 
based  on,  however.  In  the  following  sections  simple  theoretical  models  are  outlined. 


5.5.1  Quiet  Time,  Steady  State  Models 


The  study  of  quiet  geomagnetic  conditions  using  theoretical  radiation  belt  models 
serves  as  an  important  test  of  our  current  understanding  of  the  trapped  radiation 
environment.  It  is  assumed  that  under  quiet  time  conditions  the  radiation  belts  can 
be  described  in  terms  of  an  equilibrium  balance  between  sources,  internal  trans¬ 
port,  and  losses.  Symbolically  one  writes 

3  f  _/9f\  /  9  f\  /9f 

Ft  ”  \FT  /source  +  \FT /transport  +  \F"t 

The  sources  are  considered  to  be  located  on  the  boundary  of  the  trapping  region 
(and  thus  descrtbable  via  a  boundary  condition  on  the  trapped  flux  for  each  species) 
as  well  as  in  the  interior  of  the  trapping  region  (such  as  the  CHAND-source).  The 
transport  mechanisms  are  radial  diffusion  and  pitch  angle  diffusion.  Particle  losses 
are  caused  by  charge  exchange  ion  neutralization,  energy  degradation,  and  pitch 
angle  diffusion  into  the  bounce  loss  cone.  Particles  diffusing  outward  may  encounter 
the  magnetopause  and  become  lost  from  the  magnetosphere.  Radiation  belt  theoreti¬ 
cal  models  constitute  the  combination  of  these  processes. 

5.5.  1.  1  FORMULATION  OF  BOUNDARY  CONDITIONS 

For  steady  state  conditions  one  needs  in  general  boundary  conditions  on  the 
particle  distribution  function  f: 

(1)  At  L  =  L  (outer  boundary)  one  specifies  the  particles  energy  spectra 

max 

and  angular  distributions,  which  defines  f  at  the  outer  boundary. 

(2)  At  L  =  L  .  =1  (inner  boundary  formed  bv  the  surface  of  the  earth)  the 

min 

particle  distribution  function  vanishes  (f  =  0). 

(3)  At  a  =  7T / 2:  3f/3  a  =  0. 

o  o 

(4)  At  aQ  =  0  or  n:  dil'daQ  =  0. 

(5)  At  p  =  p  (or  E  =  E  ):  f  =  0  where  E  ~  100  MeY. 

rmax  max  max 

(6)  By  symmetry  f(aQ)  -  f(ii  -  a  )  so  that  one  needs  to  consider  only  the 
range  0  <  aQ  <  77  /  2 . 

The  complete  radiation  belt  modeling  can  be  simplified  in  a  way  that  retains 
much  of  the  essential  physics  but  greatly  simplifies  the  mathematics,  namely, 
separating  radial  diffusion  from  pitch  angle  diffusion.  One  treats  radial  diffusion 
for  aQ  =17/ 2  particles  only,  and  with  given  f(ao  =  tt/ 2)  the  pitch  angle  diffusion  pro¬ 
cess  at  fixed  I. -shells  only.  This  approach  is  not  strictly  correct  when  there  is 
significant  interaction  between  the  pitch  angle  and  radial  diffusion  modes,  but  it 
reduces  an  almost  unmanageable  numerical  problem  to  cases  where  solutions  can 
be  found.  The  results  should  therefore  be  treated  with  some  caution,  especially  in 
regions  where  the  time  scales  of  the  two  processes  are  similar. 


/loss 
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f>.  5.  1.  2  ELECTRON  iU(  >DEL 


Radiation  belt  electrons  are  described  by  the  diffusion  equation 


d_f 
o  t 


o 


1  a 

.sin  2  a  T( a  )  a a 
o  o  o 


a  cv 
o  o 


sin  2(7  TL7  ) 
u  o 


of _ 

a  a 

o 


-  G  (I.)m 
e 


<5.  89) 


where  G  (I.)  is  the  Coulomb  loss  factor  for  electrons. 
v  48 

I.yons  and  Thorne  separated  radial  and  pitch  angle  diffusion  by  substituting  for 
Hq.  (5.  89)  the  'no  equations 


o_f 
o  t 


2  o 


FT. 


and 


a  I' 
at 


1  o 

sir.  2  a  T( (7  1  d(7 

o  o  o 


D  sin  2  a  T(  a  ) 

a  a  o  o 

o  o 


(5.  90) 


(5. 91) 


They  solved  Eq.  (5.  91)  under  an  assumed  angular  distribution  preserving  exponential 

decay  conditions  and  obtained  the  lowest  normal  mode  pitch  angle  distribution  and 

electron  precipitation  lifetimes  r  ,.  Their  results  are  shown  in  Figures  5.  24  and 

w 

5.  25  respectively.  The  precipitation  lifetimes  were  calculated  assuming  a  mean 

ELF  whistler  mode  wave  amplitude  of  14  =  35  nT  ( 1  nt  =  1  >  =  10  3  G).  For  other 

wo  _2 

wave  amplitudes  these  lifetime  scale  as  (B  /  B  ) 

w  w  o 

Using  these  characteristic  electron  scattering  precipitation  lifetimes  it  is 
possible  to  solve  the  equation  for  steady  state  radial  diffusion  [Eq.  (5.  90)]  for 
equatorially  mirroring  radiation  belt  electrons.  The  results  are  shown  in 
Figure  5.  25. 

A  sin.plification  in  this  work  was  a  priori  approximation  of  the  Coulomb  collision 
"loss  time"'  defined  as  T  cc  =  £(dC/dt)  where  d£  dt  is  given  by  Eq.  (5.64);  this 
is  a  rather  crude  approximation  for  the  last  term  in  Eq.  (5.90),  since  T  really 
depends  on  the  distribution  function  itself  as  given  by  Eq.  (5.  66).  Much  of  the 
essential  radiation  belt  electron  physics  is  retained,  however,  except  when  T  cc  is 
the  smallest  time  scale.  Notice  how  well  the  electron  slot  region  that  separates  the 
radiation  zone  into  two  belts  is  reproduced  compared  with  actual  data  from  the 


48.  Lyons,  L.  H.  ,  and  Thorne,  R.  M.  (1973)  Equilibrium  structure  of  radiation 

belt  electrons,  J.  Geophvs.  Res.  78:2142. 
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Figure  5.24.  Predicted  Radiation  Belt  Electron 
Pitch  Angle  Distributions.  48  The  example  shown 
was  calculated  using  a  Gaussian  ELF  whistler 
mode  hiss  turbulence  frequency  distribution 
peaked  at  600  Hz  and  with  a  bandwidth  of  300  Hz 


OGO-5  spacecraft,  both  in  width  and  L-shell  location  for  different  energies.  Per¬ 
fect  agreement  in  all  details  should  not  be  expected,  since  the  theory  treats  aver¬ 
age  conditions  based  on  parametrized  wave  characteristics,  while  the  data  are  for 
a  specific  time  when  the  space  observations  were  made. 

From  such  calculations  we  have  learned  that  for  energetic  electrons,  wave- 
particle  interactions  are  the  cause  of  the  division  into  two  distinct  belts,  since 
electrons  in  the  "slot"  region  are  lost  due  to  pitch  angle  scattering  into  the  atmos- 
sphere  at  a  fast  rate.  Those  few  electrons  that  survive  the  cross-L  transport  to 
arrive  in  the  inner  radiation  zone  became  very  stable  since  the  time  scale  for  wave- 
particle  scattering  is  very  long  there. 
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Figure  5.  25.  Theoretical  and  Observed 
Radiation  Belt  Electron  Radial  Flux  Dis¬ 
tributions  at  the  Geomagnetic  Equator.  The 
theoretical  results  (dashed  curves)  are  cal¬ 
culated  from  radial  diffusion  of  energetic 
electrons  subject  to  pitch  angle  scattering 
due  to  the  (resonant)  ELF  whistler  mode 
turbulence  (hiss)  within  the  plasmasphere. 

The  data  (solid  lines)  are  taken  from  the 
OGO-5  spacecraft.  Notice  how  well  the 
theory  simulates  the  location  and  width  of 
the  observed  radiation  belt  separation  (slot) 
region  over  a  wide  range  of  electron  energies, 
leaving  essentially  only  the  finer  details  of 
the  equilibrium  distribution  for  future 
investigation^*” 


5.5.  1.3  PROTON  MODEL 

Empirically,  radiation  belt  protons  are  not  distributed  into  two  radiation  zones. 
At  constant  proton  energy  there  is  generally  a  single  peak  in  the  proton  flux  radial 
distribution  and  no  specific  wave  mode  has  yet  been  identified  that  interacts  in  a 
dominant  way  with  _  100  keY  protons  below  the  average  plasmapause  location.  On 
the  other  hand,  protons  (like  all  ions)  are  subject  to  the  charge  exchange  process. 
Therefore,  if  we  neglect  the  pitch  angle  scattering  process,  we  can  write  the 
radial  diffusion  equation  for  equatorially  mirroring  -  tt / 2 )  protons 
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(5.  92) 


-  V  ‘  Gp(L) 


where  A  ^  =  <  o^q  [H]  v  >  is  the  charge  exchange  loss  frequency  per  unit  distribu¬ 
tion  function  f.  To  simulate  the  steady  state  radiation  belts  (that  is,  -*  0)  this 
equation  can  likewise  be  solved  when  suitable  boundary  conditions  are  imposed. 

The  boundary  conditions  are  established  by  using  a, measured  outer  zone  proton 
flux  spectrum,  and  by  covering  an  energy  range  sufficiently  large  that  f(L,  ft)  -  o 

at  ii  =  H  where  p  corresponds  to  the  upper  range  of  the  radiation  belt  trap- 
ffidx  max 

ping  energies.  Figure  5.  26  shows  the  computed  radiation  belt  proton  spectra  ob¬ 
tained  with  this  method,  and  data  are  shown  for  comparison  (takenfrom  Explorer  45 
below  L  =  5.3).  Notice  that  the  spectral  maximum  found  at  L  <  5  generally  is 
displaced  towards  higher  energies  with  lower  L-shells.  The  flux  values  computed 
below  a  few  tens  of  keV  are  below  radiation  belt  energies  (convective  processes 
also  operate  in  that  range)  and  should  not  be  considered  reliable.  The  comparison 
with  data  is  quite  good,  however,  giving  support  to  the  usefulness  of  the  radial 
diffusion  theory  at  least  above  100  keV. 


Figure  5.  20.  Radiation  Belt  Proton 
Energy  Spectra  Obtained  From  a 
Theoretical  Calculation  Based  on 
Proton  Radial  Diffusion  Subject 
to  Coulomb  Collisions  and  Charge 
Exchange  Losses.  No  wave- 
particle  interactions  were  con¬ 
sidered.  Boundary  conditions  were 
imposed  at  L-6.6  with  data  from 
the  ATS-6  spacecraft,  and  the 
computed  spectra  at  lower  L-shells 
are  compared  with  available 
equatorial  data  from  Explorer  45. 
'Die  results  of  the  theoretical 
calculations  below  a  few  tens  of 
keV  energies  may  be  inaccurate 
since  convection  processes  may 
dominate  over  diffusion  at  these 
low  energies^” 


ENERGY  (  keV) 


A  complementary  view  of  the  theoretical  radiation  belt  proton  fluxes  is  given 
in  Figure  5.27,  which  depicts  the  computed  radial  distributions  for  different 
proton  energies.  Oualitativelv  similar  to  the  observations,  the  radial  flux  maxima 
are  displaced  towards  lower  L-shells  with  higher  proton  energy.  On  L-shells 
well  beyond  the  flux  radial  peak  location,  the  characteristic  time  scales  of  diffusive 
transport  are  substantially  shorter  than  those  of  the  loss  processes.  This 
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Figure  5.27.  Theoretical  Proton  Radial 
Distributions  Calculated  From  Proton 
Radial  Diffusion  Subject  to  C'ou'omb  Colli¬ 
sion  Energy  Tosses  and  Charge  Exchange 
Loss.  Notice  there  is  no  division  (or  slot) 
within  the  trapping  region.  The  curves 
now  show  that  the  radial  flux  maximum  is 
displaced  towards  lower  L-shells  with  higher 
proton  energy.  For  details  of  the  numerical 
calculations  see  Spjeidvik^8 

situation  is  called  diffusive  equilibrium.  There  is  also  an  important  connec¬ 
tion  between  the  radial  location  of  the  proton  flux  maximum  at  a  given 
energy  and  the  magnitude  of  the  radial  diffusion  coefficient.  As  D^j  increases, 
the  diffusive  equilibrium  radial  range  extends  towards  lower  L-shells,  and  the 
radial  peak  location  is  found  at  lower  L.  Also,  the  absolute  magnitude  of  the  flux 
radial  peak  increases  as  Dj^  increases.  An  enhanced  D^l  causes  a  substantial 
increase  in  the  trapped  fluxes  on  L-shells  below  the  peak  flux  location.  Proton 
models  that  treat  other  than  equatorially  mirroring  pitch  angles  in _a  comprehensive 
manner  have  not  yet  been  developed.  The  reason  is  partially  that  a  is  not 
known  for  protons,  and  partially  that  solving  the  simultaneous  proton  radial  and 
pitch  angle  diffusion  problem  subject  to  the  simultaneous  -variations  is  mathe¬ 
matically  difficult  owing  to  the  great  inhomogeneity  in  the  coefficients. 

5.5.  1.4  HEAVY  ION  MODELS 

Besides  electrons  and  protons,  the  Earth's  radiation  belts  contain  appreciable 
ion  fluxes  of  energetic  helium,  carbon,  oxygen,  and  other  ions.  Theoretical  model¬ 
ing  of  these  ion  populations  can  be  carried  out  by  a  generalization  of  the  proton 
model.  One  must  consider  a  distribution  function  for  each  charge  state  of  a  given 
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radiation  belt  species,  as  well  as  the  coupling  between  charge  states  imposed  by 
the  charge  exchange  chemistry.  For  helium  ions  one  may  write  the  coupled  equa¬ 
tions  for  equatorially  mirroring  ions: 


■  \ofl  +  '2/2  '  Al2f  1 


+  t;  tu  p 


(5.  93) 


_  A  f  +  A  f 
212  121 


+  G  2^  ^  M- 


(5.  94) 


where  the  subscripts  indicate  the  ion  charge  state. 
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These  equations  have  also  been  solved  numerically,  ’  and  some  of  the  find¬ 
ings  are  illustrated  in  Figure  5.  28.  Notice  that  there  is  a  transition  from  charge 
state  1  (He+)  to  the  state  2  (He2+)  at  energies  around  1  MeV,  so  that  the  lower 
charge  state  of  helium  ions  is  most  abundant  at  the  lower  energies  and  the  higher 
charge  state  at  the  higher  energies.  The  spectral  features  computed  for  helium 
ions  below  ~  1  MeV  stem  from  the  energy  dependence  of  the  radial  diffusion  coeffi¬ 
cient  L>l^  and  the  loss  rates. 

For  even  heavier  ions  a  multitude  of  charge  states  are  available  to  the  radiation 
belt  heavy  ions:  for  carbon  ions  there  are  six  positive  charge  states,  and  for  oxy¬ 
gen  ions  there  are  eight.  In  general,  for  a  given  ion  charge  state,  i,  of  an  elemental 
species  with  s  available  states,  the  radial  diffusion  equation  for  equatorially 
mirroring  ions  may  be  written 

(m)  ‘  L2(m){D1.LiI-"2(ft-)  }  *  oau*-'12  Q.) 
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where  A.,  denotes  charge  state  transformation  from  state  i  to  state  j,  so  that 

lJ  49 

A  _  =  0  for  i  =  j  [Spjeldvik  ]. 

49.  Spjeldvik,  W.  N.  (1979)  Expected  charge  states  of  energetic  ions  in  the 
magnetosphere.  Space  Sci.  Kev.  23:499. 
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Figure  5.  28.  Theoretical  Helium  Ion  Radial  Distributions  Calculated 
From  Radial  Diffusion  Theory  With  Coulomb  Collisions  and  Charge 
Exchange.  Solid  curves  show  equatorially  mirroring  fluxes  of  He++ 
while  the  dashed  curves  depict  He4.  Notice  that  He4  is  dominant 
below  ~  1  MeV  while  He++  is  dominant  at  higher  energies^ 


The  radial  diffusion  equation  has  been  solved  for  equatorially  mirroring  radia- 
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tion  belt  oxygen  ions  by  Spjeldvik  and  Fritz,  and  one  of  the  findings  is  that  oxygen 
ions  can  be  more  numerous  than  protons  at  multi  MeV  energies.  The  oxygen  ion 
charge  state  distribution  is  predicted  to  vary  from  a  dominance  of  Of  ions  at 
energies  below  ~  100  keV  to  successively  higher  charge  states  at  the  higher  energies; 
for  example  at  4  MeV  the  fourth  charge  state  should  be  most  abundant.  A  significant 
finding,  illustrated  in  Figure  5.  29,  is  that  the  radiation  belt  charge  state  distribution 
in  the  interior  of  the  radiation  belts  (L  *  5)  is  almost  completely  independent  of  the 
charge  state  distribution  of  the  outer  radiation  zone  source.  In  fact,  this  figure 
shows  that  regardless  of  the  source  being  ionospheric  (source  charge  state  1)  or 
solar  (source  charge  state  6  dominant),  a  steady  state  radiation  belt  distribution  over 
ion  charge  states  that  is  independent  of  the  source  characterists  evolves  in  the 
interior  of  the  trapping  region. 
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Figure  5.  29.  Relative  Charge  State  Distribution  for 
Radiation  Belt  Oxygen  Ions.  Left  panel:  an  ionos¬ 
pheric  source  ol'O*  is  assumed;  Right  panel;  a  solar 
(wind)  source  of  0®+  is  assumed.  One  finds  the 
steady  state  radiation  belt  charge  state  distribution 
in  the  interior  of  the  trapping  region  becomes  largely 
independent  of  the  source  charge  state  characteristics 


The  findings  concerning  charge  states  of  energetic  radiation  belt  ions  may  be 
summarized  as  follows: 

(1)  Higher  charge  states  are  crucially  important  for  the  overall  structure  of 
the  heavy  ion  component  of  the  Earth's  radiation  belts  above  -  100  keY. 

(2)  Charge  state  redistribution  processes  are  of  major  importance  throughout 
much  of  the  radiation  belts  and  at  all  energies. 

(3)  Radiation  belt  charge  state  distributions  become  largely  independent  of 
source  charge  state  characteristics  because  of  the  frequent  "internal"  charge 
exchange. 

(4)  Relative  charge  state  distributions  are  to  a  large  extent  independent  of  the 
diffusive  transport  rate  in  much  of  the  inner  magnetosphere. 

The  last  two  conclusions  are  valid  below  a  charge  state  redistribution  zone  adjoining 
the  particle  injection  boundary  region  in  the  outer  radiation  zone. 

5.5.  1.5  THEORETICAL  RADIATION  BELT  IONIC  COMPOSITION 

Because  of  experimental  difficulties,  only  very  limited  information  is  available 
concerning  the  actual  composition  of  radiation  belt  ion  fluxes.  For  this  reason,  we 


□  resent  a  theoretical  prediction  of  radiation  belt  protons,  helium,  ami  oxygen  ion 
fluxes  at  L  =  3.  35  (in  the  center  of  t lie  radiation  belts).  This  is  shown  in 
Figure  5.30.  It  must  be  kept  in  mind  that  these  predictions  have  vet  to  be  verified 
or  refuted  bv  experimental  observation,  particularly  in  the  100-1000  ke\  range. 

■>  i 

Certain  local  spectral  minima  are  predicted  based  on  Cornwall's-  radial  diffusion 
coefficients,  but  with  other  theoretical  coefficients  the  spectra  mav  be  different. 
The  prediction  is  that  protons  should  be  the  dominant  ion  species  from  0.  1  to 
1  MeV,  helium  ions  should  dominate  at  1  to  several  Mc\  per  ion,  and  at  higher 
energies  oxvgen  (and  possibly  also  carbon)  should  dominate.  1  hese  comparisons 
are  done  at  equal  total  ion  energy.  If  comparison  is  made  at  equal  energy  per 
nucleon,  then  protons  would  be  the  dominant  constituent  at  almost  all  radiation  belt 
energies  beyond  a  few  tens  of  keV. 


Ion  Energy  (keV) 


Figure  5.  30.  Theoretical  Model  of  the  Radiation  Belt  Ion 
Fluxes;  Comparison  of  Theoretically  Predicted  Knergv 
Spectra  of  Proton,  Helium  and  Oxvgen  Fluxes  at  I.  3.35 
for  Different  Diffusion  Coefficients^*' 

5.5.2  Geomagnetic  Storms 

Most  radiation  belt  modeling  carried  out  so  far  deals  with  equilibrium  condi¬ 
tions  appropriate  for  quiet  time  conditions.  There  is  at  the  present  time  no  fully 
comprehensive  radiation  belt  model  for  the  effects  of  geomagnetic  storms  and  other 


disturbances.  The  reason  is  partially  the  incomplete  knowledge  of  the  radiation 
belt  particle  source  mechanisms,  and  in  part  that  disturbed  time  modeling  is  diffi¬ 
cult  because  of  the  changes  in  magnetic  field  and  time  varying  electric  fields  that 
must  be  also  incorporated.  Research  is  being  conducted  in  this  area,  but  no  defi¬ 
nite  models  are  available  as  of  this  writing. 

Specific  types  of  disturbances  and  associated  time  variability  have,  however, 
been  analyzed.  Among  these  are  studies  of  the  post-geomagnetic  storm  decay  of 

energetic  particle  fluxes  following  the  (yet  not  fully  explained)  injection  events. 
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Lyons  and  Thorne  have  demonstrated  that  the  radiation  belt  electron  particle 
scattering  lifetime  can  be  used  to  explain  quantitatively  the  restoration  of  the  radia¬ 
tion  belt  electron  slot  region  after  it  becomes  filled  in  during  the  storm  injection 
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process.  Spjeldvik  and  Thorne  subsequently  demonstrated  that  the  precipitating 
electron  flux  into  the  middle  latitude  ionosphere  following  such  storms  causes  en¬ 
hanced  D- region  ionization  of  sufficient  magnitude  to  explain  V LF  radio  wave  phase 

5  1 

anomalies  observed  at  such  times.  Spjeldvik  and  Lyons  have  suggested  a  simpli¬ 
fied  prediction  model  for  these  effects. 

There  is  reason  to  think  that  magnetic  storms  differ  considerably  from  one 
another,  not  only  by  the  magnitude  of  the  ring  current  storm  index  Dgt  but  also  in 
the  injection  characteristics  of  energetic  ion  fluxes.  Once  injected,  radiation  belt 
particles  become  subject  to  the  normal  processes  in  the  trapping  region  discussed 
earlier.  For  specific  magnetic  storms,  the  post-storm  decay  of  protons  and  heavy 

ions  has  also  been  studied,  and  fair  agreement  between  predictions  and  observations 
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has  been  reported  in  limited  energy  ranges  where  the  data  were  available.'  *' 


5.5.2.  1  TIME  VARIABILITY  OF  THE  PARAMETERS 

Within  the  framework  of  radiation  belt  diffusion,  the  injection  of  particles  is 

described  by  a  (time  variable)  source  term  S(L,  p,  cv^rt).  The  boundary  conditions 

may  also  be  time  var  able,  reflecting  dynamic  conditions  on  the  outer  edge  of  the 

stable  trapping  region  f(L  ,  ,  p,  a  ; t) ,  and  the  transport  coefficients,  U  and 

_  max  o  1,1. 

a  ,  will  also  be  time  variable  reflecting  the  geomagnetic  activity  level.  Based 
°  °  54 

on  limited  data,  Lanzerotti  et  al  estimated  that  the  radial  diffusion  coefficient 
might  vary  with  the  K  geomagnetic  index  as 
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(Due  to  the  large  number-  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  pagp  127.  ) 
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It  is  not  yet  known  liow  1 J ^  j  and  D^  a  might  vary  with  geomagnetic  conditions. 

u  o 

It  is  also  known  that  the  exospheric  neutral  density  varies  with  activity  (because  of 
the  heating  of  the  upper  atmosphere',  and  that  the  plasmuspheric  densities  can  vary 
strongly.  A  first  approach  may  be  to  solve  the  appropriate  transport  equations 
using  perturbation  theory  for  the  ditferent  variables.  Some  geomagnetic  conditions 
may,  however,  be  too  drastically  altered  to  be  treated  as  perturbations,  so  caution 
is  in  order.  At  the  present  time  there  are  many  unknowns,  and  specific  models 
have  yet  to  be  developed. 

5.5.2.  2  MAGNETIC  TOPOLOGY  VARIATIONS 

During  the  early  phase  of  magnetic  storms  the  earth's  magnetic  field  becomes 
compressed  on  the  dayside.  It  has  been  observed  that  the  subsolar  magnetopause 
can  be  pushed  inward  from  an  average  location  of  ~  10  R  .  to  ~5  R  during  large 
storms.  This  implies  a  large  H-field  change.  However,  most  magnetic  storms  do 
not  cause  such  a  large  perturbation  of  the  magnetic  field.  Increase  in  the  lower 
energy  (~  10  -  500  keV)  radiation  belt  ion  fluxes  produces  an  enhanced  ring  current 
(see  Section  5.  7)  which  depresses  the  earth's  magnetic  field  earthward  of  the  en¬ 
hanced  particle  flux  region,  and  causes  adiabatic  deceleration  of  the  trapped  radiation 
belt  particle  fluxes.  On  the  other  hand,  geoelectric  fields  penetrating  into  the  trapping 
region  during  disturbed  conditions1^’  ^  can  cause  cross-L  non-diffusive  transport 
and  thus  adiabatic  acceleration  of  the  particle  population.  The  relative  influence 
of  these  processes  depends  on  particle  energy.  These  effects  must  also  be  incor¬ 
porated  into  storm-time  radiation  belt  prediction  models. 

5.6  EMPIRICAL  RADIATION  BELT  MODELS 

Since  the  discovery  of  the  earth's  radiation  belts,  the  population  of  trapped  par¬ 
ticles  has  been  measured  with  ever-improving  instrumentation,  ’’’he  early  Geiger 
counters  flown  on  the  first  few  spacecraft  had  little  or  no  particle  identification 
capabilities.  Subsequent  instruments  used  foil  techniques,  solid  state  detectors, 
magnetic  and/or  electric  deflection  techniques  and  electronic  signal  discrimination. 
This  led  to  a  clear  separation  of  electrons  and  ions  (which  were  then  assumed  to  be 
solely  protons).  During  the  later  years  of  space  exploration  the  presence  of  a 
multitude  of  different  ion  species  was  established.  The  purpose  of  this  section  is 
to  provide  a  brief  overview  of  existing  radiation  belt  data,  as  a  source  guide  for 
those  who  require  crude  numerical  estimates  of  the  radiation  environment.  It  is 
emphasized  that  the  older  empirical  proton  model  in  reality  represents  contribu¬ 
tions  from  many  ion  species,  and  that  sometimes  the  heavy  ion  contribution  may  be 
dominant. 


5.6.1  Data  Acquisition  and  Processing 


Empirical  radiation  belt  models  are  compiled  by  NASA / National  Space  Science 
Data  Center,  Goddard  Space  1' light  Center  in  Maryland  |  Vette  et  al,  °°  and  references 
therein].  These  models  represent  a  systematic  effort  to  compile  many  years  of 
data  containing  a  large  number  of  disparate  satellite  observations  into  a  few  key 
models.  These  observations  were  separated  in  space  and  time,  and  made  with 
highly  varying  instrumentation  so  that  subjective  judgments  were  necessary  regard¬ 
ing  data  quality.  The  complexity  of  the  task  is  appreciated  if  one  considers  the 
volume  of  space  to  be  covered  and  the  time  variation  in  the  particle  fluxes.  Most 
of  the  data  used  in  the  NASA  models  were  obtained  in  the  1960’s  and  early  1970  s. 

For  example,  the  present  AP-8  model  for  protons  is  a  combination  of  94  different 
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instrument  energy  detector  channels  from  instruments  flown  on  24  satellites. 

The  combined  data  sets  were  smoothed  to  obtain  omni-directional  flux  distributions 
in  the  13,  L  parameter  space.  The  fluxes  were  integrated  over  all  pitch  angle"  and, 
therefore,  directional  information  was  not  utilized.  Other  limitations  in  the  data 
sets  are  variations  in  instrument  response  and  the  lack  of  local  time  dependence 
considerations.  Also  note  that  the  flux  models  are  compiled  from  a  data  base 
obtained  over  a  brief  time  period  in  comparison  with  other  geophysical  and  solar 
time  scales.  Long  term  dynamical  changes  in  the  radiation  belts  are  not  well  under¬ 
stood  so  that  extrapolation  to  different  epochs  must  be  done  cautiously.  For  ex¬ 
ample,  much  higher  fluxes  than  the  models  give  have  recently  been  observed.  One 

example  of  very  long  term  effects  is  the  adiabatic  energization  of  inner  belt  protons 
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by  the  secular  variation  of  the  earth's  magnetic  field. 

5.  6.  1.1  PHOTONS  (IONS) 

It  has  become  customary  to  refer  to  energetic  protons  located  below  L  -  2  as 

inner  belt  or  inner  zone  protons.  This  is  a  misnomer,  however,  since  the  proton 
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fluxes  normally  have  only  one  radial  flux  peak.  Sawyer  and  Vette  have  completed 
an  extensive  model  of  the  trapped  "proton"  environment  out  to  L  =  6.6  for  energies 
above  100  keV.  The  measured  ions  were  labeled  "protons",  although  no  actual  ion 
identification  was  made.  The  model  is  composed  of  two  parts,  AP8MIN  and 
AP8MAX,  which  correspond  to  observations  made  during  the  1960-1970  sunspot 


55.  Vette,  J.  I.  ,  Teague,  M.J.,  Sawyer,  D.M.,  and  Chan,  K.W.  (1979)  Modelling 

the  earth's  radiation  belts,  Solar-Terrestrial  Prediction  Proceedings,  Vol.  2, 
(Ed.  R.  F.  Donnelly),  U.  S.  department  of  Commerce,  N’OAA,  Boulder,  Colorado. 

56.  Sawyer,  D.  M. ,  and  Vette,  J.  I.  (1976)  AP-8  Trapped  Proton  Environment 

for  Solar  Maximum  and  Solar  Minimum,  NSSDC /VVDC-A-R&.S76-06, 


minimum  and  maximum  conditions,  repectively.  The  difference  between  these  two 
models  is  thought  to  result  in  part  from  differences  in  upper  atmosphere  heating 
during  solar  active  periods  such  that  the  trapped  particle  collision  rate  (due  to  the 
exospheric  expansion)  was  increased.  AP8MAX  differs  from  AP8MIN  for  altitudes 
less  than  about  1000  km  and  for  L  values  less  than  2.  9,  however,  steep  spatial 
gradients  in  the  ion  fluxes  at  lower  altitudes  are  difficult  to  determine  accurately. 

An  equatorial  profile  of  the  ion  (proton)  fluxes  as  given  by  AP8MIN  is  shown  in 
Figure  5.  31.  Note  the  presence  of  >  400  MeV  protons  for  L  <  2.  This  is  consistent 
with  higher  energy  protons  being  produced  by  the  decay  of  neutrons  produced  in  the 
atmosphere  by  cosmic  rays  (the  GRAND  source).  The  lower  energy  fluxes  can 
arise  from  inward  radial  diffusion  as  discussed  in  the  theoretical  modeling  section. 

AP  8  MIN  equatorial  omnidirectional  radial  profile 
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Figure  5.31.  Radial  Distribution  of  AP8IW1X  <  imnidirectional  I  luxes  of  Protons 
in  the  Equatorial  Plane  With  Energies  Above  Threshold  Values  Between  0.  1  and 
400.  0  MeV56 


The  AP8  proton  models  include  no  data  after  1970  nd  ve little  data 
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100  MeV  energies.'  Although  the  proton  belt  is  considered  far  more  st.au-  ‘ban 
the  electron  belts,  significant  enhancement  and  depletions  have  been  observed.  ■> 
secondary  equatorial  proton  peak  (40-  110  MeV)  of  unknown  origin  was  observed  b> 

"  CQ 

Mcllwain.  This  peak  moved  from  I.  -  2.25  to  L  =■  2.  1  earth  radii  between  .lanuur\ 

1963  and  January  1965.  During  the  large  Ma,\  1967  magnetic  storm,  2.  2  -  2  .Me'. 

proton  flaxes  were  observed  to  increase  more  than  a  factor  of  10  at  I  -2.2 
e  n 

[Bostrom  et  al  ].  Lower  energy  (0.  26  -  0.  65  MeV)  protons  were  observed  to  be 
preferentially  enhanced  during  the  same  storm  at  L  =  3.  0  |Kothwell  and  Katz  ). 

Beyond  L  =  5  earth  radii,  order-of-magnitude  fluctuations  occur  on  time  scales  as 

54  , 

short  as  10  min.  Very  large  increases  in  lUeV  heavy  ion  fluxes  at  I.  =  2  -  a 
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were  observed  during  the  August  1972  magnetic  storm  event.  *'  l  or  this 

reason  engineering  applications  should  consider  large  deviations  from  the  mean 
models. 

Off-equatoriallv  mirroring  ions  intercept  the  upper  atmosphere  at  higher  equa¬ 
torial  pitch  angles  on  lower’  I. -shells,  l  or  example.  Figure  5.  32  shows  the  B  B 
distribution  for  three  different  L-values  just  above  the  atmosphere;  B/B  being  the 
ratio  of  local  to  equatorial  magnetic  field  induction  on  the  same  field  line.  At 
L  =  1.  17  the  equatorial  pitch  angle  distributions  have  very  steep  loss  cone  gradients 
such  that  the  particle  flux  vanishes  atQ^^f,  73°.  At  1.  =  1.50  this  "cutoff  is 
a  ~  37°  and  for  L  =  2.00,  Of  -  21°.  These  cutoffs  come  about  because  the  atrnos- 
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pheric  bounce  loss  cone  is  wider  at  the  lower  l. -shells. 

The  South  Atlantic  Magnetic  Anomaly  has  a  controlling  effect  on  the  inner  zone 

particle  fluxes  in  the  vicinity  of  the  loss  cone.  This  anomaly  arises  from  the  earth  s 

magnetic  field  being  less  intense  at  a  latitude/longitude  region  located  near  the 

coast  of  Brazil.  Trapped  particles  at  these  low  1. -shells  will  encounter  their  lowest 

mirroring  altitude  (H  .  )  and  thus  the  densest  atmosphere  in  this  longitude  region. 

-  2  - 1  - 1 

Figure  5.  33  shows  omnidirectional  proton  flux  contours  in  protons  cm  sec  MeV 
at  750  km  altitude.  These  contours  were  derived  from  5  to  7  MeV  (dashed  lines)  and 
28  to  45  MeV  (solid  lines)  data  channels  on  the  AFG1.  particle  identifier  instrument 
flown  on  AF  Satellite  72-1  in  1972.  The  data  show  that  the  location  of  the  proton  peak 
flux  in  the  low  altitude  section  of  the  trapping  region  adjoining  the  South  Atlantic  is 
dependent  on  energy.  The  5  to  7  MeV  peak  is  located  around  1.  -  2  while  the  28  to 
45  MeV  peak  is  located  around  I.  =  1.  3. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  127.) 
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The  azimuthal  drift  around  the  earth  of  radiation  belt  particles  through  the 
South  Atlantic  anomaly  region  produces  a  "windshield  wiper"  effect.  Electrons 
present  in  the  enhanced  loss  cone  created  by  the  locally  low  B- field  region  are  lost. 
These  electrons  are  subsequent'y  replenished  by  pitch  angle  diffusion  during  the 
drift  around  the  earth  outside  the  anomaly. 

•Solar  cycle  variations  that  modify  the  atmosphere  and  ion  composition  and  high 
altitude  nuclear  detonations  can  significantly  affect  the  trapped  particle  populations. 
A  15  year  time-study  of  inner  belt  55  MeV  protons  concluded  that  the  observed  flux 
variations  were  consistent  with  expected  atmospheric  loss  processes. 

5.6.  1.2  HEAVY  IONS 

The  abundance  of  trapped  heavy  ions  suggest  clues  to  the  origin  of  the  radiation 
belt  particles.  At  high  energies  the  heavy  ions  also  constitute  a  hazard  to  opera¬ 
tional  space  systems.  Over  the  last  decade  data  have  been  acquired  to  sketch  the 
spectral  and  angular  dependence  of  the  helium  and  carbon-nitrogen-oxygen  (CNO) 
fluxes. 

Helium  Ions:  Helium  ions  are  sometimes  referred  to  as  alpha  particles, 

-  2+ 

although  in  a  strict  sense  an  alpha  particle  is  only  the  totally  ionized  state  (He  ) 

of  a  helium  ion.  figure  5.  34  shows  average  equatorial  helium  ion  spectra  over  a 
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range  of  L-shells  during  quiet-time  conditions.  Notice  that  these  spectra  are 
very  steep  at  L  >  4,  become  much  harder  with  lower  L-shells,  and  are  almost  flat 
at  L  -  3.  This  is  expected  since  the  lower  energy  helium  ions  are  preferentially 
lost  through  the  ion  charge  exchange  mechanism  and  Coulomb  energy  degradation  as 
the  particles  diffuse  radially  inward.  The  observed  ratio  between  helium  ion  and 
proton  fluxes  simultaneously  observed  is  given  in  Figure  5.35.  Care  must  be 
exercised  whether  the  ratio  is  defined  at  the  same  total  ion  energy  (as  in  this 
figure)  or  at  the  same  energy  per  nucleon.  At  higher  total  ion  energies  the  heavy 
ion  flux  may  be  dominant  over  the  proton  flux  at  some  L-shells  (Figure  5.  35) 
while  at  the  same  energy  per  nucleon,  proton  fluxes  usually  dominate. 

63.  Parsignault,  D.  R. ,  Holeman,  E.  ,  and  Filz,  R.  C.  ( 1981)  Solar  cycle 

modulation  of  the  55-lUeV  proton  fluxes  at  low  altitudes,  J.  Geophys.  Res. 
86:7749. 
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64.  Fritz,  T.  A.,  and  Spjeldvik,  W.  N.  ( 1979)  Simultaneous  quiet  time  observations 

of  energetic  radiation  belt  protons  and  helium  ions:  the  equatorial  a/p  ratio 
near  1  MeV,  J.  Geophys.  Res.  84:2608. 
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Figure  5.  34.  Energy  Spectra  of  Equatorial  Radiation  Belt  Helium 
Ions  Deduced  From  Mass  Ion  Observations  on  Explorer  45  During 
the  Geomagnetically  Quiet  Period  June  1-15,  1972.  The  data  are 
given  at  L  =  2.  25,  2.  5,  2.  75,  3,  3.  25,  3.  5,  4,  4.  5,  and  5  [F  ritz 

and  Spjeldvik64] 


The  equatorial  helium  ion  pitch  angle  distribution  is  generally  more  anisotropic 

than  the  comparable  proton  pitch  angle  distribution.  For  L-values  between  2.  5  and 

4.  5  the  proton  fluxes  most  often  vary  as  sinn  a  where  n  =  4  ±  2  and  a  is  the 

o  o 

equatorial  pitch  angle,  while  for  helium  ions  n  —  10  t  4.  Beyond  L.  -  5  the  helium 
fluxes  are  quite  dynamic  and  characteristic  quiet-time  values  are  difficult  to  define. 
During  magnetic  storms  the  fluxes  of  energetic  helium  ions  can  increase  by  orders 
of  magnitude  in  the  heart  of  the  radiation  belts  (I.  -  2.5),  and  the  relative  abundance 
of  the  different  ion  species  can  vary.  At  higher  I. -shells  substorm  effects  can  be 
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Figure  5.35.  Quiet  Time  He/p  Ion  Flux  Ratios  in  the  Equatorial 
Radiation  Belts  of  the  Earth  Deduced  From  Explorer  45  Observa 
tions  During  June  1-15,  1972.  The  ratios  are  calculated  from 
flux  observations  at  equal  energy  per  ion.  The  shaded  areas 

mental  uncertainty  in  the  data  [from  Fritz  and 


significant.  For  example,  during  the  June  18,  1974  substorm,  helium  ions  were 
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more  numerous  than  protons  at  geosynchronous  altitude  for  E  2  800  keV/ion. 

6  6 

Blake  and  Fennell  also  have  noted  that  a  strictly  sinusoidal  pitch-angle  distribu¬ 
tion  did  not  completely  describe  the  98-240  keV/nucleon  ion  fluxes  at  the  geo¬ 
synchronous  altitude  (L  ~  6.  6). 

Carbon-Nitrogen-Oxygen  (CNQ):  Trapped  oxygen  ions  can  at  times  be  much 

more  numerous  than  protons  particularly  at  L  2  5  when  compared  at  equal  total 

ion  energy.  Figure  5.  36  shows  the  radial  flux  profiles  of  carbon  and  oxygen  ions 

6  7 

as  measured  by  the  ISEE-1  spacecraft.  At  equal  total  ion  energies  in  the  MeV 

range  the  carbon  to  oxygen  flux  ratio  is  of  the  order  of  0.  5,  and  at  equal  energy  per 

nucleon  the  carbon  flux  can  dominate.  This  indicates  an  extraterrestrial  source 

for  the  very  energetic  trapped  heavy  ions  since  the  ionospheric  C/O -ratio 
*  5  68 

is  <  10  [Blake  J.  The  CNO  flux  pitch  angle  anisotropy  is  even  more  pronounced 
than  that  of  helium,  having  a  value  of  the  anisotropy  n- index  typically  between  12 
and  16.  At  the  higher  total  energies  (>  800  keV)  and  at  geosynchronous  altitudes 
the  differential  CNO  flux  has  been  measured  to  be  higher  than  that  of  protons  and 
helium. 

Very  Heavy  Ions:  Ions  heavier  than  oxygen  are  also  present  in  the  earth's 
radiation  belts  in  small  quantities.  It  has  been  demonstrated  that  substantial  injec¬ 
tion  of  ions  with  nuclear  charge  ZN  >  9  can  take  place  during  some  magnetic  storms, 

and  many  orders  of  magnitude  flux  intensity  enhancements  at  0  Z  10  MeV  have  been 
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observed  lasting  for  many  months.  During  such  disturbed  conditions  the  trapped 
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fluxes  of  other  ions  (He,  O,  and  so  on)  can  also  be  greatly  enhanced. 

5.6.  1.3  TRAPPED  ELECTRONS 

Empirical  flux  models  have  been  developed  that  describe  the  inner  and  outer 
electron  radiation  belts.  The  National  Space  Science  Data  Center  at  NASA/Goddard 
Space  Flight  Center  has  developed  two  inner  belt  models:  AE-6  for  sunspot  maxi¬ 
mum,^  and  AE-5  epoch  1975  for  sunspot  minimum  conditions.  There  is  also  an 
outer  belt  electron  model  called  AE-4.  An  updated  outer  belt  model,  AE-7,  is  now 
under  development.  In  this  section  a  brief  description  of  the  models  is  given  and 
typical  electron  flux  versus  L-shell  profiles,  energy  spectra,  and  pitch  angle 
distributions  are  quoted. 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  127.) 


Figure  5.  36.  Equatorially  Trapped  Carbon  and  Oxygen  Ion  Fluxes 
Measured  From  October  197?  -  January  1978  During  Quiet  Times. 
Count  rates  may  be  ^approximately  converted  to  omni-directional 
flux  units:  ions/(crro  sec  keV)  bv  multiplying  by  50/AE  where  A  E 
is  the  appropriate  energy  passband  in  ke\  60 


Inner  Belt  electrons:  The  empirical  data  AE-5  model  is  based  on  flux  data  from 
five  satellites,  OGO-1,  OGO-3,  1963-38C,  OV3-3  and  Explorer  26  [Teague  and  Vette"  *]. 
This  data  base  covered  the  period  from  December  1964  to  December  1967,  repre¬ 
senting  a  transition  from  solar  (sunspot)  minimum  towards  maximum  conditions. 

During  this  period  the  time-averaged  Zurich  sunspot  number  R  ranged  in  value  from 
about  10  ( December  1964)  to  a  round  100  (December  1  967).  In  constructing  the  AE-5  model 
it  was  assumed  that  the  total  inner  belt  electron  flux  is  composed  of  four  components: 

(1)  quiet  day  flux  at  solar  minimum,  (2)  quiet  day  flux  at  other  times  during  the 
solar  cycle,  (3)  storm-time  flux,  and  (4)  residual  flux  (contamination)  from  the  1962 
high  altitude  .Starfish  nuclear  explosion.  By  late  1967  the  Starfish-generated  energetic 
electron  flux  at  C  <  1  MeV  had  decayed  to  insignificant  levels.  It  should,  however, 
be  noted  that  for  energies  C  >  690  keV  the  available  data  for  the  inner  belt  AE-5 
model  were  quite  limited.  For  example,  within  the  observing  range  of  the  instru¬ 
ments  the  monthly  averaged  OGO-1  and  OGO-3  data  showed  non-zero  counts  in  this 
energy  range  only  when  the  Starfish  fluxes  were  still  present  or  during  storm  times. 


7  1.  Teague,  M.  J.  ,  and  Vette,  ,J.  I.  ( 1972)  The  Inner  Zone  Electron  Model  AE-5, 
NSSDC / WDC-A- H&S  72-10. 
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The  quiet  day  sola r  evcle  variation  was  defined  by  taking  the  ratio  of  the 
omnidirectional  flux  measured  from  solar  minimum  to  a  standard  reference  epoch 
(chosen  as  October  11*671.  Insufficient  data  necessitated  the  assumption  that  the 
B-field,  and  therefore  the  particle  pitch  angle  distribution,  is  independent  of  the 
solar  cycle  as  seen  in  energetic  trapped  electrons.  Also,  the  presence  of  Starfish 
residual  energetic  electrons  restricted  the  analysis  of  the  uncontaminated  radiation 
belts  to  later  times  when  these  man-made  radiation  belt  electron  fluxes  had  decayed 
away.  It  should  be  emphasized  that  the  solar  cycle  variation  has  been  determined 
for  only  one  specific  cycle  and  that  the  results  may  not  apply  to  other  solar  cycles. 

The  effect  of  magnetic  storms  on  the  time  averaged  radiation  belt  electron  flux 
depends  on  (1)  the  frequency  at  which  magnetic  storms  occur,  (2)  the  magnitude  of 
the  stormtime  electron  flux  enhancement  over  quiet  time  values,  and  (3)  the  dura¬ 
tion  and  characteristics  of  the  magnetic  storms.  The  inner-belt  electron  flux 
increases  are  infrequent  but  substantial  and  long -lasting,  so  that  a  representative  model 
storm  effect  was  not  extractable  from  this  NASA  data  base.  However,  the  average 
long  term  impact  of  magnetic  storms  was  estimated  by  forming  the  ratio  R  which 
is  the  average  electron  flux  (June  1966  -  December  1967)  divided  by  the  quiet-time 
electron  flux  (October  1967).  The  results  are  shown  in  Figure  5.  37  and  reflect 
a  pronounced  flux  peaking  at  1  MeV.  Of  course,  the  largest  relative  storm  time 
energetic  electron  flux  enhancements  are  found  in  the  electron  "slot-region" 
located  at  L  ~  2-4,  depending  on  energy. 

An  inner  belt  electron  model  for  solar  minimum  was  derived  from  the  model 
discussed  above.  This  is  called  the  AE-5  Epoch  1975  Projected  Model  and  was 
constructed  by  numerically  removing  the  estimated  temporal  variations.  For 
example,  the  Starfish  residual  energetic  electron  flux  component 

(L  <  1.6  K^,  E  >  700  keY'  was  "removed",  from  the  data  base,  using  the  theoretical 
residual  Starfish  electron  model  of  Teague  and  Stassinopolous. 

Similarly,  a  solar  maximum  model,  AE-6,  was  constructed  using  the  AE-5 
model  at  solar  maximum  values  (epoch  1967)  and  with  the  estimated  Starfish  resi¬ 
dual  (background)  energetic  electron  fluxes  also  subtracted  out.  This  model  is 
called  AE-6  Epoch  1980. 

The  electron  models  were  later  verified  by  comparing  them  with  additional  data 
sets  from  the  spacecraft  OV3-3,  OV1-13,  OV1-19,  OGO-5,  OGO-1,  OGO-3, 
1963-038C,  and  OSO-4  (Teague  et  al^].  These  additional  data  were  also  obtained 


72.  Teague,  M.  J.  ,  and  Stassinopoulos,  E.G.  (1972)  A  Model  of  the  Starfish  Flux 
in  the  Inner  Radiation  Zone,  NASA  Goddard  Space  Flight  Center,  Maryland, 

ncr- 5  5T1T - 


73.  Teague,  M.  J.  ,  Schofield,  N.  J. ,  Chan,  K.W.,  and  Vette,  J.  I.  (197  9)  A  Study 
of  Inner  Zone  Electron  Data  and  Their  Comparison  With  Trapped  Radiation 


before  1970  and  substantiate  the  model  empirical  mean  values  for  that  particular 
epoch.  Figure  5.38  shows  a  comparison  between  these  data  sets  and  the  AE-6  and 
AE-5  1967  models  and  represents  quiet  time  energy  spectra  of  equatorially  mirror¬ 
ing  electrons  as  measured  at  L  =  1.5.  The  OGO-l  day  300  1964  results  are 
evidently  seriously  ontaminated  by  the  Starfish  detonation  residual  energetic  elec¬ 
trons.  Otherwise,  these  flux  data  show  fair  agreement  between  the  AE-6  electron 
model  and  the  different  data  sets.  Figure  5.  39  shows  comparison  of  the  empirical 
model  and  measured  equatorial  pitch  angle  distribution  at  L  =  1.4.  The  AE-5  1967 
(solar  maximum)  model  and  the  AE-5  1975  (solar  minimum)  model  bracket  the  data 
within  a  factor  ±  2  to  3,  except  at  very  low  equatorial  pitch  angles. 
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Figure  5.  37.  Estimate 
of  the  Relative  Average 
Effect  of  Substorms  on 
Inner*  Belt  Electron 
Fluxes.  ^9  It  is  cautioned 
that  different  disturb¬ 
ances  can  have  different 
effects 
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The  quiet  day  inner  belt  electron  fluxes  for1  K  <  690  keV  at  1.  3  <  L  <  2.4  can 
be  represented  by  an  analytic  formula  based  on  the  empirical  data.  The  ener¬ 
getic  electron  flux  is  parametrized  as  follows: 

j(a  ,  L.  G  >  =  A(a  ,  L)  G  exp  [-  C  /£  (a  .  L)]  (5.  98) 

J  o  o  o  o 


where  L.  is  the  Mcllwain  L-shell  parameter,  aQ  is  the  equatorial  pitch  angle  and  the 
parameters,  A  (aQ,  L),  CQ(a0»  L),  are  related  to  the  equatorial  values  (aQ  =  90°) 
bv  the  following  empirical  expressions: 

sin  (a  -  a  olc 

Ala  ,  L)  =  A  (90°,  L)  - - -  <S  >  a  S  a 

o  .  m  , ,  ^  x  o  c 

sin  «S  -  ooLC» 

(5.  99) 

=  A  (90°,  L.)  90°  S  a  SO 
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CQ(ao,  L)  =  eo  (9°°.  1A  Sin11  ao/sinn  0 


0  >  a  ~  a 
v  o  c 


=  E  (90°,  l.) 
o 


90°  Sa  So. 
o 


(5.  100) 


Values  for  0,  a  ,  n,  m,  G  (a  =  90°,  l.landAla  =  90°,  L)fora  given  1, -shell 
c  o  o  o 

are  given  in  Table  5.  2. 


In  Kqs.  (5.  99)  and  (5.  100)  the  equatorial  electron  pitch  angle  distributions  are 
described  as  being  fairly  flat  out  to  an  angle  0,  where  they  rapidly  drop  as  sin11  a. 
The  parameter  ^  is  the  minimum  allowed  equatorial  pitch  angle  (the  loss  cone 
angle)  and  corresponds  to  a  100  km  mirroring  altitude.  These  formulas  should  be 
used  with  caution,  however,  since  they  represent  extrapolations  based  on  an  imper¬ 
fect  and  incomplete  data  set.  The  dipole  approximation  gives  aQj  in  terms  of 

B  from  Table  5.  2. 

c  _ 


a 


oLC 


arc  sin 


(5.  101) 


where  B^  is  the  magnetic  induction  at  the  100  km  altitude  level  on  a  given  field  line. 
Figures  5.40  and  5.41  show  the  equatorial  integral  electron  flux  profiles  as  given 
by  AE-5  (1975  projected!  and  Atl-6. 


Figure  5.40.  Empirical  Radiation 
Belt  Electron  Observations;  Equatorial 
Flux  vs  L-Shell  for  the  AE5-MIN 
Model  as  Taken  From  the  National 
Space  Science  Data  Center  (NSSDC) 
Computer  Models.  The  flux  is  in  units 
of  electrons/(cm^  sec),  and  the 
energy  range  is  0.  04  -  4  MeV 


Figure  5.41.  Empirical  Radiation 
Belt  Electron  Observations; 

Equatorial  Flux  vs  I.  Profiles  for 
the  AE-6  Model  as  Taken  From  the 
(NSSUC)  Computer  Models.  The 
flux  is  in  units  of  electrons /(cm^  sec), 
and  the  energy  range  is  0.  04  -  4  MeV 


Outer  Belt  Electrons:  The  outer  belt  trapped  electron  fluxes  are  located  be¬ 
tween  ~3.  5  and  ~  11  R^,.  Rapid  changes  in  the  magnetic  field  and  background  plasma 
can  modify  the  electron  flux  levels  within  minutes.  Because  of  the  apparent  coupling 
between  magnetospheric  substorm  processes  and  the  outer  belt  trapped  electron 
fluxes,  time-averaged  models  have  been  developed.  There  are  discernible  changes 
in  the  average  flux  over  the  time  period  studied  so  that  the  model  is  given  for  two 
epochs.  The  NASA  models  are  called  AE-4  epoch  1964  (sunspot  minimum)  and 
AE-4  epoch  1967  (near  sunspot  maximum).  The  data  base  was  acquired  between 
1959  and  1968  from  23  instruments  on  1 1  satellites. 

Because  of  the  lack  of  azimuthal  symmetry  of  the  geomagnetic  field  in  the  outer 
radiation  zone,  studies  of  the  radiation  belt  electron  structure  beyond  L  ~  5  requires 
that  the  conventional  B-L  coordinate  system  (calculated  from  the  earth's  internal 
magnetic  field)  be  augmented  with  the  additional  coordinate  local  time,  LT.  The 
empirical  LT  dependence  of  the  outer  radiation  belt  electron  fluxes  has  been  deter¬ 
mined  to  be 

log  .1  -  C(E,  L)  cos  [7;  — ~ ■ 1  ]  (5.  102) 
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With  LT  in  hours  and  C(E,  L  ~  5)  =  0  [Vette  et  al  '  |.  But  this  is  true  onlv  in  a 
time  average  sense,  and  it  was  found  that  the  electron  flux  intensitv  levels  ran  varv 
by  orders  of  magnitude  at  any  given  local  time. 

Given  the  equatorial  flux  (o  =  nl  2  or  B  =  B  ),  the  off-equatorial  outer  belt 

o  o  y  * 

integral  electron  flux  4  tt/2)  can  be  estimated  bv: 

J  (>  E,  B,  L]  =  J  [>  E,  B  =  Bo,  L)  G[B,  LI  (5.  103) 


where 


G[B,  L|  =  (B/Bq) 


B  -  B  \  r>,+  1/2 

B  -  B  ) 

c  o  / 


;  B  <  B 


(5.  104) 


G[B,  LI  =  0 


;  B  S  B 

c 


The  parameters  m,  B^  and  B  are  all  empirical  functions  of  1.  and  are  given 
in  1  able  5.3.  The  parameter  B  is  the  magnetic  field  value  at  the  desired  location 
off  the  geomagnetic  equator,  and  B  _  is  the  value  of  the  magnetic  field  at  100  km 
altitude  on  the  same  magnetic  field  line. 


74.  Singlev,  G.  \V.  ,  and  Vette,  .1.1.  (1972)  The  AE-4  Model  of  the  Outer  Radiation 
Zone  Electron  Environment,  \SSI)C  72-06,  NASA. 


Table  5.3.  Flux  Model  Parameters  for 
Outer  Belt  Electrons.  (See  text  for  details)' 


AE-4  B/B  Model  Parameters  (Outer  Belt 
Electrons)0 


L 

(R.J 

E 

m 

B0 

(G) 

Be 

(G) 

3.  00 

1.  12 

0.  01154 

0.  580 

3.  10 

0.  87 

0. 01046 

0.  582 

3.  20 

0.  7  1 

0. 009511 

0.  585 

3.  40 

0.  66 

0. 007929 

0.  588 

3.  60 

0.  63 

0. 006680 

0.  593 

4.  00 

0.  60 

0. 004870 

0.  596 

4.  50 

0.  60 

0. 003420 

0.  599 

5.  00 

0.  60 

0. 002493 

0.  600 

5.  50 

0.  60 

0. 001873 

0.  601 

6.  00 

0.  60 

0. 001443 

0.  601 

6.  50 

0.  60 

0. 001134 

0.  602 

7.  00 

0.  60 

0. 000909 

0.  602 

7.  50 

0.  60 

0. 000739 

0.  603 

8.  00 

0.  60 

0.  000609 

0.  603 

8.  50 

0.  60 

0. 000507 

0.  6  035 

9.  00 

0.  57 

0. 000428 

0.  6035 

9.  50 

0.  .2 

0. 000363 

0.  604 

10.  00 

0.44 

0. 000312 

0.  604 

10.  50 

0.  35 

0. 000269 

0.  604 

11.  00 

0.  24 

0. 000234 

0.  604 

B  - 
o 

0.  311654 

~T? 

Figures  5.42  and  5.43  show  the  .4  E -4  equatorial  omnidirectional  trapped  model 
electron  fluxes  from  0.05  to  4.50  Me\  .  Epoch  1064  represents  solar  minimum  and 
epoch  1967  solar  maximum.  Using  the  above  expressions  with  Table  5.  3  and  Figure 
5.42  and  5.43,  flux  estimates  can  be  made  at  non -equatorial  latitudes. 

Measurements  from  the  0X1-19  satellite  have  indicated  that  the  AE4  model 
fluxes  may  be  significantly  too  low,  particularly  at  higher  energies  beyond  1  MeV.  ' ’ 
These  newer  data  were  averaged  over  periods  that  included  two  magnetic  storms  in 
1969.  The  OV1-19  instrumentation  measured  radiation  belt  electrons  in  the  53  keV 
to  5.  1  MeV  energy  range  in  24  differential  energy  bands,  which  significantly  im¬ 
proved  the  high  energy  data  coverage  over  that  which  had  been  available  to  construct 
the  AE4  models.  The  OV1-19  data  has  now  been  incorporated  into  a  new  NASA 
model  called  AE7-III. 

75.  Vampola,  A.  I..  (1977)  A  New  Study  of  the  Outer  Zone  Electron  Environment, 

A  Hazard  to  CMOS.  SAMSO-TK-77  -  127. 
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Figure  5.44  shows  the  equatorial  electron  flux  versus  L-shell  profile  as  given 
bv  the  AE7-H1  model.  Note,  however,  that  no  direct  equatorial  measurements 
have  been  included  in  the  empirical  model  for  electrons  (E  >  300  keV)  above  8000  km 
and  below  1.  =  5.  Much  of  the  data  base  was  acquired  from  satellites  orbiting  at  a 
significant  inclination  to  the  magnetic  equator,  making  equatorial  flux  representa¬ 
tions  based  on  these  data  uncertain.  Figure  5.45  shows  the  differences  between  the 
AE4  anti  tht  AE7  models  at  energies  above  1  MeV.  The  AE7 -LO  empirical  model 
is  based  on  data  taken  on  the  AZUR  satellite  and  is  shown  in  Figure  5.46.  Note 
that  C  2  1  MeV  electrons  are  most  likely  to  penetrate  spacecraft  shielding  and 
contribute  to  the  accumulated  radiation  dosage  and  damage.  It  is  partially  for  this 
reason  that  the  empirical  model  uncertainties  are  of  interest  to  spacecraft 
designers. 

Magnetic  storms  may  cause  large  energetic  electron  flux  enhancements  that 
last  for  several  weeks  in  the  central  parts  of  the  radiation  belts.  Figure  5.47  shows 
an  examole  of  an  electron  flux  enhancement  as  observed 16  with  the  OGO-5  space¬ 
craft  during  1968.  This  example  shows  that  the  electron  flux  at  1.53  MeV  increased 
by  more  than  4  orders  of  magnitude  during  one  particular  magnetic  storm.  Figure 
5.48  shows  the  mean  exponential  decay  time  of  these  electron  flux  enhancements  in 
days  for  a  wide  range  of  L-shells  and  energies.  Also  plotted  are  the  predicted  elec- 
tron  lifetimes,  (see  Section  5.  5.  1.  2  on  theoretical  electron  models)  which  show 
reasonable,  but  not  perfect,  agreement  with  the  experimental  data  in  this  comparison. 

5,6.  1.4  SHELL  SPLITTING  EFFECTS 

The  outer  belt  particle  pitch  angle  distributions  are  particularly  interesting  due 
to  a  phenomenon  called  shell-splitting.  Here  we  shall  qualitatively  describe  the 
physical  process;  for  analytic  considerations,  see  Roederer.  1  Shell -splitting  arises 
from  the  lack  of  local  time  (azimuthal)  symmetry  of  the  earth’s  magnetic  field  at 
higher  L,  Effects  of  drift  shell -splitting  on  the  trapped  particle  populations  are  considered 
important  for  L  2  4,  but  it  should  also  exist  at  lower  L-shells  where  the  earth's 
magnetic  multipoles  become  significant.  Shell- splitting  arises  because  particles 
that  drift  in  longitude,  preserving  the  first  two  adiabatic  invariants,  modify  their 
pitch  angle  and  radial  location  according  to  the  asymmetric  magnetic  field  topology. 
While  the  concept  of  dipole  L-shell  is  useful  to  describe  principal  features  of  trapped 
particles,  actual  non-dipolar  geometry  with  azimuthal  asymmetry  requires  more 
generalized  coordinates.  One  such  is  the  (strictly  non-invariant)  Mcllwain 


76,  West,  H.  I. ,  Jr.,  Buck,  R.  M. ,  and  Davidson,  G.  T.  (1981)  The  dynamics  of 
energetic  electrons  in  the  earth's  outer  radiation  belt  during  1968  as 
observed  bv  the  Lawrence  Livermore  National  Laboratory's  Spectrometer 
OGO-5,  J,  Geophys.  Res.  :  2111. 
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Figure  5.44.  Empirical  Radiation  Belt  Electron  Fluxes;  (a)  Comparison  of  AKI-7  Model  Spectra  With  a  Number 
of  Data  Set  at  L  =  4.  The  HI  model  curve  is  mainly  based  on  the  OV1-19  observations  from  \  ampola.  (b)  Com¬ 
parisons  of  AEI-7  Model  Spectra  with  a  number  of  data  set  at  L.-6.  6  [from  Yette  et  al^] 


figure  5.  4  7.  Radiation  Belt  Electron  Observations;  An  Intense 
Injection  of  Energetic  (_1.  53  MeV)  Electrons  at  L  =  3  as 
Reported  by  West  et  al ' 6  During  October  1968.  Note  the  four 
orders  of  magnitude  increase  and  the  subsequent  exponential 
decay 

splitting  can  also  result  from  azimuthatly  asymmetric  electric  fields;  for  mathe¬ 
matical  details  see  Schulz  and  Lanzerotti .  5  Figure  5.49  shows  particles  on  the  same 
L -shells  at  local  noon  in  tne  noon-midnight  meridian  plane.  When  radiation  belt  par¬ 
ticles  drift  around  the  earth  to  the  midnight  sector  thev  move  to  a  lower  L-shell  and 
smaller  equatorial  pitch  angles  preserving  their  first  adiabatic  invariant  values. 
Those  particles  starting  closer  to  the  equator  at  noon  drift  to  lower  L-shells  at 
midnight.  Conversely,  Figure  5.50  shows  the  position  of  particles  at  local  noon 
having  initially  been  on  the  same  L -shell  at  local  midnight.  Those  starting  closer 
to  the  equator  at  midnight  move  outward,  closer  to  the  magnetopause  on  the  dav- 
side.  If  they  encounter  the  magnetopause  thev  mav  become  lost  and  there  can  be 
a  preferential  depletion  near  =  90°.  This  gives  rise  to  the  so-called  outer  zone 
butterfly  distribution  which  is  a  pitch  angle  distribution  (BAD)  with  a  minimum 
around  an  equatorial  pitch  angle  <V()  =  77/2.  f  igure  5.51  shows  a  survey  of  the 
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Figure  5.48.  Radiation  Belt 
Electron  Observations;  Pitch 
Angle  Diffusion  Lifetimes  at 
Constant  L  Were  in  Part 
Derived  From  the  Data  in  the 
Previous  Figure,  and  the 
Predicted  Electron  Precipitation 
Lifetimes  are  Those  of  Lyons 
et  al2^  [from  West  et  al‘°] 


energetic  electron  pitch  angle  distributions  in  the  near  equatorial  magnetosphere  as 
7  8 

measured  bv  West.  The  butterfly  distributions  are  clearly  predominant  in  the 
afternoon  sector  after  the  eastward  (counterclockwise)  drifting  electrons  have 
interacted  with  the  magnetopause. 

78.  West,  H.  1.  ,  Jr.  (1979)  The  signatures  of  the  various  regions  of  the  outer 
magnetosphere  in  the  pitch  angle  distributions  of  energetic  particles. 
Quantitative  Modelling  of  Magnetospheric  Processes.  Ed. ,  W.  P.  Olson), 


111 


Figure  5.49.  Theoretically  Computed  Shell  Splitting  Effects  for  Particles  Starting 
on  Common  Field  Lines  in  the  Noon  Meridian.  Dots  represent  the  particles'  mirror 
points.  The  curved  lines  give  tne  position  of  mirror  points  for  constant  equatorial 
pitch  angle  |  Koederer  ] 

Shell  splitting  also  causes  a  coupling  between  pitch  angle  diffusion  and  radial 
diffusion  in  the  outer  parts  of  the  radiation  belts.  Anv  type  of  pitch  angle  diffusion 
may  be  accompanied  by  radial  diffusion  if  the  B-field  is  azimuthally  asymmetric. 

The  direction  of  the  radial  displacement  depends  on  the  longitude  at  which  pitch- 
angle  diffusion  took  place.  Particles  near  the  equator  that  scatter  to  lower  pitch 
angles  on  the  dayside  will  be  radially  displaced  further  from  the  earth  on  the  night- 
side.  Conversely,  displacement  to  lower  pitch  angles  on  the  nightside  leads  to  an 
inward  particle  flux  on  the  dayside.  It  is  estimated  that  particles  spend  2/3  to  3/4 
of  their  drift  period  on  the  effective  dayside  so  that  pitch  angle  diffusion  could  lead 
to  a  net  energy  conserving  outflow  of  particles.  f  irst  and  second  adiabatic 
invariant  conserving  inward  radial  diffusion  as  described  in  the  theoretical  modeling 
section  would  increase  particle  energy  with  inward  radial  motion.  After  undergoing 
many  cycles  of  outward  £  -conserving  diffusion  and  inward  p -conserving  diffusion 
a  significant  local  energization  of  trapped  particles  could  result,"’  but  the  effi¬ 
ciency  of  this  mechanism  needs  to  be  investigated  further. 

79.  Theodoridis,  G.C.  (1963)  Bimodal  diffusion  in  the  earth's  magnetosphere: 

1.  an  acceleration  mechanism  for  trapped  particles,  Ann.  Geophys.  24:944. 
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Figure  5.51.  Survey  of  Energetic  Electron  Pitch  Angje  Distribu¬ 
tions  Observed  in  the  Near-Equatorial  Magnetosphere'6 

5.8.2  Geosynchronous  Altitude  Region 

4 

Geosynchronous  altitude  is  3.6  X  10  km  which  corresponds  to  a  radial  distance 
of  6.6  R„.  A  satellite  at  this  altitude  in  the  plane  of  the  earth's  equator  will  remain 

hi 

fixed  over  the  same  geographical  location.  This  feature  is  highly  useful  for  research, 
communication,  and  surveillance  satellites.  The  natural  geosynchronous  charged  par¬ 
ticle  environment  impacts  the  lifetime  and  reliability  of  satellites  through  radiation 
effects  and  spacecraft  charging.  Many  studies  of  the  geosynchronous  environment 

80  81  g  *,> 

have  been  made  for  example,  Paulikas  and  Blake,  Young,  Garrett,  “  Baker  et 

83  84 

al,  and  Mullen  and  Gussenhoven. 

Near  local  midnight  the  magnetic  field  tines  at  geosynchronous  altitude  often  de¬ 
part  stronglvfrom  dipolar  configuration  during  magnetically  active  periods.  Empiric¬ 
ally  during  such  times  the  pitch  angle  distributions  of  the  particle  fluxes  can  change 
from  being  peaked  in  the  direction  perpendicular  to  the  magnetic  field  line  to  a  more 
isotropic  distribution.  This  and  other  flux  changes  have  been  used  as  diagnostic 

devices  by  Higbie  et  al,  ^  Belian  et  al,^®  Baker  et  al,!il  Baker  et  al,Hl<  and  Belian 
89 

et  al  to  study  underlying  magnetospherie  processes. 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  127.) 
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In  this  section  the  long  term  temporal  behavior  of  energetic  (>  1  Me\  )  electrons 
and  the  plasma  environment  is  emphasized.  Energetic  electrons  penetrate  space¬ 
craft  shielding  and  may  cause  radiation  degradation  of  microelectronic  components. 
The  plasma  environment,  of  which  the  ions  are  an  important  component,  modifies 
the  voltage  to  which  a  satellite  will  charge.  Result.-  the  SCATHA  (Space,  raft 

Charging  at  High  Altitude)  satellite  show  that  the  ion  .  omp  >sitj on  at  geosynchronous 
altitudes  is  a  function  of  magnetic  activity  and  local  time.  *  * 

Energetic  trapped  electron  flux  intensities  I  £  5  1.0  Me\  )  at  geosynchronous 

altitude  have  been  shown  to  be  positively  correlated  with  the  average  solar  wind 
80 

speed.  The  3.9  MeV  integral  electron  flux  (.!<  •:t.PMeV>|,  for  example,  can  varv 
by  about  a  factor  of  5  from  a  solar  wind  speed  of  400  km 'sec  to  one  of  MOO  kn  sec. 


Lower  energy  electron  fluxes  (140  -  600  ke\  ),  on  the  other-  hand  show  little  such 
correlation.  Sufficiently  long  time  averages  t^  1  year)  empirically  remove  the  solar- 
wind  speed  effects  and  reflect  the  overall  average  stability  of  the  electron  fluxes 
over  longer  time  scales.  For  details  see  i’aulikas  and  Blake. 

Energy  spectra  of  energetic  electrons  in  the  geosynchronous  altitude  region 
measured  with  the  SCATHA  spacecraft  are  shown  in  f  igure  5.  52.  These  data 
represent  75-day  averages  derived  from  data  obtained  between  1  ebruarv  1979  and 
February  1980,  and  on  the  average,  the  fluxes  may  be  represented  by  a  power  law 
spectral  dependence.  Integration  of  the  fitted  curves  rises  integral  flux  levels  that 
are  consistent  with  the  AE-4  and  AK-7  models  (Figures  5.42,  5.43,  5.45,  and  5.46). 
This  implies  that  the  long  term  temporal  averages  of  the  electron  fluxes  at  geo¬ 
synchronous  altitude  did  not  materially  change  during  the  1970's.  The  scatter  of 
the  individual  SCATHA  flux  data  measurements  about  the  mean  time -averaged  flux 
is  substantial,  however,  and  at  times  the  observed  electron  fluxes  differed  from  the 
mean  values  by  over  an  order  of  magnitude.  The  flux  models,  therefore,  should 
be  used  with  caution. 
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SC3  AVERAGE  ELECTRON  FLUX 
VERSUS  ENERGY 
POWER  LAW  FIT 
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Figure  5.  52.  Time-averaged 
Energetic  Electron  Energy  Spectra 
Measured  Near  Geosynchronous 
Altitude  From  February  1979  to 
February  1980.  These  curves 
represent  a  75-dav  average  and 
approximate  a  power  law  curve. 
Individual  data  show,  however, 
that  deviations  of  2  orders  of 
magnitude  from  these  means  are 


not  uncommon 
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It  id  useful  to  estimate  the  percentage  of  time  that  the  electron  flux  will  exceed 
a  certain  value.  Figure  5.53  shows  the  cumulative  probability  distribution  for  the 
•SC  ATI  LA  1.4  ~  2.6  AleY  energy  electron  channel  in  four  1.  intervals.  For  example, 
from  this  figure  me  would  expect  to  find  1.  1  -2.6  MeY  electrons  with  a  flux  mten- 


-2  - 1  - 1  - 1 

sitv  greater  than  10  electrons  cm  sec  sr  xe\  50  percent  of  the  time  in  the 


o.o  -6.0  L -shell  region.  On  the  other  hand,  from  the  same  bar  graph,  fluxes  above 

- 1  -l.-i 

xe\  are  expected  '  10  percent  of  the  time. 


100  electron  err.  “  sec  sr 


SC3  HIGH  ENERGY  ELECTRON 
CUMULATIVE  DISTRIBUTIONS 
0419-2603  MV) 


Figure  5.53.  Cumulative  Probability  That  the  High- 
Energy  Electron  Flux  (at  Energies  1.  419  -  2.603  MeV) 
is  Less  Than  the  Levels  Shown.  This  figure  gives  an 
estimate  of  the  "spread"  in  the  individual  measure¬ 
ments  used  to  obtain  f  igure  5.52  [Mullen  and 
Gussenhoven^] 


Ions  are  a  dynamic  component  of  the  radiation  environment  at  geosynchronous 

altitudes.  The  ion  composition  at  low  energies  varies  with  magnetic  activitv  as 

84 

illustrated  in  f  igure  5.54  [Mullen  and  Gussenhoven  ].  During  magnetically  active 
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Figure  5.54.  Average  Oxygen  (Left)  and  Hydrogen 
t  Right)  Number  Densities  as  Determined  From  the 
Lockheed  Ion  Composition  Measurements  on  the 
Sl.’ATHA  Spacecraft  (for  Energies  1  -  32  keV)  vs  Kp 
for  Various  L-Shell  Intervals  (Mullen  and  Gussen- 
hoven84] 


periods  (high  Kp)  the  O  component  becomes  enhanced  relative  to  protons.  An  ex- 

+  Q0  91 

planation  for  O  enhancement  by  Kaye  et  al  and  Fennell  et  al  is  that  during  mag¬ 
netic  disturbances  0+  (and  H+)  ion  fluxes  are  accelerated  up  along  the  magnetic  field 
lines  from  the  auroral  ionosphere,  while  other  proton  (H+)  fluxes  probably  originate 
from  the  magnetotail  plasma  sheet,  which  moves  closer  to  the  earth  during  mag¬ 
netically  active  periods.  Oxygen  ion  flux  enhancements  at  keV  energies  also  increase 
the  trapped  particle  energy  density  relative  to  that  of  the  magnetic  field.  Figure  5.55, 


90.  Kaye,  S.  M. ,  Shelley,  E.G.,  Sharp,  R.  D.  ,  and  Johnson,  R.G.  (1981)  Ion 

composition  of  zipper  events,  J.  Geophys.  Res.  v8gl:3383. 

91.  Fennell,  J.  F.  ,  Croley,  D.  R.  ,  Jr.,  and  Kaye,  S.  M.  (1981)  Low-energy  ion 

pitch  angle  distributions  in  the  outer  magnetosphere:  Ion  zipper  distributions, 
J.  Geophys.  Res.  ^6^:3375. 
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are  operative  during  active  periods  at  geosynchronous  altitude  are  substantially 
different  from  those  operative  during  quiet  times.  Garrett  and  Baker  et  al®'*  give 
empirical  accounts  for  the  geosynchronous  plasma  environment. 


6.7  RING  CURRENT 

The  gradient -curvature  drift  of  radiation  belt  parti cles  causes  differential  motion 
that  is  mass  and  charge  dependent:  electrons  drift  eastward  and  positive  ions  drift 
towards  the  west.  This  constHutes  an  electrical  current  around  the  earth  in  tne 
westward  direction,  called  the  extraterrestrial  ring  current.  During  geomagnetic 
disturbances,  such  as  magnetic  storms,  the  population  of  trapped  particles  at 
~  1  -  800  keV  energies  is  substantially  enhanced  on  1. -shells  between  1.  =  3  and  1.  -  (j. 
As  a  consequence,  the  ring  current  is  intensified  and  magnetic  disturbances  at  the 
earth  result. 

5.7.1  Electrical  Current  Relations 

The  ring  current  itself  produces  a  magnetic  field  that  is  superimposed  on  the 

earth's  magnetic  field.  Enhancement  in  the  ring  current  causes  the  magnetic  field 

depression  observed  at  mid-latitudes  on  the  surface  of  the  earth  during  the  main 

and  recovery  phases  of  geomagnetic  storms  as  well  as  magnetic  field  enhancement 

beyond  L  ~  6  -  7,  as  illustrated  in  Figure  5.  3.  For  this  reason  we  shall  consider 

the  magnetic  effects  of  the  trapped  particles. 

As  one  can  see  from  Eqs.  (5.  36)  and  (5.  37)  the  particle  gradient-curvature 

drift  velocity  is  proportional  to  the  particle  energy.  With  the  equatorial  pitch  angle, 

cv  ,  one  mav  write 
o 


<  1  +  cos^  a  )  (B  >  V  I  B  !  ) 
o 


t a.  10a) 


where  electric  field  and  gravity  effects  have  been  disregarded.  Here  B  is  the  mag¬ 
netic  field  due  to  the  main  (earth's  internal)  dipole  moment  and  M  -  M  (-)  where  (-) 
is  a  unit  vector  in  the  magnetic  northward  direction.  At  the  magnetic  equator 
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For  simplicity  in  the  presentation,  we  shall  consider  equatoriallv  mirroring 

particles  only,  that  is,  0  =  jt/2;  for  the  mathematical  treatment  of  the  general 

°  92. 

case  of  an  arbitrary  pitch  angle  a  ,  see  Dessler  and  Parker.  We  then  get 


3  C  2(p  3  £  ,2(3 

'c  qST  r  =  Cq'ST  L 


(5.  107) 


where  is  a  unit  vector  in  the  eastward  azimuthal  direction  around  the  carta. 
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Maxwell's  equations  (the  Biot-Savart  law)  the  magnetic  field  generated  bv  the 
motion  of  each  particle  is 
2tr 
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drift 
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(5.  103) 


where  i  is  the  magnitude  >f  the  single-particle  drift  current  : 

q  v.l 
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and  thus 
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The  minus  sign  indicates  that  the  particle  azimuthal  drift  generated  field  opposes 
the  main  (internal)  dipole  field  earthward  of  the  ring  current  particle  population. 
There  is  also  a  magnetic  effect  of  the  particle's  spiral  motion  around  the  field 

E 

lines.  Bach  gyro-loop  may  be  considered  a  small  dipole  moment  p  -  j ,  since  for 

a  -  "  2,  £  =  G.  1  .  and  the  associated  magnetic  field  is: 
o  J- 


=  V 


e 

IT 
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which  is  in  tit'.'  direction  of  the  int-rnu!  dipole.  The  total  perturbation  at  the  origin 
due  to  a  single  equatoriallv  mirroring  particle  is  then 


AB  =  B,  a  B 
s  d  g 


2  C 
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I  a. 
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92.  Dessler,  A.  .1.  ,  and  Parker,  B. \.  <  1 '';><•)  ilvdromngnetic  theory  of  geomagnetic 
storms,  J.  Geophys.  H-  s.  37:223  ■. 


On  the  surface  of  the  earth  at  the  equator,  the  unperturbed  imternai)  iipole 
field  is  just  B  -  so  that 


(5.  1  13) 
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and  noting  that  the  total  energy  in  the  earth's  unperturbed  dipole  field  above  the 
earth's  surface  mav  be  written  as 


L.  =  1  2  3  1  M 

m  3  BK  Re  3 

re 


(5.  1 14) 


we  may  express  the  relative  ring  current  single  particle  perturbation  as 


BP  "  3U™  ' 

E  m 


(5.  115) 


it  turns  out  that  this  expression  is  valid  for  trapped  particles  in  the  radiation  belts 
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regardless  of  the  equatorial  pitch  angle  a In  deriving  Eq.  (5.  115)  it  was 
assumed  that  the  total  energy  in  the  ring  current  is  less  than  the  magnetic  field 
energy  U  .  When  that  is  not  the  case  Eq.  (5.  115)  is  no  longer  strictly  valid  and 
may  be  in  error  by  up  to  a  factor  of  2. 

By  summing  up  the  effects  of  all  the  individual  particle  motions  in  the  geomag¬ 
netic  fi°ld,  one  arrives  at  the  total  magnetic  field  perturbation: 

77  ^  Bmax  ^max 

AB  =  /  sin  a  da  f  dL  f  dE  AB  f.(a  ,  L,  C  )  (5.116) 

*7*  J  oof  I  sio 

1  °  1  Cmin 

where  f.(a  ,  L  C  )  is  the  distribution  function  for  particle  species  i,  expressed  as  a 
1  o  » 

function  of  equatorial  pitch  angle,  L-shell  and  energy. 
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Parker  developed  a  hydromagnetic  formalism  that  provides  an  alternative  to 
this  extensive  integration.  One  may  define  the  macroscopic  particle  pressures  in 
the  direction  parallel  and  perpendicular  to  the  local  magnetic  field  direction 


93.  Parker,  E.  N.  '■  1957)  Newtonian  development  of  the  dynamical  properties  of 
ionized  gases  of  low  density,  Phvs.  Bes.  107:924. 
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(5.  117) 


P ..  /  /  F.(X,  v,  cv  )  mv“  cos“  Or  dv  da 

II  J  J  l  o  o  o 

pa  =  i^//1  i(J 


9  2 

■'.(X,  v,  a  )  mv  sin“a  dv  da 
o  o  o 


(5.  1  18) 


where  P  is  the  particle  distribution  function  for  a  particle  species  i  expressed  in 
position,  speed  and  pitch  angle  coordinates.  The  summation  is  extended  over  all 
particle  species.  The  magnetic  field  pressure  is 
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(5.  119) 


With  these  pressure  expressions  the  total  gradient-curvature  drift  current  can  be 
written  as 


!d  = 
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V  P 
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(5.  120) 
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where  c  is  the  velocity  of  light  [see  Williams  ].  The  corresponding  gyration 
current  of  the  particle  distribution  as  a  whole  is 
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(5.  121) 


The  two  terms  within  the  brackets  in  Eq.  (5.  120)  stem  from  the  magnetic  field 
gradient  and  field  curvature  respectively,  and  the  three  terms  within  the  brackets 
in  Eq.  (5.  121)  represent  currents  driven  by  the  particle  pressure  gradient,  the 
magnetic  field  gradient  and  the  magnetic  field  line  curvature.  The  total  current 
of  all  particles  then  reduces  to 


I  =  I,  +  I 
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(5.  122) 
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As  a  rule  of  thumb,  it  requires  a  total  of  4  X  10  ergs  of  particle  kinetic  energy 
to  produce  a  surface  magnetic  field  depression  of  -100  nT  (1  nT  =  ly  =  10  ®  G). 


94.  Williams,  D..1.  ( 1983)  The  earth's  ring  current:  causes,  generation  and  decay. 
Space  Sci.  Rev.  34:223. 

-  Wv 
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5.7.2  Composition  and  Sources 


The  population  of  trapped  particles  that  form  the  bulk  of  the  ring  current  is 
made  up  primarily  of  electrons,  protons,  helium  ions,  carbon  ions,  and  oxygen  ions. 
The  composition  is  found  to  yarv  substantially  with  energy,  location  (such  as  I. -shell) 
and  with  geomagnetic  conditions.  The  heavier  ions,  such  as  He  and  O  ,  may  be 
dominant  during  disturbed  conditions  on  [.-shells  in  the  range  L  =  3  to  L,  =  5.  while 
prolonged  quiet  periods  tend  to  favor  H  (protons)  above  tens  of  keV  energies.  The 
latter  is  also  a  reasonable  expectation  since  at  typical  ring  current  energies 
(-70  ke\  >  the  charge  exchange  lifetime  of  if  is  longer.  At  lower  energies,  below 
a  few  keV,  the  H  lifetimes  are  shorter  than  that  of  He  and  O  lifetimes,  and  the 
opposite  may  be  true.  Figure  5.56  (Reference  1*5)  shows  a  relative  comparison  of 
ring  current  ion  flux  observations  during  four  different  time  periods.  These  results 
pertain  to  G  5  20  keV  energies  while  the  ionic  composition  at  higher  energies  re¬ 
mains  to  be  investigated  observationally. 

Ionospheric  or  atmospheric  ions  probably  form  a  significant  fraction  of  the  ring 
current  population.  This  is  inferred  from  the  observations  reported  by  Shelley  et 
al^  and  Sharp  et  al. '  ^  Low-orbiting  polar  satellites  detected  the  precipitation  of 
oxygen  ions  (L  =  6.8)  during  disturbed  times  and  also  detected  field-aligned  upward 
moving  accelerated  ions  from  the  auroral  ionosphere.  Sufficient  pitch  angle  scatter¬ 
ing  at  higher  altitudes  could  cause  these  upward  moving  ions  to  become  stably  trapped 
and  form  part  of  the  ring  current. 

Based  on  total  ion  (no  mass  resolution)  observations,  it  has  become  clear  that 
the  greatest  contribution  to  the  ring  current  comes  from  -  20-200  ke\  ions  where 

08 

the  mass  composition  is  yet  unknown.  This  is  illustrated  in  Figure  5.57  |\\illinms  ]. 

A  practical  measure  for  the  overall  strength  of  the  extraterrestrial  ring  current 

is  the  D  -index  which  measures  the  middle  latitude  spatially  averaged  decrease  in 
st 

the  horizontal  component  H  of  the  earth’s  surface  magnetic  field  as  reported  by  a  num¬ 
ber  of  magnetic  observatories:  D  =  <AH  •.  Under  this  definition  the  quiet  time  ring 

st 

current  corresponds  to  l)  =  0.  Hourly  values  of  the  1)  ,  index  are  published  by 
r  s  t  s  t 

95.  Lundin,  R. ,  Lyons,  L.  H.  ,  and  Pissarenko,  N.  ( 1980)  Observations  of  the  ring 

current  composition  at  L-values  less  than  4,  Geophvs.  Res.  Letts,  425. 

96.  Shelley,  E.G.,  Johnson,  R.G.,  and  Sharp.  R.  U.  ( 1974)  Morphology  of  energetic 

D+  in  the  magnetosphere,  Magnetospheric  Physics,  (Ed.,  B.  M.  lUcC'ormac) 

D.  Heidel,  Hingham,  Massachusetts. 

97.  Sharp,  R.  C.  ,  Johnson,  H.  G.  .  and  Shelley,  E.G.  ( 1 976)  The  morphology  of 

energetic  0+  ions  during  two  magnetic  storms,  temporal  variations, 

J.  Geophvs.  Res.  JL1 : 3288. 

98.  Williams,  D.  J.  (1981)  An  overview  of  radiation  belt  dynamics.  Proceedings  of 

the  Air  f  orce  Geophysics  Laboratory  Workshop  on  the  Earth's  Radiation 
Belts:  January  21-27,  1981,  Eds.,  R.  G.  Sagalvn,  W.  N.  Spjeldvik  and 
\V.  J.  Burke,  A FG I.- Tit- 9 1 -081  I,  AD  A  113959. 
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Figure  5.56.  Observed  Energy  Densities  of  II  .  and  lie  ,  and  o'*  Ions  in  the 
Radiation  Belts  at  King  Current  Energies  (0.2-  17  ke\  ion'  vs  inbital  I’ar.un- 
eters.  The  data  were  derived  from  the  l'KOMICS  experiment  on  i.ie  I’KOGNd)/.- 
spacecraft.  The  dashed  curves  indicate  apparent  energy  densities  calculated 
from  a  total  ion  (i-J/q)  spectrometer  (at  0.  1  -45  ke\  '  under  the  assumption  that 
onlv  protons  were  measured.  I  he  results  show  the  importance  of  die  heavier 
ions  at  different  times  and  locations.  The  four  panels  represent  lata  for  four 
different  periods*  and  the  Dst-index  historv  is  also  shown  as  a  guirle  to  the  ring 
current  activitv:>0 


NASA  /  National  Space  Science  Data  Center,  Goddard  Space  I  light  t  enter,  Ma.-vl  .in 

Magnetic  storms  generally  have  U  depressions  on  the  order  of  100  to  200  nT  uvi 

large  storms  may  exceed  D  1  =  300  n  T>,  and  the  I)  index  may  .Iso  fluctuate 

substantially  for  other  geomagnetic  conditions  for  which  I)  gener.llv  remains 

s  t 

loss  than  50  nT.  Figure  5.58  shows  an  example  of  the  I)  ^  index  plotted  for  -’aru  - 
December  1072,  and  the  occurrence  of  four  magnetic  storm  periods  in  dune, 
August,  September*  and  October  'November*  is  evident.  The  !)  index  \r>  i  he :  ^  a  > :  ■  • 
ve  rv  useful  to  identify  rr.ajnetii*  storms  from  surfa»  m  a  Lf no  tog  rn  •  re^'mls. 
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Figure  5.  57.  Integral  Representation  of  the  Ring  Current  energy 
Density  as  a  Function  of  Ion  Energy.  The  contribution  to  the  ring 
current  energy  density  is  greater  where  this  curve  is  steepest. 
This  is  also  where  the  ionic  composition  currently  is  unknown  ^ 
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Figure  5.  58.  An  Example  of  the  Time  V  ariation  of  the  Horizontal  Magnetic  Field 
Component  (Dg()  at  the  Equator.  The  large  rapid  drops  in  Ds .  correspond  to  build¬ 
ups  of  the  ring  current  during  magnetic  storms  and  is  followed  b\  its  subsequent 
decay 


5.7.3  Adiabatic  Effects  Produced  by  the  Ring  Current 

Much  of  the  time  the  magnitude  of  the  elect  rival  current  set  up  bv  t  he  a  7  imuthall  v  In 
ing  radiation  belt  particles  changes  slowlv  in  comparison  with  'he  ion  drift  perio  i.  Thor, 
fore,  the  third  adiabatic  invariant,  which  is  proportional  m  'lie  enclosed  magneti  •  flu:-: 
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threading  the  drift  path,  is  must  often  conserved  following  injection  or  ..coele ration. 

iiii 

Siiraas  and  IXivis  have  shown  that  significant  adiabatic  effects  will  take  place 
even  for  moderate  values  of  D  .  To  separate  me  adiabatic  and  nonadiuhatic  fea¬ 
tures  one  can  transform  to  a  D  0  reference.  If  i,(E  I .  '  is  the  equatoriallv 

S  t  'll! 

mirroring  flux  for  E  -  £  j  II.  =  I.  ,  D  t  -  0)  and  j.,(£  1..,)  is  tlie  corresponding 

displaced  particle  flux  for  D  ?  0  then  one  has 

st 

j2tE2>  1.2>-  (E2'E,1  jjlEj.l-,).  15.123) 

This  follows  from  Liouville's  theorem,  which  states  that  the  phase  space  density 

■> 

is  constant  along  dvnamical  particle  trajectories  (f  =  const.  =  j/p“).  for  equatori¬ 
allv  mirroring  non -relativistic  particles,  the  energy  £  is  mapped  from  the 
unperturbed  E  j  energy  bv  conservation  of  the  first  adiabatic  invariant: 

e  =  £  (  (Bg/B  )  (5.  1241 


where  B.;  is  the  post-perturbation  magnetic  field  and  B^  is  the  value  of  the  quiet 
time  magnetic  field  induction. 

For  a  dipole  field  where  =  0.  312  G, 


(5.  1251 


and 


BE 

B2  =  ~j-  +  AB(r)  . 
L2 


(5.  126) 


The  magnetic  flux  enclosed  by  the  drift  path  (the  third  adiabatic  invariant)  is  given 

by 
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-  277  B  2 

<j>2  =  — -j- -  +  2  77  J  A  B(r)  rdr  , 


(5.  127) 


(5.  128) 


99.  Sbraas,  F.  ,  and  Davis,  L.  H.  11968)  Temporal  Variations  of  the  10  keVto 

1700  keV  Trapped  Frotons  Observed  on  Satellite  Explorer  26  During  First 


where  AH(r)  is  the  magnetic  fie  hi  change  induced  bv  the  ring  current  enhancement. 
Equation  (5.  128)  uses  the  fact  that  the  integral  from  0  to  I.  includes  the  return 
magnetic  flux  through  the  earth  so  that  the  net  dipole  magnetic  flux  through  die 
entire  equatorial  plane  is  zero.  This  fact  allows  the  0  to  Lj  integration  interval 
to  be  replaced  bv  one  for  Lj  to*.  The  equatoriallv  mirroring  ions  will  now  be 
located  at  L9,  where  bv  equating  4>j,  and  <J>>: 

L 

BE/L2  =  BE/L1  +  f  ^B(r)rdr  (5.129) 

o 


and 


A  B(r)  =  0.  7  D  .  f(r)  .  (5.  130) 

st 
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The  function  f(r)  is  shown  in  Figure  5.  59  [Soraas  and  Davis).  For  a  given  Dg^, 
can  be  found  and,  hence,  B^.  Knowing  B.>t  the  resulting  energy  £  ^  can  be  deter¬ 
mined.  In  Figure  5.  60  we  show  an  illustrative  example  as  presented  bv  Soraas  and 

„  .  99 

Davis. 


Figure  5.  59.  The  Radial  Dependence 
of  the  Ring  Current  Magnetic  Field 
Used  in  the  Calculation  of  Adiabatic 
Effects  on  Trapped  Protons  by 
Soraas  and  Davis.  ^  Notice  that 
the  decrease  is  greatest  near  L-4 
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Figure  5.  60.  The  Radial  Proton  Integral 
Knergy  Intensity  Profile  as  Measured  on 
Day  109  of  1965  When  Dst  =  46  nT,  Together 
With  the  Transformed  Profile  Correspond¬ 
ing  to  Dst  =  0  Assuming  the  Three  Adiabatic 
Invariants  of  Motion  Conserved.  The  radial 
dependence  of  the  ratio  between  the  mag¬ 
netic  field  after  and  before  the  buildup  of 
the  ring  current  and  the  radial  movement 
AR  of  the  particles  are  shown  in  the  lower 
part  of  the  figure^ 


5.8  RADIATION  EFFECTS  ON  SPACE  SYSTEMS 


Among  the  known  effects  of  particle  radiation  on  space  systems  are  spacecraft 
charging  phenomena  and  effects  of  penetrating  radiation  on  materials.  Specific 
effects  include  detector  malfunction  and  degradation,  optical  svstem  degradation, 
memory  system  alteration,  and  control  svstem  malfunction  or  failure.  1  or  manned 
space  operations,  biological  effects  are  a  major  concern. 

A  crude  measure  for  damage  done  bv  penetrating  energetic  radiation  is  ra  ii  at  jot 
dosage  which  is  measured  in  rads.  This  unit  is  defined  as  an  energy  deposition  of 
100  ergs  (6.25  X  10  MeV)  in  1  g  of  a  material  substance.  This  definition  does  not 
distinguish  between  the  different  kinds  of  incident  radiation  or  the  different  effects 
on  the  material.  Radiation  dosage  is  thus  onlv  an  overall  measure,  and  it  is  dten 
necessary  to  examine  specific  interaction  cross  sections  when  studying  ra  nation 
effects. 


Energy  is  deposited  through  chemical  (molecular  bond  changes,  b  >und  electron 
excitation  and  ionization)  and  nuclear  (element  transmutation,  nuclear  excitation 
and  induced  radioactivity)  interactions.  The  macroscopic  effects  are  evident  in 
device  failure  after  a  critical  level  of  radiation  exposure  is  reached.  Most  often 
this  critical  Level  depends  directly  on  the  nature  and  energy  characteristics  of  the 
incident  radiation. 

A  major  concern  is  the  on-orbit  lifetime  of  microelectronic  devices  that  are 

4  fj 

designed  to  a  specific  level  of  radiation  "hardness  '  (such  as  10  -  10'  rad).  There 

is  in  manv  cases  a  trade -off  between  orbit  choice  and  svstem  lifetime  that  must  be 
determined. 

A  lowest  order  approximation  to  the  expected  radiation  exposure  effects  can  be 
estimated  bv  combining  the  energy  deposition  rate  versus  incident  energy  curves 
of  Janni  and  Radke**^  wdth  the  expected  radiation  belt  flux  intensity  deduced  from 
previous  observations,  as  in  Section  5.6,  or  from  theoretical  modeling.  A  sir:  ole 
(but  very  crude)  approach  to  estimating  the  radiation  dosage  follows; 

It  is  assumed  that  shielding  is  equal  in  all  directions  so  that  a  spherical  shield 
approximation  can  be  used.  The  shielding  is  also  assumed  to  be  aluminum  or  close 
to  it  in  density.  The  incident  omnidi roctional  particle  fluxes  are  normalized  (or 
scaled)  to  unity  at  a  selected  energy  so  that  dosage  need  only  be  calculated  as  a 
function  of  spectral  shape.  The  dose  rate  for  a  given  energy  spectrum  is  found  bv 
multiplying  the  resulting  dosage  by  the  model  (or  measured)  omnidirectional  flux 
at  the  selected  energy.  Total  dosage  is  determined  bv  integrating  over  the  expended 
exposure  time  of  the  satellite. 

Figure  5.61  shows  a  number  of  exponential  spectral  radiation  curves  normalized 
to  unitv  at  1  MeV  for  energetic  electrons  up  to  6  MeV.  By  picking  the  one  curve 
that  most  nearly  approximates  the  actual  expected  radiation  energy  spectrum,  one 
has  a  one -parameter  spectral  representation,  the  spectral  e-folding  energy  g  . 
One  then  proceeds  to  calculate  the  radiation  exposure  for  the  normalized  spectrum 
behind  a  certain  thickness  of  shielding  hv  using  the  curves  in  Figure  5.  62.  This 
process  has  to  be  averaged  over  the  expected  radiation  conditions  for  the  expected 
spacecraft  orbit  during  the  period  of  the  desired  orbital  operations. 


100.  .lanni,  .1..  and  Radke,  G.  (1970)  The  radiation  environment  and  its  effects 
on  spacecraft,  Quantitative  Modeling  of  IUagnetospheric  Processes, 

Ed .  ,  \V.  P.  Olson,  Geophysical  Monograph  21,  AGF. 
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Figure  5.61.  Model  Klectron  Spectra  for 
[Engineering  Applications;  Normalized 
Integral  Omnidirectional  Klectron  Flux  as 
a  Function  of  Klectron  Knergv.  Hv  matching 
model  energv  spectra  to  the  curves  the 
appropriate  K0  can  be  quickiv  found.  The 
flux  at  1  MeY  is  normalized  to  1  electron/cm^ 
[Janni  and  Hadke^^] 


A  similar  technique  also  applies  to  ions,  f  or  protons.  Figure  5.  63  shows  a 
similar  set  of  exponential  spectral  radiation  curves  extending  to  300  MeV,  and 
Figure  b.  64  gives  the  radiation  exposure  dosage  as  function  of  the  aluminum  shield¬ 
ing  thickness.  Notice  that  the  shielding  is  generally  less  effective  in  reducing  the 
radiation  dosage  due  to  the  very  energetic  ions  in  the  radiation  belts. 
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DOSE  VERSUS  SPHERICAL  SHIELD  THICKNESS 
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figure  5.62.  Kadiation  Dosage  Experienced  Behind  Various 
Shielding  Thicknesses  of  Aluminum  Absorber  for  Several  2 

Incident  Electron  Energy  Spectra  Normalized  to  1  Electron/cni 
at  1  RleV.  Multiplication  by  the  1  MeV  flux  gives  the  expected 
radiation  dose  rate  IJanni  and  Kadke100] 

5.8.1  Detector  Malfunctions 

Single  particle  upsets  occur  when  an  incident  single  particle  creates  enough  free 
electrons  in  the  target  material  to  simulate  a  device  logic  state  change.  Single  par¬ 
ticle  effects  are  particularly  severe  in  small  (<  10p.m)  sensitive  regions  of  micro¬ 
electronic  devices.  Figure  5.65  illustrates  the  incidence  of  a  cosmic  ray  (very 
energetic  heavy  ion)  in  a  single  memory  cell  commonly  used  for  onboard  informa¬ 
tion  storage.  Notice  that  the  volume  where  the  ionization  takes  place  is  at  least 
comparable  to  the  sensitive  cell  region  itself.  It  is  presently  not  established  whether 
reduction  in  cell  size  will  always  increase  the  soft  error  or  single  event  upset 
rate;  it  is  conceivable  that  with  very  small  memory  cells  the  ionization  volume 
could  encompass  many  cells. 
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SCHEMATIC  REPRESENTATION  OF  SENSITIVE  REGION 


IN  A  SINGLE  MEMORY  CELL  / 

/ 

COSMIC  RAY 


Figure  5.  65.  Illustration  of  Radiation  Effect;  Electron-Hole 
Generation  Near  a  Sensitive  Region  as  a  Result  of  Local 
IonizatiorL Produced  bv  a  Traversing  Cosmic  Ray  or  Energetic 
Particle  0 


On  the  average,  the  effect  of  single-particle  incidence  on  detector  materials, 

such  as  aluminum  and  silicon,  is  the  generation  of  one  electron-hole  pair  per 

3.6  eV  of  energy  deposition.  Thus,  1  rad  of  incident  radiation  in  1  g  of  material 
1 3 

creates  1.  74  X  10  electron-hole  pairs.  Even  a  moderately  energetic  radiation 
belt  particle  (for  example  1  MeV)  will  create  a  large  number  of  free  charge  carriers 
in  the  detector  material  and  may  lead  to  false  signals.  Prolonged  exposure  to 
energetic  particles  degrades  the  detector  performance  bv  the  accumulation  of  mater¬ 
ial  mierostructural  damage.  For  example,  solid  state  detectors  of  the  Al-Si-Au 
variety  are  found  to  have  a  factor  of  10  increase  in  useful  lifetime  when  the  aluminum 
side  is  facing  the  radiation  exposure  (as  compared  to  the  gold  side).  Very  energetic 
and  very  heavy  cosmic  ray  ions  have  a  particularly  devastating  effect  on  detector 
systems;  for  further  details  see  Adams  and  Partridge**'*  and  McNulty  et  al. 

Figure  5.66  shows  an  example  of  a  nuclear  interaction  occurring  near  a  sensi¬ 
tive  region  of  a  radiation  particle  detector  or  other  solid  state  device.  An  incident 

28 

proton,  for  example,  will  stimulate  a  Si  nucleus  to  emit  an  alpha  particle,  which 

28 

has  a  short  range.  The  recoiling  Si  nucleus  stops  in  even  a  much  shorter  distance. 


101.  Blake,  J.  B.  (1981)  Personal  communication. 

102.  Adams,  J.  H.  ,  Tr.  ,  and  Partridge,  K.  (1982)  Do  Trapped  Heavy  Ions  Cause 

Soft  Upsets  c  t  Spacecraft'*  NRI.  Memorandum  Report  4846. 

103.  McNulty,  P.  J.  ,  Wyatt,  R.  C.  ,  and  Farrell,  G.  E.  ,  Filz,  R.  C.  ,  and 

Rothwell,  P.  L.  (1980)  Proton  upsets  in  LSI  memories  in  space. 

Space  Systems  and  Their  Interactions  With  Earth's  Space  Environment, 

Elds.  H.  B.  Garrett  and  C.  P.  Pike,  Progress  in  Astronautics  and  Aeronautics, 
Vol.  Tl.Publ.  by  American  Institute  of  Aeronautics  and  Astronautics.  ""  ~ 
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Figure  5.66.  Illustration  of  the  Production  of 
Numerous  Secondary  Particles  From  the 
Nuclear  Interaction  of  the  Primary  Particle 
in  a  Sensitive  Volume  of  a  Solid  State  Detector 
Device 


The  combined  effect  can  deposit  tens  of  MeV  in  a  small  (10  X  20  X  20)  p.  m  volume 

103 

element  creating  a  soft'  (data)  error.  Accelerator  data  indicate  that  the  proton 

—  fi  9 

induced  (E  >40  MeV)  soft  error  rate  is  less  than  10  soft  errors/ (protons-cm  ). 

An  upper  limit  to  the  expected  error  rate  (error/sec),  therefore,  can  be  found  by 

_  0 

using  the  proton  flux  models  for  E  >  40  MeV  times  10  .  However,  if  trapped  heavy 

104 

ions  are  sufficiently  abundant  they  could  dominate  the  soft  error  rate. 


5.8.2  Memory  Alteration 

Certain  microcircuitry  used  in  current  spacecraft  instrumentation  has  proven 
very  susceptible  to  the  effects  of  energetic  heavy  ions  in  the  radiation  belts  and  in 
the  cosmic  radiation.  Memory  chips  and  microprocessors  are  frequently  found 
to  have  their  logical  states  and  information  content  severely  altered  by  the  localized 
energy  deposition  process.  .Similar  effects  can  result  from  alpha-particle  emission 
from  nuclear  interactions  and  from  natural  and  induced  radioactivity  in  the  devices 
themselves. 


104.  Adams,  J.H.,  Jr.,  Silverberg,  R.  ,  and  Tsao,  C.  H.  (1981)  Cosmic  Ray 
Effects  on  Microelectronics,  Part  1;  The  Near-earth  Particle  Environ¬ 
ment,  NhL  Memorandum  Report  4506. 
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If  the  effects  are  infrequent  in  occurrence,  engineering  design  emphasizing 
redundancy  of  the  critical  components  could  circumvent  the  problem.  However, 
when  the  effects  are  frequent  and/or  persistent  this  approach  may  not  be  feasible. 

Figure  5.67  shows  the  energy  deposited  in  a  10pm-  thick  sensitive  region  by 
different  ions  over  a  range  of  incident  energies.  The  vertical  scale  on  the  right 
denotes  the  number  of  electrons  produced.  As  a  circuit  becomes  smaller  and  more 
complex,  less  deposited  energy  (charge)  is  needed  to  trigger  errors. 


0.1  1.0  10 

KINETIC  ENERSY  (MeWNUCLEON) 


Figure  5.67.  Energy  Deposited  in  10  pm  of  Silicon  by  Different 
Ions.  The  scale  on  the  right  shows  the  number  of  free  electrons 
released.  The  low  energy  part  of  the  curve  occurs  when  the 
ion  penetrating  range  is  less  than  10  pm  so  that  it  deposits  all 
its  energy  [E.  Petersen,  NRL1^^] 


105.  Petersen,  E.  (1981)  Private  communication. 
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5.8.3  Control  System  Failure 

Radiation  induced  errors  in  eleetronie  rireuitrv  ran  be  partirularlv  damaging 
when  thev  occur  in  rritiral  cireuilrv  surh  as  control  systems  or  in  derision  making 
logic.  While  other  non -rritiral  circuits  mav  continue  to  function  wdth  false  informa¬ 
tion,  control  systems  can  latch-up,  that  is,  be  switched  into  an  undesired  mode 
from  which  there  mav  be  no  reset  option.  Certain  circuitry  switching  may  cause 
burnout  of  electrical  systems  or  even  worse  effects,  particularly  when  propulsion, 
attitude  or  weapons  systems  mav  be  involved.  For  these  reasons  it  is  imperative 
that  proper  safeguards  and  redundancy  design  be  considered  in  the  earlv  stages  of 
spacecraft  engineering. 

5.8.4  Biological  Effects 

There  is  an  extensive  literature  on  space  biology  [see  for  example  Bacq  and 
106 

Alexander  |.  Here  we  shall  only  point  out  that  the  quiet  time  radiation  belts  at 
some  locations  present  a  lethal  radiation  dosage  to  a  man  in  a  space  suit  or  even 
within  a  vehicle.  But  even  outside  the  main  trapped  radiation  zone,  there  are 
intermittent  high  fluxes  of  solar  energetic  particles.  For  example,  it  is  believed 
that  the  energetic  particle  fluxes  associated  with  the  August  1972  solar  flare/mag¬ 
netic  storm  event  would  have  been  extremely  harmful  to  humans  almost  anywhere 
in  the  earth's  space  environment.  The  method  presented  above  can  also  be  used  to 

estimate  human  radiation  exposure  behind  different  shielding  designs.  A  definitely 

107 

lethal  dosage  is  about  500  rads,  although  lesser  amounts  of  radiation  can  be 
harmful  too. 

5.9  MAN’S  IMPACT  ON  THK  RADIATION  BELTS 

The  activity  of  mankind  can,  to  a  significant  degree,  influence  the  earth's  radia¬ 
tion  environment.  Examples  are  nuclear  detonations  (fission  and  fusion),  accelerator 
particle  beams  (neutral  and  charged),  release  of  chemical  substances,  injection  of 
metallic  powders,  and  electromagnetic  wave  energy  production.  The  effects  of  some 
of  these  modification  sources  have  not  vet  been  studied,  but  for  others  a  substantial 
bodv  of  knowledge  is  available. 


106.  Bacq,  Z.  M.  ,  and  Alexander,  P.  (1961)  Fundamentals  of  Radiology. 

Pergamon  Press,  New  York. 

107.  Desrosier,  N.W.,  and  Rosenstock,  H.  M.  (1960)  Radiation  Technology  in  Food, 

Agriculture  and  Biology,  The  Avi  Publishing  Co.  ,  Westport,  Connecticut. 
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5.9.1  Nuclear  Detonations 

A  great  pulse  of  electromagnetic  energy  and  a  vast  number  of  free  neutrons 
and  other  particles  are  released  in  nuclear  explosions.  In  the  nuclear  fission  pro¬ 
cess,  the  fission  fragments  also  carry  significant  kinetic  energy.  The  product  of 
nuclear  fusion  is  generally  a  stable  particle  (He),  which  may  be  ionized.  As  a  rule 
of  thumb,  '  10  fast  neutrons  (each  of  which  decay  into  a  proton-electron-neutrino 
triplet)  are  released  per  megaton  nuclear  explosive  yield. 

The  size  of  the  nuclear  fireball  depends  not  only  on  the  explosive  yield,  but 

also  on  the  medium  in  which  the  detonation  occurs.  In  field-free  empty  space  the 

fireball  will  expand  without  limits,  but  in  the  presence  of  material  substances  or  a 

magnetic  field  the  fireball  is  effectively  restrained.  In  a  dense  gas  (such  as  below 

'  100  km  in  the  earth's  atmosphere,  collisions  between  the  explosion  products  and 

the  atmospheric  constituents  dissipate  much  of  the  detonation  energy  as  heat.  About 

half  of  this  energy  is  radiated  away  and  the  thermalized  remainder  is  typically  at 
i  Ob 

ti 0 00  to  8000  k.  At  an  altitude  of  60  km  in  the  earth's  atmosphere,  a  1  megaton 
fissional  detonation  will  have  a  fireball  radius  of  ~  4  km,  and  for  the  same  nuclear 
explosive  yield  this  radius  will  be  smaller  close  to  the  ground.  The  fireball  itself 
may  accelerate  to  velocities  of  several  km /sec  due  to  buovancy  and  shock  processes. 
High  altitude  and  space  detonations  (more  than  100  km  above  the  earth)  have  the 
fireball  size  limited  by  the  magnetic  field.  This  occurs  because  the  explosion 
generates  electrically  charged  fragments  which  are  thus  susceptible  to  the  magnetic 
force,  q  v  X  B,  where  q  is  the  particle  charge  and  v  its  velocity.  A  nuclear  detona¬ 
tion  of  1  megaton  can  have  fireball  expansion  to  ~  1000  km  across  the  magnetic  field 
when  the  B-field  has  a  value  B  ~  0.  5  G.  The  expansion  is  not  magnetically  limited 
along  the  field  lines,  except  for  the  mirror  force  in  converging  magnetic  field 
topology. 

Depending  on  the  location  of  the  nuclear  detonation,  a  certain  fraction  of  the 
neutrons  may  decay  within  the  earth's  magnetic  field  trapping  region,  and  the  decay 
products  will  thus  constitute  artificially  created  trapped  radiation,  l'or  nuclear 
fission,  the  fission  fragments  also  emit  particles  (such  as  electrons  and  a-particles) 
before  reaching  a  nucleonic  configuration  as  a  stable  isotope.  This  process  further 
contributes  to  the  trapped  radiation,  and  the  characteristic  electron  energy  is  1  to 
8  MeV  from  this  source. 


108.  Zinn,  J.  ,  Hoerline,  H.  ,  and  Petschek,  A.G.  ( 1966)  The  motion  of  bomb 
debris  following  the  Starfish  test.  Radiation  Trapped  in  the  Earth's 
Magnetic  Field,  D.  Reidel,  HollancT 


122 


Studies  of  nuclear  detonation  effects  have  shown  that  even  small  high  altitude 
explosions  (in  the  kiloton  range)  affect  the  radiation  belts  considerably  (see  for 
example  reviews  by  Hess*^  and  Walt.  * 1  ^  Table  5.4  gives  an  overview  of  the 
known  radiation  belt  effects  of  the  Teak,  Orange,  Argus- 1,  Argus-2,  Argus-3,  and 
Starfish  nuclear  detonations  carried  out  at  high  altitudes  by  the  United  States,  and 
the  USSR-1,  USSR-2,  and  USSR-3  high  altitude  nuclear  detonations  by  the  Soviet 
U  nion. 

Energetic  charged  particles  exiting  the  upper  atmosphere  along  the  geomagnetic 
lines  of  force  are  generally  within  the  atmospheric  bounce  loss  cone.  In  the  ab¬ 
sence  of  significant  pitch  angle  scattering,  such  particles  will  follow  the  field  lines 
and  precipitate  into  the  conjugate  hemisphere.  Empirically,  however,  a  significant 
fraction  of  the  nuclear  detonation  particles  become  trapped  in  the  radiation  belts. 
This  implies  that  significant  pitch  angle  scattering  must  take  place  from  the  angular 
source-cone  region  (aQ  <  «oL(J  to  stabiv  trapped  particle  orbits  ( a  >  aQy  c-),  and 
this  pitch  angle  scattering  must  take  place  on  the  time  scale  of  a  single  half-bounce 
T  '2  |whichis  of  the  order  of  seconds]  (see  Figure  5.5). 

Following  an  artificial  injection  of  particles  into  the  radiation  belts,  the  normal 
radiation  belt  radial  and  pitch  angle  diffusion  mechanisms  will  operate.  The  initial 
narrow  injected  radial  distribution  will  broaden,  and  the  charge  exchange  <for  ions), 
and  Coulomb  energy  degradation  mechanisms,  and  transient  plasma  wave  inter¬ 
actions  will  modify  the  characteristics  of  the  injected  distributions.  Depending  on 
the  location,  the  artificial  radiation  belts  may  last  for  days  or  years.  1 '  '  ' 

The  effects  of  accelerator  beams  injected  into  the  radiation  belt  region  are 
likely  to  be  similar  to  those  of  the  nuclear  detonation  particles,  bat  the  yield  (in 
terms  of  number  of  particles)  is  likely  to  be  much  smaller.  On  the  other  hand, 
since  the  beam  particles  may  be  generated  over  a  very  wide  range  of  energies 
(thermal  to  relativistic)  a  more  precise  study  of  their  effects  is  warranted. 


109.  Hess,  W.  X.  (1968)  The  Radiation  Belt  and  .Magnetosphere,  Blaisdell  Pub¬ 

lishing  Company,  Waltham,  Massachusetts. 

110.  Walt,  M.  (1977)  History  of  artificial  radiation  belts.  The  Trapped  Radiation 

Handbook,  (Eds.,  J.  B.  Cladis,  G.  T.  Davidson,  and  F.  I..  Newkirk), 
L.ockheed  Palo  Alto  Research  Laboratory,  DN'A  252411,  Revision, 

January  1977. 

111.  Walt,  M.  ,  and  Newkirk,  L.  L.  (1966)  Addition  to  investigation  of  the  decay  of 

the  Starfish  radiation  belt,  J.  Geophvs.  Res.  7  1 :  1966. 

112.  Stassinoupolous,  E.G.,  and  Verzario,  P.  (1971)  General  formula  for  decay 

li feti mes  of  Starfish  electrons,  J.  Geophvs.  Res.  76 : 1841. 
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6.9.2  Releaae  of  Chemicals 


Chemical  releases  for  research  purposes  have  been  carried  out  at  high  altitudes. 
In  most  cases  barium  or  lithium  was  released  to  trace  magnetic  field  lines  locally 
and  to  assess  the  magnitude  of  electric  fields  and  upper  atmosphere  winds.  Chemi¬ 
cal  releases  into  outer  regions  of  geospace  are  also  planned.  Such  programs  may 
modify  the  environment  locally  (for  example,  by  altering  plasma  wave  dispersion 
characteristics),  but  are  not  expected  to  impact  the  radiation  belts  seriously  unless 
large  quantities  of  chemicals  are  used. 

Extensive  operations  with  rocket  propulsion  or  special  ion  engines  could,  how¬ 
ever,  drastically  alter  the  different  particle  populations  and  could  lead  to  profound 
changes  in  the  radiation  belt  structure.  To  date  no  comprehensive  environmental 
impact  analysis  has  been  carried  out. 

6.9.3  Transmission  of  Radio  Waves 

It  has  been  suggested  that  electromagnetic  wave  energy  from  tropospheric 
thunderstorm  activity  and  whistler-mode  waves  from  VLF  radio  transmitters  can 
perturb  the  energetic  electron  component  of  the  earth's  radiation  belts.  Correlative 

studies  indicate  that  energetic  electron  precipitation  not  only  occurs  from  natural 
113 

sources  but  is  also  associated  with  strong  terrestrial  radio  transmitter  radio 
operations.  117  Precisely  to  what  extent  man's  electromagnetic  wave  generation 
influences  the  overall  radiation  belt  structure  is  not  known,  however. 

6.9.4  Effects  of  Space  Structures 

Proposed  operations  of  large  man-made  metallic  and  electrically  insulated  space 

structures  will  produce  local  "singular"  regions  in  the  magnetosphere.  Associated 

with  space  shuttles,  space  platforms  or  space  power  arrays  will  be  hydromagnetic 

wakes  in  which  the  wave  and  particle  behavior  will  go  through  a  sudden  change.  It 

it  not  known  whether  or  not  these  cavity  phenomena  may  have  a  significant  effect  on 

118 

the  radiation  belts  themselves.  For  some  details  see  Garrett  and  Pike1  and 
references  therein. 

As  of  this  writing  it  is  not  easy  to  predict  to  what  extent  man's  endeavor  in  the 
terrestrial  magnetosphere  will  alter  the  natural  space  plasma  and  radiation  environ¬ 
ment.  As  geospace  becomes  more  heavily  utilized  these  questions  should  be  properly 
addressed. 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  127. ) 
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